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Abstract 

Accretion disk theory was first developed as a theory with the local 
heat balance, where the whole energy produced by a viscous heating was 
emitted to the sides of the disk. One of the most important new invention 
of this theory was the phenomenological treatment of the turbulent vis- 
cosity, known the " alpha" prescription, where the {r<j}) component of the 
stress tensor was approximated by (ap) with a unknown constant a. This 
prescription played the role in the accretion disk theory as well important 
as the mixing-length theory of convection for stellar evolution. Sources of 
turbulence in the accretion disk are discussed, including hydrodynamical 
turbulence, convection and magnetic field role. In parallel to the optically 
thick geometrically thin accretion disk models, a new branch of the opti- 
cally thin accretion disk models was discovered, with a larger thickness for 
the same total luminosity. The choice between these solutions should be 
done of the base of a stability analysis. The ideas underlying the necessity 
to include advection into the accretion disk theory are presented and first 
models with advection are reviewed. The present status of the solution for 
a low-luminous optically thin accretion disk model with advection is dis- 
cussed and the limits for an advection dominated accretion flows (ADAF) 
imposed by the presence of magnetic field are analysed. Related problems 
of mass ejection from accretion disks and jet formation are discussed. 



1 Introduction 

Accretion is served as a source of energy in many astrophysical objects, includ- 
ing different types of binary stars, binary X-ray sources, most probably quasars 
and active galactic nuclei (AGN). While first developement of accretion theory 
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started long time ago (Bondi and Hoyle, 1944; Bondi, 1952), the intensive de- 
velopement of this theory began after discovery of first X ray sources (Giacconi 
et al, 1962) and quasars (Schmidt, 1963). Accretion onto stars, including neu- 
tron stars, terminates at an inner boundary. This may be the stellar surface, or 
the outer boundary of a magnetosphere for strongly magnetized stars. We may 
be sure in this case, that all gravitational energy of the falling matter will be 
transformed into heat and radiated outward. 

The situation is quite different for sources containing black holes, which 
are discovered in some binary X-ray sources in the galaxy, as well as in many 
AGN. Here matter is falling to the horizon, from where no radiation arrives, 
so all luminocity is formed on the way to it. The efficiency of accretion is not 
known from the beginning, contrary to the accretion into a star, and depends 
strongly on such factors, like angular momentum of the falling matter, and 
magnetic field embedded into it. It was first shown by Schwarzman (1971), that 
during spherical accretion of nonmagnetized gas the efficiency may be as small 
as 10~^ for sufficiently low mass fluxes. He showed that presence of magnetic 
field in the accretion flux matter increase the efliciency up to about 10%, and 
account of heating of matter due to magnetic field annihilation in the flux rises 
the efficiency up to about 30% (Bisnovatyi-Kogan, Ruzmaikin, 1974; Meszaros, 
1975). In the case of a thin disk accretion when matter has large angular 
momentum, the efficiency is about 1 /2 of the efficiency of accretion into a star 
with a radius equal to the radius of the last stable orbit. Matter cannot emit 
all the gravitational energy, part of which is absorbed by the black hole. In 
the case of geometrically thick and optically thin accretion discs the situation is 
approaching the case of spherical symmetry, and a presence of a magnetic field 
playes also a critical role. 

Here we consider a devclopcmcnt of the theory of a disk accretion, starting 
from creation of a so called "standard model" , and discuss recent trends, con- 
nected with a presence of advection. The flow of papers (see e.g. references 
in Narayan, Barret, McClintock, 1996: Mcnon, Qiiataert, Narayan, 1999) con- 
sidering advection dominated accretion flow (ADAF) as a solution for many 
astrophysical problems should be treated with some caution, because of its 
vague physical backgroiuid. The suggestions underlying ADAF: ignorance of 
the magnetic field annihilation in heating of accretion plasma flow, and electron 
heating ONLY due to binary collisions with protons (ions) were critically anal- 
ysed in several papers (Bisnovatyi-Kogan and Lovelace, 1997, 2000; Quataert, 
1997, Blackman, 1999), and the above mentioned suggestions have been strongly 
doubted. 

Competition between rapid accumulation of observational data, mainly from 
Hubble Space Telescope, X-ray satellites, and rapid as well construction of the- 
oretilal models for their explanation creates sometimes a sutuation when the 
model dyes during the time of its publication. Some ADAF models give a good 
examples of this situation. One is connected with an explanation of the unusual 
spectra of the galaxy NGC 4258, with a nonmonotoneous curve. It was claimed 
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(Lasota et al., 1996), that such spectrum may be explained only by ADAF 
model. Observations (Herrnstein et al., 1998; Cannizzo et al., 1998) gave data 
in radio, in optical/UV region showing existence of spectral maximum on the 
place where suspected minimum had been explained by ADAF, and posing con- 
straints on ADAF model. Restrictions on ADAF model had been obtained from 
radio observations of elliptical galaxies (Wrobal and Herrnstein, 2000). 

Another example is connected with an attempt to prove the existence of the 
event horison in black holes by ADAF model. The figure 7 from the work of 
Menon, Quataert and Narayan (1997) (fig. 2 from Narayan, Garcia and McClin- 
tock, 1997) was considered as a proof of the existence of the event horizon of 
black holes and at the same time as a triumph of the ADAF model. Regretfully, 
the data on this figures occurcs to be incomplete and a full set of observational 
data altered this picture (Chen et al., 1998). 

It is, of course, dangerous to justify any physical model by astronomical 
observations without a firm physical ground. Here we analyse physical processes 
in an optically thin accretion flow at low accretion rates, connected with a 
presence of a small-scale magnetic field. We show, that accurate account of 
these processes strongly restrict the boundaries of the ADAF solution. Namely, 
the efficiency of the accretion flow cannot become less then about 1/3 of the 
standard accretion disk value. It makes senseless the attempts to connect ADAF 
models with existence of the event horizon of a black hole, or to explain why 
we do not see a large luminosity from supermassive black holes in the nearby 
galactic nuclei. 

It seems now, that explanation of unusual observational features, in par- 
ticular, of " under luminous galactic nuclei" is connected with more complicated 
processes than a simple ADAF solution. Energy losses connected with magneto- 
hydrodynamical acceleration of matter from accretion disk, formation of jets 
may be very important. These problems are considered in the last section of 
the review. 

2 Foundation of accretion disk theory 

2.1 Development of the standard model of the disk accre- 
tion 

Matter falling into a compact object tends to form a disk when its angular 
momentum is sufficiently high. It happens when the matter falling into a black 
hole comes from the neighboring ordinary star companion in the binary, or 
when the matter appears as a result of a tidal disruption of the star which 
trajectory of motion approaches sufficiently close to the black hole, so that 
forces of selfgravity could be overcomed. The first situation is observed in many 
galactic X-ray sources containing a stellar mass black hole (Cherepashchuk, 
1996, 2000). A tidal disruption happens in quasars and active galactic nuclei 
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( AGN) , if the model of supermassive black hole surrounded by a dense stellar 
cluster of Lynden-Bell (1969) is true for these objects. 

The models of the accretion disk structure around a black hole had been in- 
vestigated by Lynden-Bell (1969), Pringle and Rees (1972). The modern "stan- 
dard" theory of the disk accretion was formulated in the papers of Shakura 
(1972), Novikov and Thorne (1973) and Shakura and Sunyaev (1973). It is im- 
portant to note, that all authors of the accretion disk theory from USSR were 
students (N.I. Shakura) or collaborators (I. D. Novikov and R.A. Sunyaev) of aca- 
demician Ya.B.Zeldovich, who was not among the authors, but whose influence 
on them hardly could be overestimated. 

The equations of the standard disk accretion theory were first formulated by 
Shakura (1972); some corrections and generalization to general relativity (GR) 
were done by Novikov and Thorne (1973), see also correction to their equations 
in GR made by Riffert & Herold (1995). The main idea of this theory is to 
describe a geometrically thin non-self-gravitating disk of the mass M^, which is 
much smaller than the mass of the black hole A/, by hydrodynamic equations 
averaged over the disk thickness 2h Fig.l shows the disk geometry. 



2.2 Basic Equations 

Let us start with Navier-Stokes equations in Cartesian coordinates, describing 
the motion of the viscid compressible fluid with variable viscosity rj and bulk 
viscosity tje coefficients (Chapman and Cowling 1939) 
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or in vector notation. 
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Here, the usual tensor definitions are used with 6ik and d the unit tensor and 
gi the gravitational accelleration. Summation over repeated indices is assumed. 
The last term in (|l|) is the divergence of the viscosity stress tensor aik , which 
is related to the symmetric and divergentloss deformation tensor Ttk by the 
relation 
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The bulk viscosity 775 is connected with finite rate of relaxation processes (ion- 
ization, chemical or nuclear reactions, etc). It is important only for rapid pro- 
cesses (e.g., shocks) which are not considered here. Thus we take rjB — 
everywhere below. 

For accretion flows, it is natural to express equation (|l|) in cylindrical coor- 
dinates (r,(j},z). For making transformation to cylindrical coordinates we need 
to write in these coordinates some differendial operators on the scalar a, vector 
a = (flr, a^, a^) and tensor 
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a) gradient of the scalar: 
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f) Laplacian of the scalar: 
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g) Laplacian of the vector: 
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For the transformation of equation (3) to cyhndrical coordinates we have 
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The second term in the brackets on the left-hand side of equation (Q) is a scalar 
product of the vector u and tensor Vu, which is the same in Cartesian and 
curvilinear coordinates. The last term on the right-hand side of equation (2) 
is the divergence of the tensor 2riTik, which is calculated using equations (9) 
and (|l^). As a result we get the components of equation in cylindrical 
coordinates: 
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The continuity equation for the disk flow is 
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2.3 Derivation of Thin disk Equations 

For axisymmetric disk flows the 0— partial derivatives of the different scalar 
quantities vanish so that equations (|l4|)-(|l6|) become 
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Many problems concerning thin accretion disks can be solved using the above 
equations in the equatorial plane or averaged over the thickness of the disk. The 
vertical structure of the disk is determined by the equilibrium of the pressure gra- 
dient and the vertical gravitational force y^, which may include self-gravitation 
for sufficiently massive disks. 

Problems connected with convective and meridonal motion in disks need 
more detailed description with the full set of equations (p^)-([l^). A complete 
description requires including the magnetic force in the Navier-Stokes equations 
and the Maxwell equations for the field evolution. 

In the equatorial plane of a thin disk, we have 
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where mirror symmetry about the equatorial plane is assumed. However, note 
that duz/dz ^ even in this case. Neglecting u^, and all z derivatives, we 
obtain the approximate equations 
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Combining equations (| 
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dt r dr 
and (|2^) we get the radial equation in the form 
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Following Shakura (1972), we write the viscosity coefficient rj using the a— 
prescription, 
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where u^o is a sound velocity in the equatorial plane (z = 0), and zq is half- 
thickness of the disk. Integrating the equations ( p5| ) and (^7|) over the thickness 
of the disk, we get an approximate set of equations describing the thin disk 
behaviour, neglecting viscosity in the radial equation ( p^ ) which is taken in the 
symmetry plane of the disk 
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The average values, surface density E, mass flux over the disk M, and the 
values of angular velocity O and specific angular momentum j are determined 
as follows 

pdz, M = 2'Kr I pUrdz — 2Trr'E,Ur, 

= i = ru^ (33) 

The equations (|30|)-(|3^ ) determine the variables (S, M, j), where (u^, ^) are 
expresssed from ( p3[) with 

Ur, M < (34) 



The solution of the equation of vertical balance (g6|) permits to obtain the 
variables 

/ Qp \ 1/2 

Wzo = I ^-^ I , ^^0, Pa (35) 

as a functions of (r, S), when the equation of state is given in the form P{p) 
(e.g. a polytropic). Here po, Pq are the density and the pressure in the equato- 
rial plane. 

When pressure P and velocity Ur derivatives are neglected in (30), the bound- 
ary conditions consist of the value of the mass flux and Keplerian angular ve- 
locity at infinity; and of a condition on the angular velocity derivative near the 
accreting body (black hole, neutron star or other), determining the total angu- 
lar momentum flux through the disk and zero density and pressure on the outer 
boundary of the disk. 

f CM \ 

M(oo) = Afoo, f2(r^(X3) = f — j , p(±zo) - 

n{n,,) = f^,„, or (r„) = (36) 

where M is the mass of the star-disk body. If all terms remain in the equation 
( ^0|) then an additional boundary condition is necessary. It follows from the 
demand, that the solution must go smoothly through the sonic (singular) point 
in the radial motion of the accretion disk to the black hole, or be adjusted 
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continuously to the conditions at the outer boundary on the stellar equator. 
The stellar angular velocity may be much less than the Keplerian velocity of 
the accretion disk. In that case a boundary layer is established between the disk 
and the star, which needs to be treated separately. 

In the case of general equation of state P{p, T) a new variable the tempera- 
ture T appeares, which is connected with a heat flux F (ergs/cm^/s), radiative 
or convective. Two additional equations of the energy balance over the disk 
and heat transfer to its outer boundary must be added together with boundary 
conditions for obtaining the full solution. These equations will be considered in 
later sections. 

For solution of nonstationary accretion problems initial distributions are 
needed 

S(r), j{r), M{r) (37) 

in addition to the boundary conditions. In most cases the external force is 
represented by the gravity of the star, which may be treated as a point mass 
with a potential 
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For a thin disk we can use the expansion 
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2.4 Structure of Polytropic Accretion Disks 

For a polytropic equation of state 

P = Kp^+^, (41) 

taking the main term in the gravitational force Qz from ( ^o|) we get the solution 
of the equation (^) (Hoshi, 1977) 
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where the density in the equatorial plane po is connected with r and zq by the 
relation 
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Some physical quantities from ( |33| ) -(|35|) can be expressed in terms of r and po 
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Formally an isothermal disk has infinite thickness zq, but the density falls ex- 
ponentially with the characteristic height Zi from (^). Instead of (^)-(^6|) we 
have 
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Restrict ourself by stationary accretion disks with -St = 0. Then we have from 



M = const 

and (|3l|) is integrated, giving (Shakura, 1972) 
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The integration constant jo after multiplication by M gives the total (advective 
plus viscous) flux of angular momentum within the accretion disk. A positive 
value of jo corresponds to a negative total flux (due to a negative M) , which 
means that the central body accretes its angular momentum. When jg is neg- 
ative the central body decreases its total angular momentum, while increasing 
its mass. 

For a thin disk neglecting and V and taking only the main term for 
from (^ ) we get from (^) 
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The solution ( |54| ) is used up to the inner boundary of the accretion disk at 
f = Tin- The integration constant can be scaled by the Keplerian angular 
momentum jKi^in), so that 
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Substituting (|55|) into (|5^) we get with account of 
density in the equatorial plane 
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In the formula for viscosity coefficient ( |29| ) zq must be replaced by zi from (^ 
in the isothermal case, 
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Then instead of (|55|)-(|64|) we get a solution 
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remain the same in all 



The formulas for the values of Zi in 
approximations for isothermal case. 

For accretion into a black hole in presence of a free innner boundary with 
^ = 1 the zero order approximation is not valid in the vicinity of ri„ , where other 
terms in ( ^0[ ) must be taken into account. For accretion onto the slowly rotating 
star with angular velocity smaller than the Keplerian velocity on the equator, 
the drop of the angular velocity from Keplerian in the disk to stellar equatorial 
happens in a thin boundary layer, which must be considered separately with 
proper account of the pressure term in (|30|). 



2.5 Standard Accretion Disk Model 
2.5.1 Equilibrium equations 

The small thickness of the disk in comparison with its radius h <^ r indicate to 
small importance of the pressure gradient VP in comparison with gravity and 
inertia forces. That leads to a simple radial equilibrium equation denoting the 
balance between the last two forces occuring when the angular velocity of the 
disk f2 is equal to the Keplerian one f2^, 

/r<M\ 

^ = ^K=[^) . (68) 
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Note, just before a last stable orbit around a black hole, and of course inside 
it, this suggestion fails, but in the "standard" accretion disk model the relation 
( |68| ) is suggested to be fulfilled all over the disk, with an inner boundary at the 
last stable orbit. 

The equilibrium equation in the vertical z-direction is determined by a bal- 
ance between the gravitational force and pressure gradient 

dP GMz 

For a thin disk this differential equation is substituted by an algebraic one, 
determining the half-thickness of the disk in the form 

. 1 / P^ 



The balance of angular momentum, related to the (f) component of the Euler 
equation has an integral in a stationary case written as 

M{j ~ jin) = -27rr^ 2htr4,, = ^'^^'^ 

Here j — v^r — flr^ is the specific angular momentum, is a component 
of the viscous stress tensor, M > is a mass flux per unit time into a black 
hole, jo is an integration constant having, after multiplication by M, a physical 
sence of difference between viscous and advective flux of the angular momentum, 
when jin itself is equal to the specific angular momentum of matter falling into a 
black hole. In the standard theory the value oi jin is determined separately, from 
physical considerations. For the accretion into a black hole it is suggested, that 
on the last stable orbit the gradient of the angular velocity is zero, corresponding 
to zero viscous momentum flux. In that case 



= ^Kvln. (72) 

corresponding to the Keplerian angular momentum of the matter on the last 
stable orbit. During accretion into a slowly rotating star which angular velocity 
is smaller than a Keplerian velocity on the inner edge of the disk, there is 
a maximum of the angular velocity close to its surface, where viscous flux is 
zero, and there is a boundary layer between this point and stellar surface. In 
that case (|7^) remains to be valid. The situation is different for accretion discs 
around rapidly rotating stars with a critical Keplerian speed on the equator (see 
Section 3). Note, that in the pioneering paper of Shakura (1972) the integration 
constant jin was found as in ([72|), but was taken zero in his subsequent formulae. 
Importance of using jin in the form ( [72| ) was noticed by Novikov and Thorne 
(1973), and became a feature of the standard model. 
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2.5.2 Turbulent Viscosity and Instabilities 

The choice of the viscosity coefficient is the most difficult and speculative prob- 
lem of the accterion disk theory. In the laminar case of microscopic (atomic or 
plasma) viscosity, which is very low, the stationary accretion disk must be very 
massive and very thick, and before its formation the matter is collected by disk 
leading to a small flux inside. It contradicts to observations of X-ray binaries, 
where a considerable matter flux along the accretion disk may be explained only 
when viscosity coefficient is much larger then the microscopic one. In the paper 
of Shakura (1972) it was suggested, that matter in the disk is turbulent, what 
determines a turbulent viscous stress tensor, parametrized by a pressure 

tr4, = -apvl = -aP, (73) 

where Vs is a sound speed in the matter. This simple presentation comes out 
from a relation for a turbulent viscosity coefficient r]t « pvtl with an average 
turbulent velocity vt and mean free path of the turbulent element /. It follows 
from the definition of in (^l|), when we take I « h from ( [70|) 

tr,p = pvthr— fa pvtVs = -apv^, (74) 

where a coefficient a < 1 is connecting the turbulent and sound speeds vt = avg- 
Presentations of trcf, in ( [73| ) and ( [z^ ) are equivalent, and only when the angular 
velocity differs considerably from the Keplerian one the first relation to the 
right in (^4|) is preferable. That does not appear (by definition) in the standard 
theory, but may happen when advective terms are included. 

Developement of a turbulence in the acctretion disk cannot be justified sim- 
ply, because a Keplerian disk is stable in linear approximation to the develope- 
ment of axially symmetric perturbations, conserving the angular momentum, 
and against the Reyleigh- Taylor effect. It was suggested by Ya.B.Zeldovich, 
that in presence of very large Reynolds number Re = ^ the amplitude of 
perturbations at which nonlinear effects become important is very low, so in 
this situation turbulence may develope due to nonlinear instability even when 
the disk is stable in linear approximation. Another source of viscous stresses 
may arise from a magnetic field, but it was suggested by Shakura (1972), that 
magnetic stresses cannot exceed the turbulent ones. 

Magnetic plasma instability as a source of the turbulence in the accretion 
discs has been studied extensively in last years (see review of Balbus and Hawley, 
1998). They used an instability of the uniform magnetic field parallel to the axis 
in differentially rotating disk, discovered by Velikhov (1959). It could be really 
important in absence of any other source of the turbulence, but it is hard to 
belive that there is no radial or azimuthal component of the magnetic field in 
matter flowing into the accretion disk from the companion star. In that case the 
field ampliflcation due to twisting by a differential rotation take place without 
necessity of any kind of instability. Explanation of viscosity in accretion disks. 



16 



based on developement of plasma instabilities, have been considered by Coppi 
(2000). 

It was shown by Bisnovatyi-Kogan and Blinnikov (1976, 1977), that inner 
regions of a highly luminous accretion discs where pressure is dominated by 
radiation, are unstable to vertical convection. Developement of this convection 
produce a turbulence, needed for a high viscosity. Other regions of a standard 
accretion disk should be stable to developement of a vertical convection, so other 
ways of a turbulence exitation are needed there. 

For Keplerian angular velocity the angular momentum per unit mass j = 
10 r^l"^ is groving outside. In this respect it is similar to the viscid flow 
between two rotating cylinders (Taylor experiment), when the inner cylinder is 
at rest. There are however important differences between these two problems. 
There are two boundaries on cylinder surfaces with adhesion conditions, and 
free boundaries in the case of accretion disk. The radial equilibrium in the fluid 
between cylinders is supported by the balance between centrifugal, and pressure 
gradient forces, while the last is substituted by the gravity in the accretion 
disks. Finally, the fluid between cylinders in the laboratory experiments is 
usually strongly subsonic, determining a non-compressible medium, while the 
Keplerian motion in the thin accretion disk is highly supersonic. Nevertheless, 
it is hoped that the fundamental theory checked in the laboratory will answer 
the astrophysical questions. 

Phenomenological analysis of the Taylor experiment, and the onset of tur- 
bulence in the "stable" case of the inner cylinder at rest had been done by 
Zeldovich (1981). There are two characteristic specific energies in the problem. 
One is the energy, necessary for performing the ring exchange with conservation 
of the angular momentum of each ring — 2^ '^'•^^ (Ar)^; here the exchange 
is supposed to happen for rings with small distance Ar between them. This 
enegry is positive in the " stable" , and negative in the case unstable to Reyleigh- 
Taylor effect. Another energy Et produced after merging of these two rings 
with conservation of the total angular momentum Et = \r (^) (Ar)^ is al- 
ways positive. The ratio of these two energies is defined by Zeldovich (1981) as 
Taylor number 



^g/= — =4 ■ (75) 



^ ^ d[(c^r^)2]/d r 
Et r^{duj/drY 

The boundary Ty = separate regions stable and unstable to Reyleigh- Taylor 
effect. In the case of close cylinder radiuses Rout — Rin <C i? = \{Rout + Rin) 
we get with account of ([75|) 



UJout ~ W,;„ ~ CO, Rout ~ Rin ~ R, Rout — Rin = t, -p = 7 , Ty W 8 — . (76) 

dr t R 

This case at i corresponds to the flat Couette flow, which is formally 
stable at all Re numbers (Schlichting, 1964) when is constant {v is the fluid 
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velocity), in contradiction with the experiment. Account of the curvature of the 
curve v{r), even small, is very important here, and gives the instability, similar 
to the case with the free outer boundary, also corresponding to Ty = 0, when the 
instability starts at Re = ^ ~ 2000. Note, that account of a compressibility 
leads to the instability even in the case of the flat Couette flow Glatzel (1989). 
The maximum of instability region corresponds to Mach number Ma « 4.9, 
when unstable mode appeares at Re « 84. 

Taylor experiments with rotating cylinders has shown that at Ty > the 
flow is becoming turbulent, but the Reynolds number Re, corresponding to 
the boundary of stability is increasing with increasing Ty. For the case of 
the Keplerian disk ui ^ r^^/^ we have Ty = formally corresponding to 

= I = 0.222, for which the experimental number of the critical Re « 2 x 
10^ (Zeldovich, 1981). Experimental dependence of the critical value Re = 
f{Ty) is characterized by the Couette-flow critical value i?e(0) = 2000, and 
can be approximated by the dependence Re « lO^Ty^ for Ty > 0.6, > 
0.075. The last dependence was interpreted theoretically by Zeldovich (1981) 
by introduction of the split regime of the turbulence between cylinders. 

In the laminar viscid flow between cylinders the exact solution (Schlichting, 
1964) gives for ujin = the rotational velocity distribution as 



v(r) = —2 -5- r w uJout — [r - Rm) = ^^ont- = ^ontx (77) 



The local values of Re' and Ty' numbers, defined by fluid parameters near the 
inner cylinder, have been introduced by Zeldovich (1981) as 

Re'=P^=x''Re Ty'^8-^xTy. (78) 
77 R 

It was suggested by Zeldovich (1981) that turbulence between rotating cylinders 
is installed, when one of the inequalities is fulfilled 

a. Re f{Ty) or b. Re' = {Ty'). (79) 



If the relation (|79|b) is fulfilled first the global condition of instability is not 
reached, and the turbulence is started in the local region near the inner cylinder, 
and was called as a split regime. Assuming that ([79|b) is fulfilled first we get 



Re' f{Ty') J{y) 
Re^ — = 5— =^2/^-' y = xTy. (80) 

Suppose now, that the function (f>(y) — has a minimum at 
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, /(j/m) /oi\ 

y = ym, <Pm = — 2 — ■ 

Than, taking into account we get the relation Ty > ym, so 

cj>[Tv) > ^(y^), Re = Ty^^{y^) < f{Ty), (82) 

the global criterion of the turbulence is not fulfilled, and the split regime is 
realized. The turbulence is started in the region x < Xm = ym/Ty. So, the 
criterion of the beginning of the split regime of the convection has a quadratic 
dependence on Ty in accordance with the experimental results. 

It is noticed by Zeldovich (1981) that a boundary of the split regime of 
turbulence does not necessary mean the beginning of the turbulence in the 
accretion disk, where the condition (|79|a) is probably needed. Two examples of 
different functions are considered by Zeldovich (1981). For the function 

/(y) = ^fi--) ' (83) 
V yoj 

the absolute instability (|79|a) exist only at Ty < yQ, and the minimum is reached 
at 



y=ho, ^ra='-^^i^. (84) 



f{ym) ^ 27 A 

yl, 4y2- 

From the condition of the plane Couette frow it follows A = 2000, and from 
Taylor experiments, assuming there the split regime turbulence and assuming 
(p2[), we get = 10*, t/o = 1.16. Note, that for the Keplerian accretion disk 



we have Ty = f > 1.16. 

Another function from Zeldovich (1981) is 



f{y)=AeM-). (85) 

yo 



where A — 2000, similar to (^3|), and y^ — 0.3 from the Taylor experiment. Here 
Re = /(16/9) = 2000 exp(l. 78/0.3) w 8 x 10^ for the beginning of the absolute 
instability to the turbulence developement. 

The evaluations above show, that the problem of the hydrodynamic insta- 
bility of the Keplerian accretion disk is far from a full understanding, and there 
are arguments, both experimental and theoretical, supporting the hydrodynamic 
origin of the accretion disk turbulence. The cirtain answer could be obtained 
from the full stability analysis of the Keplerian disk, similar to the one made 
for the viscous accretion tori in cylindrical approximation, made by Kleiber and 
Glatzel (1999) for a constant specific angular momentum distribution. 

With alpha- prescription of viscosity the equation of angular momentum 
conservation is written in the plane of the disk as 
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Mij~J^n)^^7Tr^aPoh. (86) 



When angular velocity is far from Keplerian the relation ( |7l| ) is valid with a 
coefficient of a turbulent viscosity 

2 

V = -^apoVsoh, (87) 
where values with the index "0" denotes the plane of the disk. 

2.5.3 Heat Balance 

In the standard theory a heat balance is local, what means that all heat produced 
by viscosity in the ring between r and r + dr is radiated through the sides of 
disk at the same r. The heat production rate Q+ related to the surface unit of 
the disk is written as 

Heat losses by a disk depend on its optical depth. The first standard disk model 
of Shakura (1972) considered a geometrically thin disk as an optically thick in 
a vertical direction. That implies enegry losses Q- from the disk due to a 
radiative conductivity, after a substitution of the differential eqiation of a heat 
transfer by an algebraic relation 

4 acT" 

Here a is the usual radiation energy-density constant, c is a speed of light, T is 
a temperature in the disk plane, k is a matter opacity, and a surface density 

S = 2ph, (90) 

here and below p, T, P without the index "0" are related to the disk plane. The 
heat balance equation is represented by a relation 

Q+ = Q-, (91) 

A continuity equation in the standard model of the stationary accretion flow is 
used for finding of a radial velocity Vr 

MM 

(92) 



Anrhp 27rrE 



Equations (|68|) , (|70|) , (|86|) , (|90|) , (|9l|), completed by an equation of state P{p,T) 
and relation for the opacity k — k{p,T) represent a full set of equatiions for a 
standard disk model. For power low equations of state of an ideal gas P = Pg = 
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pRT {TZ is a gas constant), or radiation pressure P = = and opacity 
in the form of electron scattering k^, or Karammer's formula Kk, the solution 
of a standard disk accretion theory is obtained analytically (Shakura, 1972; 
Novikov, Thorne, 1973; Shakura, Sunyaev, 1973). Checking the suggestion of 
a large optical thickness confirms a self-consistency of the model. One of the 
shortcoming of the analytical solutions of the standard model lay in the fact, 
that solutions for different regions of the disk with different equation of states 
and opacities are not matched to each other. 

2.5.4 Optically thin solution 

Few years after appearence of the standard model it was found that in addition 
to the opically thick disk solution there is another branch of the solution for 
the disk structure with the same input parameters M, M . a which is also self- 
consistent and has a small optical thickness (Shapiro, Lightman, Eardley, 1976). 
Suggestion of the small optical thickness implies another equation of energy 
losses, determined by a volume emission 

Q-P^qph, (93) 

where due to the Kirchoff law the emissivity of the unit of a volume q is 
connected with a Plankian averaged opacity Kp by an approximate relation 
q ~ acTQKp. Note, that Krammers formulae for opacity are obtained after 
Rosseland averaging of the frequency dependent absorption coefficient. In the 
optically thin limit the pressure is determined by a gas P = Pg. Analytical 
solutions are obtained here as well, from the same equations with volume losses 
and gas pressure. In the optically thin solution the thickness of the disk is larger 
then in the optically thick one, and density is lower. 

While heating by viscosity is determined mainly by ions, and cooling is deter- 
mined by electrons, the rate of the energy exchange between them is important 
for a thermal structure of the disk. The energy balance equations arc written 
separately for ions and electrons. For small accretion rates and lower matter 
density the rate of energy exchange due to binary collisions is so slow, that in 
the thermal balance the ions arc much hotter then the electrons. That also 
implies a high disk thickness and brings the standard accretion theory to the 
border of its applicability. Nevertheless, in the highly turbulent plasma the 
energy exchange between ions and electrons may be strongly enhanced due to 
presence of fluctuating electrical fields, where electrons and ions gain the same 
energy. In such conditions difference of temperatures between ions and elec- 
trons may be negligible. Regretfully, the theory of relaxation in the turbulent 
plasma is not completed, but there are indications to a large exhancement of 
the relaxation in presence of plasma turbulence, in comparison with the binary 
collisions (Quataert, 1997). 
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2.6 Accretion disk structure from equations describing 
continuously optically thin and optically thick disk re- 
gions 

In order to find equations of the disk structure valid in both hniiting cases 
of optically thick and optically thin disk, and smootly describing transition 
between them, use Eddington approximation for obtaining formulae for a heat 
flux and for a radiation pressure (Artemoma et al., 1996). Suppose that disk is 
geometrically thin and has a constant density along z- axis. Defining Sr as the 
energy density of the radiation, F^ad as radiation flux in z-direction, Prad as 
radiation pressure we write momentum equations for averaged Kp in the form 
(Bisnovatyi-Kogan, 1989) 



^ ^ -P»,ar. ( A - l) , (94) 



„ -KepFrad- (95) 

dz 

Consider the case when scattering opacity is much larger then absorption 
opacity Hp, and suggest that heat production rate is proportional to the mass 
density p. Then, neglecting the flux in the radial direction we get 

Frad = 2^PZ: (96) 
^0 

where Fq is the flux from the unit surface of the disk at z = h. Substituting 
(H) into (H) we get 

Using (H) in (H) we get 

dPrad „ ^0 2 /-no^ 

C— = -2k,.—p Z. (98) 

dz Zjq 
Introducing the scattering optical depth 



Kepdz = Kep{h — z) = To — HePZ, 

To = Keph = ^KeSo, (99) 



we rewrite (Bg) in the form 



^dPrad^^F^r^^ (100) 
dr So i^e 
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Solve ( |100[ ) with the foUowing boundary condition 



resulting in 



cPrad — Fq 

In the symmetry plane of the disk at r = tq we have 



2 _ 

3 ^ 2to 



cPrad,0 ^ PqIt: + 



Using (103) in (|97|) we get in the symmetry plane, where T = Tq 



where absorption optical depth 



Kppdz — Kpp{h — z) — TaO — KpPZ, 



Introducing effective optical depth 



1/2 



(101) 



(102) 



(103) 



(104) 



(105) 



(106) 



we get finally the expressions for the vertical energy flux from the disk Fq and 
the radiation pressure in the symmetry plane as 



^0 



_ 2acT(f 
3to 

Prad,0 



3^"^3^j ' 



4 
3to 



3 1+3^^+3% 



(107) 



(108) 



At To ^ ^ 1 we have ( p9[ ) from (107). In the optically thin limit ^ tq ^ 1 
we get 



Fq = acTQTaQ, Prad,0 = -^CLcTQTaQ- 



(109) 



Using Fq instead of and equation of state P = pTZT + Prad,o, the equations 
of accretion disk structure together with equation 



Q+^Po 



(110) 
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Figure 1: The dependences of the optical depth tq on radius, r* = r/r^, for 
the case Mbh = 10* Mq, a = 1.0 and different values of m. The thin solid, 
dot-triple dash, long dashed, heavy solid, short dashed, dotted and dot-dashed 
curves correspond to m = 1.0,3.0,8.0,9.35,10.0,11.0,15.0, respectively. The 
upper curves correspond to the optically thick family, lower curves correspond 
to the optically thin family. 



with (5+ from (p8|), have been solved numerically by Artemova et al. (1996). It 
occures that two solutions, optically thick and optically thin, exist separately 
when luminosity is not very large. Two solutions intersect aX m — mi,, where 

Mc^ , 47rcGM 
m- = 7 , ^Edd = , (111) 

-t^Edd He 

and there is no global solution for accretion disk at m > rrib (see Figure 
It was concluded by Artemova et al (1996), that in order to obtain a global 
physically meaningful solution at m > rhb, account for an advective term in 



(110) is needed. 



2.7 Accretion discs with advection 

Standard model gives somewhat nonphysical behaviour near the inner edge of 
the accretion disk around a black hole. For high mass fluxes when central re- 
gions are radiation-dominated (P w P^, k ~ ^e), the radial dependence follows 
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relations (Shakura, Sunyaev, 1973, Bisnovatyi-Kogan, 1999) 

p ^ j-2 ^ oo, T-r~^/^, (112) 
/i-J->0, E - r3/2 J^-i ^ oo, - J- ^ 0, 

where limits relate to the inner edge of the disk with r = r^^, 

J = = 1- (113) 

3 \ r 

At smaller M, when near the inner edge P ^ Pg, k ^ k^, there are different 
type of singularities 

^ ^ ^21/20^1/5 J. _ ^-3/5^3/5 ^ _ ^-2/5^-3/5 _^ 

This results from the local form of the equation of the thermal balance ( |9l| ) . It 
is clear from physical ground, that when a local heat production due to viscosity 
goes to zero, the heat brought by radial motion of matter along the accretion 
disk becomes more important. In the presence of this advective heating (or 
cooling term, depending on the radial entropy S gradient) written as 

the equation of a heat balance is modified to 



g+-Oadv = 0-. (116) 

In order to describe self-consistently the structure of the accertion disk we should 
also modify the radial disk equilibrium equation, including pressure and inertia 
terms 

r{n^ -n\)^-—-Vr — . (117) 
p dr dr 

Appearence of inertia term leads to transonic radial flow with a singular point. 
Conditions of a continious passing of the solution through a critical point choose 
a unique value of the integration constant jin. 

Simplified solutions with inclusion of the advective terms into the verti- 
cally averaged equations describing accretion disks was obtained by Paczyhski 
& Bisnovatyi-Kogan (1981). This approach with some modifications have been 
used by many researchers to study transonic accretion flows around black holes 
(Muchotrzeb & Paczyhski 1982; Matsumoto et al. 1984; Abramowicz ct al. 
1988; Beloborodov 1998). The importance of transonic nature of accretion 
flows on the disk structure has been emphasized by Hoshi & Shibazaki (1977), 



25 



Liang & Thompson (1980) and Abramowicz & Zurek (1981), and later studied 
in more details by Abramowicz & Kato (1989), see also Kato et al. (1998). 

The problems at finding the numerical solution of advective disk structure 
are connected with the possible non-uniqueness of global solutions at a > 0.01 
and the non-standard behavior of a type of a singular point. It was reported by 
Matsumoto et al. (1984), Abramowicz et al. (1988) that in the case of viscosity 
prescription ( |68| ) the singular point changes its type from a saddle to node when 
increasing a. The presence of the nodal-type singular point leads to creating of 
possibility of multiple solutions as the authors have claimed. 

In the paper of Artemova et al. (2001) it was shown that the mentioned 
problems have been created by several inconsistencies in the preceding studies. 
Some problems are connected with an inaccurate averaging of the equations 
over the disk thickness, another ones appear due to incomplete investigation 
of the singular points (Abramowicz et al. 1988). It was found that in the 
case of viscosity prescription ( |73| ) a set of equations describing the vertically 
averaged advective accretion disks has two singular points, independing on a 
and accretion rate. Calculations of the advective disk structure are performed 
in Paczyhski-Wiita potential (Paczyiiski & Wiita 1980) 

(118) 

where M is the black hole mass and 2rg = 2GM/c^ is the gravitational ra- 
dius. The disk self-gravity is neglected. Note, that the multiplicity of singular 



points in solutions for accretion flows in Paczyhski-Wiita potential (118) was re- 
vealed by Fukue (1987), Chakrabarti & Molteni (1993) and Chakrabarti (1996) 
in somewhat different context. Another type of accretion flows with multiple 
singular points was found by Dullemond & TuroUa (1998) and Turolla & Dulle- 
mond (2000) , who considered accretion disks experienced a phase transition. It 
was shown by Artemova et al. (2001), that at a < 0.01 the inner and outer (with 
respect to the black hole location) singular points are of a saddle type, and only 
one integral curve ('separatrix') which crosses the inner point simultaneously 
crosses the outer one. This separatrix corresponds to the unique global solution 
which is determined by two parameters, a and rh = Mc^ / LEdd, for the given 
black hole mass. At larger a > 0.1 the inner singular point changes its type to a 
node, while the outer point remains of a saddlc-typc. There is still one integral 
curve which goes continuously through both the singular points providing the 
unique global solution. 

In the case of viscosity prescription in the first relation of (Q) 

tr</, = w-T = P^^~r (119) 

dr dr 

it was found that there is only one singular point which is always a saddle, 
and only one physical solution which passes through this point exists. Solu- 
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tions which correspond to both forms of viscosity ( |73| ) and ( |119| ) are very close 
quantitatively at low a limit, a < 0.1. 

A numerical method was developed by Artemova et al. (2001) to solve a 
set of equations describing the vertically averaged advective accretion disks. 
The method is based on the standard relaxation technique and explicitly uses 
conditions at the inner singular point and its vicinity. These conditions had 
been obtained by expanding the solution into power series around the singular 
point. Such a modification of the method allowd to construct solutions which 
smoothly pass the singular points and satisfy the regularity conditions at these 
points with high computer precision in wide range of parameters a and m. 

The equations describing the radial disk structure are written for the mid- 
plane density p, pressure P, radial velocity v and angular velocity Q. They 
consist of the mass conservation equation taken the form, 



M = inrhpv, (120) 

M > 0, and h is the disk half-thickness, which is expressed in terms of the 
isothermal sound speed Cg — \/ Pj p of gas, 

h 



(121) 



The equations of motion in the radial (117) and azimuthal direction 



= 0, 



(122) 



M dj 

At: dr ' dr' 

where fix is the Keplerian angular velocity, flj^ = GM/r{r — 2rg)^, j = Or^ is 
the specific angular momentum and tj.^ is the (r, (/))-component of the viscous 
stress tensor. Other components of the stress tensor are assumed to be negligi- 
blly small. The vertically averaged energy conservation equation is written in 
(115), ( pT^ ) 1^ with the advective term in the form 



^^adv n 

Zwr 



dE pd 

dr dr \p 



(123) 



and 



ht 



dn 



rcfj' 



Q- 



2aT^c 



(124) 



dr ' 3k/9/i ' 

are the viscous dissipation rate and the cooling rate per unit surface, respec- 
tively, T is the midplane temperature. The equation of state for accretion matter 
consisted of a gas-radiation mixture is 



^The vertical averaging in equation (116) have been made differently by different authors 
(compare e.g. Sha k ura fc Sunyaev 1973, and Abramowicz et al. 1988). Our choice of the 
coefficients in (123), (124), may be not the optimal one A p osteriory analysis had shown that 



using the factor 4 instead of 2 in the denominator of (123) would be more consistent choice, 
but this change has a little influence on our numerical results. 
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ci = UT 



2>~ 



where TZ is the gas constant. The specific energy of the mixture is 

2 p 



(125) 



(126) 



Both prescriptions of viscosity ( |73[ ) and (119) have been used in calculations 
with V in the form 



V = -acsh. 



(127) 



Note that in the limit Q, — > both the prescriptions ( |73[ ) and (119) coincide. 
Integration of equation (122) gives 



^htr4, = -^(i - Jin) 



(128) 



where the integration constant j\a has meaning of the specific angular momen- 
tum of accreting matter near the black hole horizon. The value of jin is chosen 
to obtain the global transonic solution with the subsonic part at large radii and 
the supersonic part in vicinity of the black hole horizon. In the case of viscosity 
prescription ( [73| ) the expression (128) results in an algebraical equation, and 
the radial structu re o f the accretion disks is described by two first order differ- 
ential equations (117) and (116). In general formulation, these two equations 
require to fix two parameters as boundary conditions to determine the solution. 



In the case of viscosity prescription (119) expression ( |l2q ) results in additional 
first order differential equation, so we have to fix three parameters as boundary 
conditions. 

We show here two examples of numerical results obtained by Artemova et 
al. (2001). Each calculated model is characterized by the unique value of 
jin. Figures 3 and 4 show the dependences of jm on the accretion rate m = 
Mc^/L^dd for three va lues of a = 0.01, 0.1 and 0.5 in the case of viscosity 
prescription (|73| ) and (119) respectively. At low to < 1 the value of jin is 
independent of to and weakly varies with a. In the low a case, a — 0.01 and 
0.1, the values of jin are close to the minimum value of the Keplerian angular 
momentum, (j_R-)min = 3.6742. At high to > 0.1 the values of jm deviate from 
(jif )min to larger or smaller values depending on a. In the case of a = 0.01 and 
0.1 one can see only minor differences between models with different forms of 
viscosity. But, for large a = 0.5 the difference in values of jin increases. 

Figures 5 and 6 show locations of the inner singular points (rs)in as a function 
of TO for different values of a and both viscosity prescriptions. Similar to the 
case of jin discussed above the models at low to show a weak dependence of 
('"s)in on m. In the low a models (squares and circles in Figures 3-6) the values 
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of {rs)in are close to the location of the black hole last stable orbit at r = 6rg. 
At high m > 0.1 the values of (rs)in are decreasing functions of m in the case 
of low a = 0.01 and 0.1, and non-monotonically behave in the case of a = 0.5 
(triangles in Figures 3-6. 

In the case of viscosity prescription (|7^) the solutions have the outer singular 
point in addition to the inner one. The change of value of /3 = form 1 to 
corresponds to the change of a state from the gas pressure to radiative pressure 
dominated one. The thin disks with /3 ~ 1 are locally stable, whereas parts of 
the disk in which /3 ~ are thermally and viscously unstable (Pringle, Rees, 
& Pacholczyk 1973). At large m > 100 the instability can be suppressed by 
the advection effect (Abramowicz et al. 1988). It was shown by Artemova 
et al. (2001), that significant deviations of jin and {rs)in in Figures 3-6 from 
the constant happens when f3 beguns to deviate from the unity. The analysis 
of the critical points had shown that they are either saddle or nodal-type. In 
Figures 3-6 the saddle-type points are indicated by the solid squares, circles and 
triangles. The nodal-type points are represented by the corresponding empty 
dots in the same figures. In the case of viscosity prescription ( [73|) the inner 
singular points can be saddles or nodes depending on values of a and rh. Note 
that the change of type from a saddle to nodal one does not introduce any 
features in the solutions. The outer singular points are always of a saddle-type. 
In the case of viscosity prescription (|ll9| ) the solutions have only one (inner) 
singular point which is always of a saddle-type. 

The models with m < 16 and low a < 0.1 have values of r-g and jin which 
are very close to the location of the last stable orbit. Tin = 6rg, and value of 
jin = {jK)min assumcd in the standard model (Shakura & Sunyaev 1973). The 
radial structures of the models are also very close to the ones for the standard 
model in the same range of rh and a. Such a good coincidence means that 
the advective terms in equations ( |116| ) and ( |117| ) are negligiblly small in the 
considered models. However, the high a models show quite significant deviation 
from the standard model for all m (see Figures 3-6). 

At high accretion rates, rh > 16, the effect of advection becomes significant, 
and at m = 100 the advective part of the energy flux is ~ 2 times larger than 
the luminosity. 

2.8 Convection in accretion discs 

2.8.1 Vertical structure and convective instability in the radiation- 
dominated disk region 

Let us consider in more details the vertical structure of the accretion disk. Fol- 
lowing Shakura(1972), Shakura and Sunyaev(1973) let us suggest that viscous 
heat production per unit mass in the accretion disk is constant, what is used 
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also in section 1.6. For variable vertical density the equation of vertical heat 
flux is different from (BQ), and has a differential form with Fq = Q+ from (Bq): 



dz 



Eo 47r H j So 



what gives after integration 



(129) 



(130) 



Here E is the mass in the layer (0,z), so that S(/i) — So/2, dS = pdz. The 
equation of a heat conductivity in the optically thick disk in vertical direction 
is written as 



^rad " 



dT'^ 



ac dT 



dz 3. .^S • (1^1) 

We consider here only regions where electron scattering dominates, so k — k^. 
Solving (131) together with (130) we obtain 



ToMl 



■^0^0 



(132) 



At the surface of the disk 



So/2, 



Fa 



acT: 



like in black-body radiation, so we obtain from (132) 

3kSo ' 
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(133) 



(134) 



Assuming high optical depth Sok ^ 1, and using ( 134 ) in ( |132| ) we obtain the 
temperature distribution over the thickness of the disk in the form 



y4 ^ y4 



1 



2S 

'^0 



(135) 



Note, that equation ( p9[ ) follows from (134) at high optical depth, after using 
black-body connection between Tg and Fq. 

Consider now a vertical density dependence in the accretion disk. In the 
radiation-dominated region with P Pr ^ we obtain from the vertical 
equilibrium equation (p9) 



a dT^ 
3yO dz 



GMz 



(136) 



Using (|130D,(|131|),(|136D, we obtain 
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L — — z. 



(137) 



After differentiating, we obtain that in the radiation-dominated disk the density 
does not depend on the vertical coordinate z 



^~ 2r3KFo ■ 

In the constant density layer with falling temperature the entropy 
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(139) 



4aT3 

is falling outside in the direction of the gravitation force. It corresponds to 
convectively unstable situation, what was noticed first by Bisnovatyi-Kogan 
and Blinnikov (1976, 1977). 



For outer region of the accretion disk with Pq ^ Pr, k 



there is no 



analytic solution for a vertical disk structure. We may however to check the 
convective instability in another way. In the gas-pressure dominated region the 
entropy 
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and with P = P„ = TZpT it follows from the vertical equilibrium equation (p9| 



We obtain from (|135|) 
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Using (|l4lD,(pI^) in ( |140D we obtain 
dS 5 dT GMz 
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V r^P ^ 4 T4 



(141) 
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It is easy to prove that entropy is growing with height in the region near the 
equatorial plane z <C /i. Here we have E « pi^z ^ Eg, and it follows from 
(p5l),(|l4|),(|43|) 



dS_ 



(GMTl 



EdE 

^0 



(144) 



Calculating the last expression in brackets, using the accretion disk solution (see 
e.g. Shakura & Sunyaev (1973)), we obtain 
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dS ^.SdS ^ 

^ = 

what means stability relative to convection. Similarily we prove a stability 
relative to convection in the region close to the surface of accretion disk, where 
T, — T — Ts, z ^ h, P = Ps- From the theory of stellar atmospheres 
(see e.g. Bisnovatyi-Kogan (2001)) we have a condition at the photosphere, 
corresponding to the optical depth t = 2/3, in the form P, — where 

GMh , , 

9 = (145) 

is the gravitational acceleration at the disk surface. Using these relations and 



(134) in (|143|) we get 



dS fGM^oh 5r4\ 



7^ V 2r3p^ AT^J S 
fGMh 15\ , /3 15\ , 



\r^PsK 32 J \2 32^ 

So, the accretion disk in a standard model is convectively unstable in the ra- 
diation dominated region, which is the nearest to the black hole, and is stable 
in the adjusting region. The most outer gas-pressure dominated regions with 
Kp ^ Ke are also stable relative to convection until the gas is fully ionized. In the 
colder regions with incomplete ionization a behaviour of the accretion disk be- 
comes more complicated, with a non-unique solution, and convective instability 
(Cannizzo, Ghosh & Wheeler, 1982). 

2.8.2 Structure of convective accretion disk region 

In a standard model with a radiative heat conductivity in vertical direction the 
radiation-dominated region exist only when (Shakura and Sunyaev, 1973) 

Existence of convection changes the structure of the radiation-dominated region, 



and condition (146). Let us estimate a part of the energy flux carried out by 
convection, using a mixing length model of the convection (Bisnovatyi-Kogan 
and Blinnikov, 1976). 

The convective heat flow Qconv (ergs/cm^/s), the convection velocity d, and 
the excess AVT of the temperature gradient above the adiabatic gradient (we 
shall take the mixing length I to be the half-thickness zq of the disk) are given 
by the equations 

Qconv = Cppvizo/2)AVT, (147) 
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where 



= 7^ 
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Hence we find that the excess AVT in the temperature gradient does not exceed 
20% of VT, while the convection velocity is w = 3 • 10* cm/sec for r = lO^g, 
what is close to the velocity Vs of the sound in this region. The heat flow is 
carried mainly by convection (Q « Qconv), because of a high matter density, 
and the radiative flux is ^ 20% of the total flux. Convection is smoothing out 
the entropy in vertical direction, so with a precision ^ 20% we may consider 
the entropy as a constant in the vertical direction. 

Bisnovatyi-Kogan and Blinnikov (1977) had obtained an analytic solution 
for the isentropic convective accretion disk in the radiation-dominated region. 



When Pr 
to 



^ Pg the surface pressure at the photosphere with r = 2/3 is equal 
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The effective temperature at the disk surface is connected with the surface 
radiative pressure in (|150[), and with the surface radiative flux in (p^) as 
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The disk thickness is obtained from comparison of two expressions in (151) as 

h = f 1 - (152) 

Sttc V j J 

In the radiative-dominated region the equation of state at constant entropy, 
using (139), is writetn as 



P = 



4 \ 1/3 



9^/3 = i^,p4/3. 



3 V256a, 

Integration of the equation of vertical equilibrium ( |69D gives 

1/3 _ 1/3 GMz^ 
P — Po 



(153) 



(154) 
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Using the surface boundary condition 
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we determine the equatorial density po from (154). Assuming Pq ^ Ps and 
Po^ Ps, we obtain 
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The equatorial pressure is obtained from the angular momentum equation ( [71] 
and (152) as 
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and from the boundary condition (133) we obtain 
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(156) 
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what justi fies the previous assumption. We used here the relation following from 
and (|l5|) 



^0 = '^h. 



(158) 



Note, that we used here the exact angular momentum equation (71) in the plane 
of the disk, instead of the approximate equation obtained by averaging over the 
thickness, used by Bisnovatyi-Kogan and Blinnikov (1977). In this approach 
the equations (156)-( p"5^ ) are valid always in the radiation-dominated region, 
both radiative and convective. The convection does not change the pressure and 
temperature distribution over the thickness of the disk, and the thickness inself, 
changing only its density distribution. The temperature distribution in the 
radiative disk (135) at constant density, and corresponding pressure distribution 

are identical to their distribution in the convective disk (155). 

The density a nd s urface density in the radiative disk are found from (133), 



( |134D , (|156| ) and (|l58D as 
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In the isentropic convective disk the density distribution is given in (155), the 
equatorial density po, surface density Eq, constant K and entropy S are found 
from (|15|),(|156|) as 
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It follows from (159),(16C), that in the equatorial plane the density of the con- 
vective disk is 4 times larger, than in the radiative one, and in the surface density 
the difference is only 64/35 times. Increase of the convective disk density at the 
same temperature narrows the radiative-dominated region, which may exist now 
at luminosities ~ 2 times larger than in the radiative one (146). 



2.8.3 Numerical simulations of convection in accretion disks 

2-D numerical simulations had confirmed in general the conclusion about the 
isentropic vertical structure of the disk in convectivly unstable regions. Fujita 
and Okuda (1998) had obtained that the excess AVT in the temperature gradi- 
ent does not exceed ~ 1%, what is even closer to the adiabatic structure, than 
followed from estimations of Bisnovatyi-Kogan and Blinnikov (1976). It proba- 
bly happend, because in the computations convective cells generally stretch from 
the disk mid-plane to the disk surface, increasing the efficiency of the convective 
mixing in vertical direction, relative to the model of isotropic convection used in 
estimations. 2-D calculations made by Agol et al. (2000) also had shown a satis- 
factory correspondence with the analytical solution of the isentropic convective 
disk structure of Bisnovatyi-Kogan and Blinnikov (1977). 

2-D numerical simulations of accretion of highly viscous nonradiating gas had 
been performed by Igumenshchev (2000), and Igumenshchev and Abramowicz 
(1999,2000). Ideal gas with an adiabatic index 7 was considered, so that pressure 
was connected with the density and specific internal energy e in the form 



P=(7-l)p£. (161) 

Calculations have been done in a range of viscous parameters a = 0.01 — 1, and 
7 = 4/3, 3/2, 5/3. Convective instabihty or large scale circulation occured at 
smaller a and 7. At 7 = 4/3 all models with a < 0.3 are unstable; at 7 = 3/2 
and 5/3 the instability takes place at a < 0.1, and a < 0.03 respectively. 
The calculations have been performed also with account of a turbulent heat 
conductivity, characterized by non-dimensional Prandtl number 

Pr=-, (162) 
X 

where x = -^/p is the thermo metr ic conductivity, and kinematic turbulent vis- 



cosity coefficient v is given in (127). Stabilization of the flow had been obtained 



at Pr = 1 for all 7 with a > 0.1, and for a = 0.03 at 7 = 5/3. 
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The high- viscosity models at a = 1 form powerful bipolar outflows. The pure 
inflow and large-scale circulation patterns occure in the moderate viscosity mod- 
els a = 0.1 — 0.3. All models with bipolar outflows and pure inflows are steady. 
The models with large scale circulation could be either steady or unsteady, 
depending on values of a and 7. All convective models are unsteady. The self- 
similar solution for accretion flows with bipolar outflows (advection-dominated 
inflows-outflows ADIOs) proposed by Blandford and Begelman (1999) have not 
been conflrmed in the numerical simulations of Igumenshchev and Abramowicz 
(2000). 

3-D numerical simulations of viscous non-radiating accretion flows made 
by Igumenshchev et al. (2000) have confirmed qualitatively and quantita- 
tively those obtained in 2-D simulations. It was shown that convective eddies 
are nearly axisymmetric and transport angular momentum inward, tending to 
smooth out the angular momentum distribution. This is in contrast with a small 
scale turbulence, acting as a viscosity, transporting angular momentum outside, 
and tending to smooth out the angular velosity distribution. 

2.8.4 Hot corona and Cyg X-1 variability 

In the convective disk the flow Qac of acoustic energy in the vertical direction 
Is given by Bierman and Liist (1960) as 

Qac ^ pv^{v/vsf ^ 102i-22erg • sec-1 • cm^, (163) 

and is a quantity of the same order as the total energy flux for characteristic 
parameters of Cyg X-1 (Bisnovatyi-Kogan and Blinnikov, 1976). 

Acoustic waves generated in the convection zone escape into optically thin 
layers, induce variable soft X rays in the photosphere, and also arc responsible 
for variable heating of the corona. Comptonization of the photospheric radia- 
tion by hot electrons under conditions where both temperature and density are 
variable should lead to flux variations in the hard range, hu > 5 kcV. This X- 
ray fluctuation mechanism, involving the emergence of waves into transparent 
layers, evidently is of the same wave nature as the variability of the ultraviolet 
excess in stars experiencing intensive convection, such as T Tauri and UV Ceti. 

The waves escaping into the transparent layers and producing variable ra- 
diation in the photosphere and corona occupy a rather narrow frequency band. 
Physically, the reason for this circumstance is the following. A media with 
a high radiation pressure and a non-uniform distribution of plasma along the 
z coordinate (across the disk) serve as efficient filters, picking out a charac- 
teristic frequency range from the broader spectrum generated by convection 
and turbulence. Waves of low frequency and a wavelength exceeding the scale 
height of the atmosphere will not escape outside but will induce oscillations 
of the coronal atmosphere as a whole. On the other hand, under conditions 
where radiation pressure predominates, high-frequency waves experience a se- 
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vere damping because of radiative friction, and their role in heating the corona 
will be insignificant. 

The propagation of acoustic waves through a medium with strong radiation 
pressure, followed by their escape into the atmosphere had been investigated 
in details by Bisnovatyi-Kogan and Blinnikov (1978a, b; 1979), and by Agol 
and Krolik (1998) for MHD waves. Waves emerging into the transparent layers 
will have a phase and group velocity equal to the velocity of sound in gas, 
Vg = ilPgl pY^"^ = ilT^Ty/^, where 7 =5/3 or 7 = 1 according as scattering 
or absorption predominates, TZ is the gas constant, and the temperature T of 
the equilibrium atmosphere is approximately equal to the temperature Tg of 
the photosphere. The characteristic frequency lOc of waves emerging into the 
atmosphere is given by the expression 



-c=(— j .^l-^jsec- (164) 
For the accretion disk model, the gravitational acceleration g at radius r is 



given by ( |145| ), and characteristic thickness h is given by (152). The quantity 
F/ Fc = kF/ gc represents the ratio of the radiation pressure force to the gravity; 
K is the opacity, including both absorption and scattering. In a spherically 
symmetric star of luminosity L and radius R we have g = GM/R^, F — L/AttR^, 
and F/Fc = L/Lc, where Lc — AttcGM/k is the Eddington limiting luminosity. 



Characteristic frequencies ujc from (164) of fluctuations of Cyg X-1 in the con 



vective accretion-disk model of Bisnovatyi-Kogan and Blinnikov (1977), with a 
black hole mass M = 10 M©, luminosity between 0.1 Lc and 0.3 Lc, correspond- 
ing to the photospheric temperature Ts, lay in the range 5-40 msec. If a corona 
with Tcor ~ 10^ Ts is present, waves whose frequency is w = Wc or even little 
lower will be able to escape. This frequency range is in a good accord with the 
observed time-scales of variability (Bolt, 1977). 

Due to strong radiative damping of the sound waves, only small part of 
acoustic energy flux of order 



Qcor - {Pg/Pr)Qac (165) 

is expended in heating the outer layers with an optical depth t < 1. An- 
other mechanism is acting to heat these layers (Bisnovatyi-Kogan and Blinnikov, 
1977). A particle in a region near the surface of the disk that is transparent 
to radiation will be subject to the influence of radiation from the whole disk, 
and not only to the local radiation pressure gradient. The force of the ra- 
diation pressures will accelerate particles in the transparent region above the 
disk. Calculations of the equations of motion for particles subject to radiative, 
centrifugal, and gravitational forces show that in the region with r < 1 the par- 
ticles (protons and electrons) will acquire vertical velocities corresponding to 
the temperatures (1 — 8) • 10* K of protons for L « O.lLc- Turbulent relaxation 
of the particle motions in the corona will tend to equalize the mean energies 
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of the ions and electrons and to form a quasi-Maxwellian particle velocity dis- 
tribution. Thus the combined action of the two heating mechanisms we have 
described will produce around the accretion disk a hot corona with Te ~ 10^ 
K for L « O.lLc. The existence of an analogous corona has been postulated 
phenomenologically by Price and Liang (1977). 

Let us estimate the density Pcor of the corona and the amount of material 
it contains. Since the gas pressure varies continuously with transition from the 
photosphere to the corona, we readily find that at the base of the corona 

Pcor ^ PsTs/T,or ^ Ps ~ g/cm^ (166) 

where p^ , Tg are the density and temperature in the photosphere of the disk for 
the parameters of the source Cyg X-1 in the region of maximum energy release. 
The surface density of the corona is determined by the condition t^s ~ 1 and is 
~ 2g/cm^, which is more than an order of magnitude lower than the density of 
the opaque disk. 

The existence of a hot corona can easily explain the peculiarities of the Cyg 
X-1 spectrum during its variations (Bisnovatyi-Kogan and Blinnikov, 1976). 
The soft X-rays at hv < 7 keV, which comprise w 70% of the total flux, are 
formed in the photosphere of the opaque disk. Some of the radiation (w 10%) 
passing through the hot corona, experience comptonization and transforms into 
a hard radiation up to hv SkT^ « 200 keV, which amounts up to « 30% of 
the total flux. In Cyg X-1 the total luminosity varies around L = O.lLc- 

The changes in the spectrum as the luminosity of Cyg X-1 varies (Holt et 
al., 1975) exhibit the following characteristic behavior. As the total energy flux 
rises, the radiation in the soft part of the spectrum {hv < 7 keV) increases, but 
in the hard range [hv > 10 keV) it remains almost constant, or perhaps may 
even decrease slightly. We shall assume that the variations in the luminosity 
are associated with variations in the power of accretion. As the mass flow M 



rises in the region with Pr ^ Pg the fraction of the acoustic flow (165) used for 
heating of the corona decreases: 

Qcor^{Pg/Pr)Qac-^ M-'- (167) 

For L sa O.IL,,, when acoustic heating predominates, the rise in M may cause 
some decrease in the heating of the corona and in the amount of hard radiation. 
At the same time the flux in the soft range is determined by the radiation of 
the disk photosphere and increases ~ M. In strong bursts of luminosity, when 
L reaches about 0.3ic, the heating begins to be governed by radiation-pressure 
forces, so the temperature of the corona and thereby also the power of the hard 
radiation should increase along with the rise of M and total energy flux. Similar 
explanation of spectral transitions in Cyg X-1 was considered by Poutanen et 
al. (1997). 

Even if a corona with T ~ 10^ K is present in the disk accretion model, 
radiation at > 200 keV cannot be produced, and more energetic electrons 
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are needed. Radiation with energies Ej up to 5 MeV was oserved by CRGO 
from Cyg X-1 (McConnell et al., 2000). It was noted by Bisnovatyi-Kogan and 
Blinnikov (1976), that explanation of hard X-ray and gamma radiation from 
Cyg X-1 would need a presence of a magnetic field in the accretion disk. Two 
ways to form very energetic electrons had been considered. Both required the 
presence of a magnetic field in the disk. A magnetic field could exist in the disk 
either through twisting of the lines of force by differential rotation (Shakura and 
Sunyaev, 1973), or through infall onto the black hole of magnetized material 
having a small angular momentum (Bisnovatyi-Kogan and Ruzmaikin, 1976). 
In the latter case a poloidal magnetic field would be generated. 

In the binary system containing the source Cyg X-1, some of the material 
flowing from the giant star is dispersed in space near the system. The attraction 
of the black hole will not only produce an accretion disk. A small proportion 
of material having a low angular momentum will fall into the black hole and be 
decelerated in the disk. If this deceleration takes place at radii of (10 — 30)rg and 
if a thin coUisionless shock wave is formed (as is very likely in the presence of an 
azimuthal magnetic field) wherein the kinetic energy is transformed into thermal 
energy and « Tj, then hot electrons with T = lO^'^ — 10^^ K will appear. 
The inverse Compton mechanism of interaction of the disk radiation with these 
electrons can lead to the generation of hard radiation with hi/ w 200 — 2000 keV 
or even higher energies. 

Another mechanism for producing fast particles is analogous to the pulsar 
process. If magnetized matter with low angular momentum falls into the black 
hole (in addition to the disk accretion) , a strong poloidal magnetic field will arise 
(Bisnovatyi-Kogan and Ruzmaikin, 1976). By analogy to pulsars (Goldreich and 
Julian, 1969), rotation will generate an electric field of strength E w —{v/c)B in 
which electrons are accelerated to energies e « R{v/c)Be « 3-10'*[i?/ gauss)] 
Mev where v/c k, 0.1 and R k, 10^ cm is the characteristic scale. In a field 
B w 10^ gauss, such electrons will generate synchrotron radiation with energies 
up to ~ 10^ keV. Just as in pulsars, it would be possible here for e^e^ pairs to 
be formed and to participate in the synchrotron radiation. 

3 Boundary layers between disks and rotating 
Stars 

3.1 Boundary Layers 

3.1.1 Structure of the boundary layer. Fitting with the disk solution 

Inside the boundary layer (BL) variables change considerably over the small 
thickness of the layer iJb <C ri„. Matter has no room to accelerate in radial 
direction, so the radial velocity term in (^0|) is negligible, but the pressure 
term is comparable with the gravitational and centrifugal forces (Regev, 1983; 
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Papaloizou and Stanley, 1986; Regev and Hougerat, 1988). The thickness Hb 
of the BL is smaller then its vertical size zq, which also remains small. The 
adopted inequalities for the boundary layer parameters 

iJb < zo < (168) 

will be confirmed by the results. The radius of the star r» differs from the radius, 
at which ^ = 0, by the very small value Hi,- In the asymptotic consideration 
of the BL we use as an inner boundary for the disk solution 

n-n = n. (169) 

The variable x 

r = r^+6x, — <1 (170) 

is used inside BL instead of r. The inner solution within BL is looked for in the 
region < x < oo, while the outer solution (p5|)-(|60|) is valid in < r < oo. 
According to the method of matched asymptotic expansion (MAE) (see Nayfeh, 
1973) the inner and outer solutions are fitted so that values for the outer solution 
at r = are equal to the corresponding values of the inner solution at a; = oo. 
This condition is valid asymptotically at (5 — s- 0. 

Consider a stationary accretion disk BL where equations (p0|), (p2|) , (|5^) are 
valid with the integration constant for ^ in (^5|). The thickness of the disk 
remains small in the boundary layer, so relations (^ij) -(|4^) for the polytropic 
case and (^8|)-(^l|) for the isothermal case are valid. Accounting only for the 
main terms in the asymptotic expansion inside BL, we get the equations from 
(|0l), (H) (Regev, 1983) 

^^-nl^HbPoil^oo'), (171) 



dx 



dx 2TiYiVi,r1 

where 



(6-^), (172) 



n 2 GM 



^ = 7^, ^l*^^- (173) 



The equation (171) was obtained from the exact equation of radial equilibrium 
(^0|), instead of the approximate equation averaged over the disk thickness, used 
originally by Regev (1983). In the MEA approximation these two equations are 
identical (Shorokhov, 2000), in the case of an advective disk the difference is 
essential, and the exact equation (^) should be used. The viscosity coeffisient 
77 is expressed in BL through the coefficient of kinematical viscosity i/f, 

V = Pon- (174) 
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While the radial extend of the BL Hh is much less then its vertical size zo the 
formula cannot be used for a viscosity coefficient approximation. Inside BL 
we thus use the a approximation in the form 

V], = abUsoHb. (175) 

3.1.2 Polytropic case 

This case was investigated by Bisnovatyi-Kogan (1994). It is convenient to use 
the variables (S,w) in the equations ( |l7l| ) ,( |l72| ) . We get from (||) -(|4|) 

PQ = dpo^^, Po = dpoS^^, (176) 
Mso=4oS^, (177) 
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Write for convenience the dimension of d^o 



[dsa] = [v][pr] 



(180) 



From the matching conditions in MAE we need to get Keplerian angular velocity 
r2 = ^K*i uj — 1 from the inner solution at a: — > oo in order to fit the outer 
solution ( |5^ ) at the inner boundary r = . This implies i^b = 1 for the constant 
in (|172|) . We get after transition to the variable E in (171), (172) 



dx 



dpo 2n + 1 „2 



dio 
dx 



dpo 2(n+ 1) 
M 



2tt abvldso 



Dividing (181) on (182) we get 



„ dYj , . - + l 

Db— = -(l + cj)E — 
du! 



(181) 
(182) 

(183) 



with 
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Dh = - 



2(n + l) d 



PO 



M 



2n + l dpodso 2-KabHb rlVt]^^ 



(184) 



The solution of (183) which fit the boundary condition on the steUar surface 
E = aX uj = Lo.^,. Taking into account, that the surface density rapidly grows 
inter the star we may put with sufficient accuracy = cxd and obtain the 



solution, using (179), (|178D and (184) in the form 



where 



K* 



M 
atHb 



(w - CJ,) 1 + 



(185) 



ds — Vn+Tr 



-{2n) 



r(r 



(186) 



r(n + i) 

For fitting the inner BL solution and the outer solution for the accretion disk 
we must make the surface densities (185) at a; = 1 and ( ^8| ) at r = ri„ equal. 
This fitting uniquely determines the outer integration constant ^. After some 
algebraic calculations with account of (|62|),(186), we get 



_2ri + 3 / J. 



2n+l 



where 



[(l-^.)(3 + ^.)]" 



(187) 



471 + 3 _ 2n + 9 4n + 3 , ; 

2-^^{2t:) ~n — {n + 1): 



T(n- 



r(7i+ 1) 



(188) 



The value of Hb is still not determined. It must be found from equation (182), 
where S is substituted from the solution ( |185| ). Before doing this we estimate 
by order of magnitude the values of ^ from (^3|) and ^ from ( |5^ ) 



UK = ilKr, ^ — , ^a{ — 

UK r UK ^ r 




(189) 



The equation (182) contains a nonphysical logarithmic divergency because of 
using the MAE method. For an approximate estimation of the value of iJf, we 
use a characteristic thickness over which the lo variation happens, and for this, 
using the definition (17C) , the following relation may be written 



-r- + = l 

M 



(190) 
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With the help of (p79|),(pi5|),(pi6|) we get from (p^) 



dff = 2(27r) _^ 1) — 



[(l-c..)(3 + c..)]-\ 

r(n + §)- 



r(n + 1) 



Using (|191D,(|192|) in (|187|) we finally get the expression for (1 - ^) 



2n+l 



= d„d^ = 2V7r(27r) 4(„+'i) „4(„+i) (^_|_ i) 2(„+i) 



r(n + i) 



(191) 
(192) 

(193) 
(194) 



Estimating the values by order of magnitude, with account of ( |189| ) we get, 
combining ( 191 ) and ( |193|) 



Qffc r 



1 /£0 
1 — V r 



(195) 



(196) 



Thus we get a complete analytical solution for the accretion disk structure with 
the boundary layer near the star for a viscosity inside the layer from ( |175| ) and 
a polytropic equation of state everywhere. Analytical solutions for an accretion 
disk with boundary layer for some simplifying assumptions about the viscocity 
coefficient and temperature distribution were obtained directly without expan- 
sions by Colpi at.al.(1991) and Glatzel (1992). The MAE method, used by 
Regev (1983) for the solution of the accretion disk BL problem permits to ob- 
tain a self-consistent analytical solution for a disk with boundary layer for a 
viscosity coefficient in the usual a representation. 



3.1.3 Isothermal case 



For an isothermal disk with the same viscosity (175) we get, using 
( pTll) and (|l7|) 



dpQ 
dpo 



1 

duj 
dx 



dx 



"'^*i/,I](l-c.2), 



K 

M(l- w) 



2t: abvK^rl 



(197) 
(198) 
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Dividing the last two equations, and using the approximate boundary condition 
E* = oo, we get the solution (Shakura and Sunyaev, 1988) 



E= - 



MVK 
2TTai,Hp 



(w-w*) 1 + 



(199) 



From (190),(199),(198) we obtain the characteristic scale, which we identify with 



2K 



[(l-c^.)(3 + c.,)]- 



(200) 



Matching ( |66| ) and (199) with account of (20C) like for plytrope, we get the 
expression for the integration constant 



1-e 



/2K 



(201) 



Comparing ( |200| ),( pOlD with ( |195| ), (|196| ) we see that the order of magnitude 
estimates for the polytropic case become exact the for isothermal case apart a 
numerical coefficients close to unity. 



In the isothermal case there is a simple solution of equation (|l98| ) when (199) 
is substituted for E, 



(l-o)) 2(a)-cj*)i+-. (2 + w + tj*) 1+- 



exp 



K 



K 



,x{l — a;*)(3 + uj^,) 



(202) 

We can see from (202) that the scale (200) is naturally appears under the expo- 
nent. It is clear that this solution makes physical sence only over this character- 
istic scale (thickness of BL) where Q < x < r^. It is proved (Nayfet, 1973), that 
the region of applicability of the inner and outer solutions, obtained by MAE 
method, do overlap and the formula, constructed from the inner (i) and outer 
(e) solutions, gives a good interpolation also for the intermediate region. The 
formula, describing the function / in the whole region, has a structure 



/ fi~^fe {fi)e 



(203) 



Here {fi)e is the value of the inner solution on the outer edge [equal to the 
value of the outer solution on the inner edge {fe)i]- The solution for E(r) is 
constructed using E^°'(r) from the relation (^) for the outer solution, and the 
solution of the equation ( p^ with account of (p70|) , (p^) , (p^ for S,(r). The 
choice of ^ from ( |193| ) gives the equality (Ei)e=(I]e)i. The solution for ri(r) is 
constructed using i^Kir) from (M) for the external, and the solution of (182), 
with account of (17C), ( |173| ), ( |185D , for the internal solution ^i{r). It is clear 
from the above consideration that {Qi)i, = {^e)i—^K*- 
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The existence of the overlapping regions, where both solutions (external and 
internal) are vaHd was first proved by Prandtl (1905, cited by Nayfet, 1973) 
for the problem of the boundary layer near the body in the flow of the rapidly 
moving viscid fluid, for which he had invented the MAE method. The interpo- 
lation formulae of the type ( |203| ) had been introduced by Vasil'eva (1959, cited 
by Nayfet, 1973), and is considered as a uniformly valid approximation in the 
whole region, due to existence of the overlapping regions. The direct proof of 
the validity of ( ^03| ) , for the cases where general analytic solutions exist, can be 
found in the books of Nayfet (1981) and Zwillinger (1992). These books, as well 
as Nayfet (1973), contain examples of the application of MAE method, including 
quasi- linear equations of the type (|2|), (|5^), and also references to many other 
applications of the MAE method and its detailed mathematical investigation. 

Using the relation M = 27rr* YiUr for the estimation of the radial velocity Ur 



in the bou n dary layer, we get, from ( |l76[ )-( |l78| ), ( |185[ ), ( |l9l|) for polytropic [or 



from (|199|), (200) for isothermal] cases, the estimation 



Ur/uso ^ ah. (204) 

This means that the solution obtained, where radial velocity was neglected, is 
valid only for sufficiently small viscosity, with aj, -^1. 

3.1.4 Boundary layers with account of thermal processes. Heat pro- 
duction in the boundary layer 

Rapid braking of rotational velocity in the boundary layer during accretion into 
a slowly rotating star leads to strong heating and large energy release. It is 
expected, that matter in the boundary layer is hotter than in the disk, and is 
emitting radiation with a larger effective temperature. The equations of a heat 
balance in the boundary layer are different from the similar equations in the 



accretion disk ( |115| ), (116), (124). The component of the stress tensor should 
be written on the form (119) only, because in the boundary layer D, deviates 
strongly from the keplerian value. The turbulent viscosty coefficient should be 
scaled by a smaller value of boundary layer thickness iJf,, instead of the disk 



thickness h, like in (174), (175). Due to thinness of the boundary layer the radial 



heat flux inside it should be taken into account, and is usually more important 



than the losses through the sides of the disk (5_ in ( |124| ). From ( |124D , ( |119D , 
and using the momentum conservation equation ( |128| ), we obtain the expression 
for the viscous heat production rate in the half-thickness of the boundary layer 
per unit disk surface in the form 

dn , / dnV GMM , ^dto 
Q.-hU,r-^h,(r-J -^(--1)^. (205) 

Here, like in the previous section, we take in the boundary layer = 1, r = r*, 
to = fl/Clx*- It is convenient in the boundary layer to define the value 
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=47rrg+«4^ng+ = ^^(cc>-l)^, (206) 

dr 

determining the heat production by the viscosity in the whole ring of the bound- 
ary layer per unit of its radial thickness in the accretion disk. Let us consider the 
radiative heat conductivity along the radius as the main term in the heat bal- 
ance of the boundary layer. Defining the radial heat flux through the optically 
thick boundary layer $ as 

4acT3 dT ^ , 

$ = inrh, (207) 

Sup dr 

we write the heat balance equation as 

^ GMM, ^^duj 

Integrating this equation, we obtain 



We consider here the layer at w = UJA^ where the temperature has a maximum, 
so the radial flux is zero in this point, and had opposite direction at both sides 
of this layer. A position of this layer should be found from time-dependent 
evolutionary calculation of the boundary layer structure. If, for example, ac- 
cretion is started into the cold neutron star, that the heating wave is started 
from the boundary layer, outside in the accretion disk, and inside the neutron 
star, heating it. In this case the level with lo — uja is moving with time to the 
boundary of the star, and for larger times when the whole star becomes almost 
isothermal, this level coincides with the star boundary. If accretion is going on 
into the hot neutron star, the heat production in the boundary layer may be 
comparable with the heat flux from the star, which is equal to on the star 
surface with oj ~ lo^,. Than we obtain instead of (^09[) , the solution 



GMM, , bj bj^ 

$ = <I>. H 



-(a.-c.,)x(l-|-^). (210) 



It follows from (210), taking into account that the angular velocity is keplerian 
at the outer boundary of the boundary layer Wout = 1, that the total energy flux 
from the boundary layer due to heating by viscosity Qfez is equal to (Popham 
and Narayan, 1995) 

Q,, = ^(l_c..)^ (211) 
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The relation may be obtained also from the conservation laws. When a 

mass dm crosses the boundary layer the change of its energy is 



(212) 



From the angular momentum conservation law the angular momentum lost by 
the accreting matter 



dJm — dmilfc*r'^{i — IjJ* 
is absorbed by the star, so that 

dJ.^, = I^^^lx^dij^, = dJm 



(213) 



(214) 

Obtaining from two last equations I^^dto^ — dmr1(l — w*), we obtain the kinetic 
energy absorbed by the star in the form 



dE^ = I^ilj^^ui^duj^ — dm^.'^^rloj^^l — a>*). 
The heat produced in the boundary layer is defined as 

dEbi = dE„i - dE^ = ^dmilj^^r'^{l - cj*)^, 



what gives (211) with account of 

dEu 



Qb 



dt 



dt 



GM 



K* 



(215) 



(216) 



(217) 



3.1.5 Explicit models of the accretion disk with a boundary layer 

Set of equations, describing stationary accretion disk together with a bound- 
ary layer was solved numerically by several authors. The equations of radial 
and azimuthal motion, and continuity equation coincide with the corresponding 
equations of the advective disk structure (117), (|7l|) with - derivative descrip- 



tion of the viscosity, and a continuity equation (120). Contrary to the advective 
disk structure at accretion into a black hole with a transonic radial motion, 
the accretion disk around a star with a boundary layer has always a subsonic 
radial velocity (204), with a Mach number ~ at, where ab is characterizing the 
viscosity in the boundary layer (174), (175). The viscosity coefhcient describing 
continuously the disk and boundary layer should include a transition from the 
thickness scale h in the disk to the smaller value of the boundary layer thick- 
ness Hb in the layer itself. The expression having these properties was used by 
Popham and Narayan (1995) 



1 



+ 



{2dPo/dr) 



{ahf {iabPof 



'1/2 



(218) 
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The coefBcient (|) is included 2 times into (218) to make a continuous transition 
to the (127) in the disk where h <$:. Ht, — -jp-jj-, and to (175) in the boundary 



dPo /dr ' 

layer where h 3> Hi,. The equation of a thermal balance should include the 
radiative heat transfer $ in radial direction from ( 207 ), which exceeds all othe r 
cooling sources inside the boundary layer. Taking (5+ and from (124), 
advective heating Qadv from (115), and $ from (207), we obtain the equation 
of the heat balance valid in the whole disk, including boundary layer 



Q+ = Q- 



47rr ar 



which may be written as 



ar dr 



inph 



-inr. 



(219) 



(220) 



The thickness of the disk h is everywhere determined by vertical equilibrium 
equation (|70|). Popham and Narayan (1995) used a bit more complicated equa- 
tions of a heat transfer. Instead of Q_ and $ from above they considered 
expressions following from simplified solution of the radiative transfer equation 
permitting to distinguish between absorption and scattering opacities, similar 
to the description in Section 1.6, which is valid also in the optically thin case. 
The equation of state P = Pg+Pr takes into account deviations of the radiation 



pressure Pr from the equilibrium value, similar to (108). 

The boundary conditions used by Popham and Narayan (1995) required 
angular velocity to be equal to the stellar one f2* at the inner boundary at 
stellar radius r*, and to the keplerian angular velocity at the outer boundary, 
specified at Rout = lOOi?*. The boundary conditions for the vertical heat flux 
had been equivalent to the Q- energy losses from (124). The radial heat flux 



from the star > was specified on its surface, what is equivalent to the 
situation described by (21C). 

Calculations have been performed for cataclismic binaries where a disk ac- 
cretion goes onto a white dwarf star. Other similar calculations (Popham et al., 
1996; Collins et al., 1998) considered white dwarfs or pre-main-sequence stars. 
Boundary layers in the accretion disks around a neutron star had not been con- 
sidered in explicit calculations, and had been treated either analytically, or by 
MAE method (Regev, 1983). Calculations of Popham and Narayan (1995) re- 
vealed the existence of a dynamical boundary layer, where there is a rapid drop 
of the angular velocity and of a velocity of the radial motion, inside the thermal 
boundary layer. There is a strong radial acceleration of matter when the flow 
approaches the dynamical boundary layer, and the density is rapidly increasing 
inside this layer. The heat flux from the boundary layer is propagating radially 
into the accretion disk and is radiated vertically on the scale of the order of the 
disk thickness. 

Time-dependent calculations for accretion disks with boundary layers around 
pre-main-sequence stars had been performed in 1-D approximation by Godon 
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(1996a, b), using nonflux thermal boundary condition dT/dr = at r = 



what is corresponding to = in (210). The results of 1-D calculations are 



qualitatively in accordance with the stationary model of Popham and Narayan 
(1995). The thermal boundary condition used above does not take into account 
a possible important heat flux from or inside the star. It seems more consistent 
to use the thermal boundary condition in the following simplified form 

where C, is the corresponding heat capacity of a star, and <&* is a heat flux on 



the star boundary from (207), which may be positive or negative, depending 
on the temperature T*. In the case of a neutron star or a white dwarf the 
heat capacity C* may be estimated from the existing models of their structure, 
and the heat capacity of the isothermal star could be considered as a rough 
approximation. 

Time-dependent 2-D model had been calculated by Kley (1991). These 
calculations show oscillations in time of luminosity which amplitude depends 
on the choice of input parameters. These oscillations could be connected with 
instabilities, but numerical effects are also cannot be excluded. 

3.1.6 Models of accretion disk with a boundary layer constructed 
by MAE method 

Matched asymptotic expansion (MAE) method was used for description of 
boundary layers in hydrodynamics (Prandtl, 1905), and was first applied for 
accretion disks by Regev (1983). The structure of the accretion disk and bound- 
ary layer with a polytropic equation of state is solved analytically by MAE 
method (Bisnovatyi-Kogan, 1994), and is described in Section 3.1. Account of 
thermal processes with the equation of state P(p, T) makes the problem more 
complicated, and analytic solution was not obtained. The equilibrium equa- 
tions describing the outer solution for the disk may be taken from the standard 
model without advection with = valid until the inner boundary of the 
outer solution at r = r*. Contrary to the standard model, where the integration 
constant jm was taken as jin = £,^k* tI with ^ 1, in the outer MAE solu- 
tion for the disk the constant ^ is slightly less than unity, and should be found 
from matching of the inner and outer solutions, like in the polytropic case (see 
Section 3.1.2). 



The thermal equations in the disk may be taken in the form (|220| ), where 



advective term may be neglected, and ( 207 ). The term with $ in this equation 



should take into account the heat flux coming from the boundary layer into 
the accretion disk, which is given by the boundary value $(r*) — ^u- The 
constant should be found from the matching to the inner solution. The 



second boundary condition for solving two differential equations (220), ( pO 



determines approach of the solution at large r to the standard one with a very 
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small $(r). It is possible to write the thermal balance equation in the disk in a 
more simple way. It is evident from the physical model, and also is confirmed 
by calculations of Popham and Narayan (1995), that the heat flux from the 
boundary layer coming into the disk is radiated from the disk on the scale 
of the order of the disk thickness h^, near the boundary layer. Instead of solving 
the differential equation, we may model the radiation flux from the boundary 
layer by an additional heating term, which is distributed over the distance /i» 
from the inner boundary by the following profiling 

^ ~ ^ exp ( ~ ^ (222) 
dr \ J 

Integration of ( |22^ ) over the radius between and oo gives the additional 
energy emitted by the disk due to the heat flux from the boundary layer equal 
to and this flux is emitted in the region 6r ^ h^, attaching the boundary 
layer. In this description there is one additional matching condition defining 
in comparison with the polytropic case. 
To obtain the inner MAE solution we should solve the equilibrium equations 



(171) and (172) together with the equations of the heat production and heat 
transfer. The constant ^6 = 1 here, like in the polytropic case from the matching 
condition for uj = 1 a,t x — oo, equal to the inner keplerian value of the outer 
solution. Matching the densities of the inner and outer solutions determines 
the inntegration constant ^ of the outer solution. Application of MAE to the 
solution of the boundary layer with thermal processes meets with a definite 
problem. Using ( |20g| ) in ( ^07|) , we write the equation of a heat transfer in the 
boundary layer in a form 

$ = -Anrh- — = — [(1 - LUAf - (1 - io)^]. (223) 

Snp dr Zr^, 

It follows from ( |223| ) that the heat flux in the boundary layer consists of the 
constant value $0 in the first term to the left, and variable part ^ (1 — w)^. Exis- 
tence of the constant flux leads to non-physical behaviour of the temperature at 
large x in MAE method, leading to unrestricted fall of the temperature. There- 
fore, in the papers of Regev (1983), Regev and Bertout (1995) the solution was 
obtained at $0 = 0, what determines that all heat produced in the boundary 
layer heats the star and does not go into the accretion disk, what is opposite to 
the condition used by Popham and Narayan (1995). As was mentioned above, 
the maximum of temperature may correspond to the intermediate level in the 
boundary layer, and heat flux from it may heat simultaneously the star and 
the boundary layer. In the direct calculations it corresponds to changing of 
boundary condition, but was not realized. In MAE method solution cannot be 
obtained directly at non-zero $0. Additional modification of MAE method are 
needed for its account. 
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Note, that large deviations of the ^ value from unity, ^=0.61, obtained by 
Regev (1983) for lj* — 0.3, could indicate the step into the region, where the 
asymptotic expansion begins to be violated. When the stellar angular velocity 
approaches a critical value, the boundary layer becomes thick (see (191) and 
the MAE method fails. 



3.2 Accretion Disks Around Rapidly Rotating Stars 

Investigations of low-mass-X-ray binaries (LMXB) containing a neutron star 
have lead to the conclusion that the accreting object can rotate rapidly with a 
surface angular velocity close to Keplerian value. The question of accretion onto 
a rapidly rotating star also arises in cataclysmic variables, in which the primary 
is a white dwarf. When the stellar angular velocity reaches its critical value 
the maximum of the angular velocity and the boundary layer disappear, and a 
regime of accretion is changing. Disk accretion into a rapidly rotating star, with 
the Keplerian angular velocity at the equator, had been investigated by Popham 
and Narayan (1991), Paczynski (1991), Colpi et al. (1991). Self-consistennt 
regime was found by Bisnovatyi-Kogan (1993a), in which star remains to rotate 
rapidly during the accretion process. 



3.2.1 Self-consistent model 

When a star reaches the critical rotation rate for which the the surface angular 
velocity becomes Keplerian, the flux of angular momentum into the star leads to 
differential rotation inside the star. It is clear that disk properties far away from 
the star do not change when the star reaches the critical rotation rate, so that 
the disk remains Keplerian. The properties of the intermediate zone between 
the star and the disk are determined by the stellar rotational velocity law, but 
we do not need to know exactly the properties of this transition zone in order to 
build a self-consistent model. We only have to know that during accretion onto 
a rapidly rotating star, the transition layer is smooth, all values are monotonous 
and there is no enhanced heat production in this layer, in contrast to the case 
of boundary layers between the accretion disk and a slowly rotating star. 

The evolution of the star under mass accretion may be characterized by the 
variations of two values: the mass M and the total angular momentum J and by 
one function, j{r), determining the angular momentum distribution within the 
star; j{r) is determined by the viscosity law inside the star. For a given mass 
flux M, determined by conditions far away from the star, the star gradually 
increase its mass and angular momentum. When the stellar rotation rate is 
much smaller than the Keplerian limit the speciflc angular momentum of matter 
accreting onto the star is approximately equal to MvkcRsc, (vk = {GM/r^y^'^ 
is the Keplerian velocity, Rsc- the stellar equatorial radius, index "e" referes 
to values at the stellar equator). When the stellar rotational velocity at the 
surface is smaller than the Keplerian value, the angular velocity of the disk has 
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a maximum near Rsc] at this point, the viscous flux is zero and the momentum 
flux is determined only by the convective term Mjc, j — v^pr. The maximum of 
the rotational velocity disappears, when the star rotation is close to the critical 
velocity and the (negative) viscous momentum flux becomes important. 

The demand of self-consistency during accretion onto a rapidly rotating star 
may be formulated as the condition that the star absorb the accreted mat- 
ter with a specific angular momentum such that the star remains in 
the state of critical rotation. Depending on the viscosity law, a critically 
rotating star may have different distributions j{r) and in general, the structure 
of the star gradually changes with time. If viscosity inside the star is high, 
the star rotates almost rigidly; this case has been considered in the above men- 
tioned papers. If the star-disc system can be described by a polytropic equation 
of state, there is a self-similar solution for such a system which does not depend 
on mass. 

3.2.2 Self-similar solution for a polytropic star-disc system in the 
case of a rigidly rotating star 

The equations describing the structure of a rotating star with a polytropic equa- 
tion of state (|4^) the equilibrium equations in spherical coordinates are (Chan- 
drasekhar 1989) 

-|^ + ?-"Ml-M^) = (224) 
p or or 

i^ + ^ + f7V^ = (225) 
p Ofi op. 

and the equation for the gravitational potentional $ is 

/ ^^ _^2)a$\ ^^^^ (226) 

r^ dr \ dr J dp \ dp J 

Here p = cos 9 and H. is the angular velocity of the stellar matter. In general 
0, is not constant throughout the star, and for a polytropic star in equilibrium 
depends only on the cylindrical radius rs'mO. Excluding $ from Eqs.(2)-(4), we 
get 

1 d fr^dP\ 1 d fl-p^dP 



r-^ dr \ p dr J r'^ dp \ p dp 

ATrGp--^{x^n^)^ (227) 
X ax 



where x = rsinO = r^ 1 — p"^, = il,{x). Introduce the nondimcnsional poly- 
tropic variables 

e-r/n, e" = p/pc, Lo = n/n,, (228) 
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where pc is the central density of the star and r* and f2* are given by 



r* 



{n + l)Kp, 



, = v^27rG/9c (229) 



In the nondimensional variables (6)-(7) Eq.(5) for Q= const, can be written in 
the form 

The solution of Eq.(8) is uniquely determined by the boundary conditions 

e = i, ^ = at e = o, 

BP) 

Vl-M^ =0 at /i = 0, M = ±1, (231) 

and by the value of lu: one also requires that the solution must satisfy the original 
equilibrium Eqs.(2),(3). The latter condition results from the fact that Eqs.(2) 
and (3) were differentiated in order to obtain (8), so that a solutions of (8) need 
not satisfy (2), (3). It is in fact necessary to find a gravitational potential $ 
from (4) which, when written in nondimensional form </> = ^/AirGpcT^ satisfies 
the scaling conditions of the self-similar solution. For compressible matter with 
n > a solution exists only if w < uJc] for uj = ujc the centrifugal force on 
the equatorial boundary of the star exactly balances gravity (see the proof in 
Bisnovatyi-Kogan & Blinnikov 1981). Numerical calculations have been done, 
which gave the structure of rotating polytropic stars for different u up to the 
critical value uj^ (James 1964; Blinnikov 1975; Ipser & Monagan 1981; see also 
Tassoul 1978). 

For a given n the structure of a critically rotating star does not depend on 
its mass. By definition, the mass M and the total angular momentum J are 
given by 



/I /•4out(A') 
^dnj^ e"(^,M)€2^^, (232) 



/I /•Jout(M) 
{l-IJ?)diij^ e"(e,/x)^^de (233) 

Introducing nondimensional values of the mass M.n and of the total angular 
momentum J7„, defined as 

Mn = ^J ^dfij^ e"(^, ^^)ed^, (234) 
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= 2 \jy - I m)^''^ (235) 

M and J may be written as 

M = Anp.rlMn, (236) 



J = iirp.rl ^27rGpcJn ■ (237) 

The average specific angular momentum of tfie star js and the derivative along 
the critical states ja = {dJ/dM)c are written, using (7), (14), as 

5- 2n 



j, = J/M = rl^2nGp^Jn/Mn ia = ^ is- (238) 

O 71/ 

Using (14) one may write pc and r*, which appear in (7), as a finction of M 



AttG \~ f M \~ ( A-kG \~ f M 

n = 



{n + l)KJ \'i-KMnJ ' \{n+l)KJ \'i-KMr. 

(239) 

The specific angular momentum of matter je at the stellar equator may be 
written as 

= = rl^2TTGp^iluj, = r^/r. = €out(7r/2). (240) 
From a comparison of (16) and (18), it follows that the ratio 

ia ^ 5 - 2n js ^ 5 - 2n J„ ^^^^^ 
je 3 - n jc 3 - n Mn^lu) 

does not depend on the stellar mass M, and is a function of n and to only. 
The nondimensional values of the angular velocity Wc, the equatorial radius 
^out(^ = 7r/2) = ^o, the mass M.n, the momentum of inertia around the rota- 
tional axis In = Jn/i^c, the average specific angular momentum of the star , 
the nondimensional derivative Ca 

c; _ Or? 

Cs=jJ{rlV^^c)=Jn/Mn, Ca = ^ Cs, (242) 

3 — n 

and the angular momentum of matter on the stellar equator 



Ce = je/{rl \/27rGpc) = '^c (243) 

for several polytropic indices n in the state of critical rotation are given in ta- 
ble 1, which is based on calculations of James(1964), Blinnikov(1975) and Ipscr 
and Managan(1981). The ratio of the equatorial and polar radii £,e/£,p and the 
nondimensional parameters for nonrotating polytropes (outer radius ^on) mass 
M-on and momentum of inertia around the symmetry Chandrasekhar 
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Table 1: Parameters of non-rotating and critically rotating polytropic stars 



n 










5o 


«p 










0.5 


2.7528 


3.7871 




0.389 
(0.367) 




2.51 










0.6 








0.362 
(0.354) 




2.29 










0.808 








0.326 


4.7652 


1.917 


5.0248 


24.43 


1.585 


0.3304 


1.0 


IT 


7r 


8.104 


0.289 


4.8265 


1.792 


4.289 


18.73 


1.262 


0.2805 


1.5 


3.65375 


2.71406 


7.413 


0.209 


5.36 


1.626 


3.2138 


12.18 


0.792 


0.176 


2.0 


4.35287 


2.41105 


7.074 


0.147 


6.307 


1.555 


2.6518 


9.62 


0.533 


0.0894 


2.5 


5.35528 


2.18720 


7.013 


0.09965 


7.7623 


1.522 


2.30563 


8.49 


0.367 


0.0 


3.0 


6.89685 


2.01824 


7.234 


0.0639 


10.123 


1.535 


2.0743 


8.125 


0.250 


— oo 



1957; James 1964; Blinnikov 1975) are also given. 

In nondimensional variables it is obvious that the derivative of the angular 

momentum of the star along the sequence of critically rotating states is always 
smaller than the equatorial value Qc] Cc is equal to the Keplerian value of a 
critically rotating non-spherical star. The viscous flux of specific angular mo- 
mentum in the accretion disk around a critically rotating star jv is therefore 
negative; near the surface of the star is equal to 

jv = ia - jo (244) 

as one requires that the total angular momentum flux into the star is equal to 
Mja for self-consistency. In the self-consistent model of accretion disk around 
the rapidly rotating polytropic star the value of the constant ^ in the formulae 
of the Section 1.4 should be taken equal to from the Table |^. 

3.2.3 General (nonpolytropic) case of disc accretion into a rapidly 
rotating star 

In realistic cases, thermal processes are very important in the disc, polytropic 
approximation is not good, and equations of the thermal balance from Section 
1.5 should be taken into account. Transition to the rapidly rotating star is 
accompanied by a substantial change of the integration constant in (|7l|), (re- 
frefl.ll), what leads to impotrant observational consequences. 
The total luminosity of the accretion disc is equal to 

r°° S CM 
L = / Ai^Frdr ^ - P)M • (245) 

When /3 = 1 the star accrets matter with Keplerian angular momentum, and 
only half of the gravitational energy of the accreted matter is radiated from the 
disc, according to the virial theorem. The fate of the remaining half is different. 
In the case of a black hole, it is absorbed, increasing the mass and angular 
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momentum of the black hole without any significant radiation; however matter 
falling from the last stable orbit to the horizon still emits radiation (Bisnovatyi- 
Kogan & Ruzmaikin 1974). 

In the case of accretion onto a (slowly rotating) neutron star there are several 
possibilities. If the radius of the neutron star Vg is smaller than 3rg, then 
Ti = 3rg; the disk luminosity is the same as in the case of a black hole, and 
the remaining gravitational energy, including the part gained during free fall 
onto the neutron star surface, is emitted close to the neutron star surface. If 
Tg > 3rg, then n = ''s, the disc luminosity is GMM /2rs and an equal amount 
of energy is emitted in the boundary layer near the stellar surface, where the 
accreting matter converts its kinetic energy into heat. 

The situation changes gradually while the central object absorbs angular 
momentum. The black hole moves from Schwarzschild to extreme Kerr solu- 
tion; the radius of the last stable orbit decreases from 3rg to rg, and the radius 
of the horizon changes from rg to rg/2. The efficiency of the disc energy emis- 
sion increases from 5.7% to 42% (Novikov & Thorne 1973). The changes go 
in an opposite direction when a neutron star accelerates its rotation. Let us 
consider only the case Vg > 3rg. The luminosity of the disc decreases because 
the stellar radius increases. The fraction of energy emitted by the neutron star 
also decreases. It changes from GMM /2rs in the case of a nonrotating star to 
the difference between rotational energy of matter in the disc and at the stellar 
surface. 

A rapid change in the efficiency of energy release and disc luminosity happens 
when the star rotation reaches its maximum value. The energy emitted near the 
stellar surface tends towards zero, but the luminosity of the disc suffers a drastic 
change. The star does not absorb any longer all of the angular momentum 
coming from the disc, but absorbs only the fraction required to maintain the 
star in a critically rotating state. The specific angular momentum of matter 
accreted by the star is equal to {dJ/dM)c and the remaining part is carried 
away in the disc by viscous stresses (see sec. 2). Viscosity carries not only angular 
momentum, but also energy, so that the energy production and luminosity of 
the disk rapidly increase, from the valu e cor responding to /3 = 1 to the value 
corresponding to f3 — ja/jc according to (245); for polytropic stars the values of 
f3 are given in the last column of table ^ . This implies a rapid increase of the 
total luminosity by a factor 2-;- 3, depending on the stellar structure. This rapid 
increase is easy to understand if one remembers that when a star accelerates 
its rotation, part of the gravitational energy is converted into rotational energy 
without heat production. When a star reaches the limiting rotation, the growth 
rate of its rotational energy strongly diminishes and a correspondingly larger 
fraction of gravitational energy is transformed into heat. The minimum of 
luminosity of an accreting neutron star is reached when the stellar angular 
velocity is slightly below the critical value. A rapid increase in luminosity must 
be accompanied by a corresponding hardening of the emitted spectrum because 
the energy release is increased, raising the effective temperature in the inner part 
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of the disc. Such events may be expected in objects hke LMXB or cataclysmic 
variables. The young T Tauri stars are born rapidly rotating, so their value 
of /3 is small from the beginning and close to that of a polytrope n = 1.5 
with P = 0.176 given by Table |l| Note, that whereas the standart solution 
of accretion disk (Shakura, 1972; Shakura and Sunyaev, 1973; Novikov and 
Thorne, 1973) is not valid in regions close to the inner edge in the case of a 
black hole (see Paczynski & Bisnovatyi-Kogan 1981), the self-consistent solution 
for accretion onto a rapidly rotating star, may be used down to the very stellar 
radius, which is the inner edge of the disc. 



4 Accretion flows in presence of magnetic fields 

4.1 Magnetic Field Amplification in Quasi-Spherical Ac- 
cretion 

Matter flowing into a black hole from a companion star or from the interstellar 
medium is very likely to be magnetized. Due to the more rapid increase of a 
magnetic energy in comparison with any other kind of energy, the dynamical 
and any other role of the magnetic field is becoming more and more important 
when matter flows inside. Schwartzman (1971) was the first to emphasize the 
importance of the magnetic field for an energy release during accretion into 
a black hole. He showed that due to more rapid growth of magnetic energy 
its density Em approaches a density of a kinetic energy Ek, and proposed a 
hypothesis of equipartition Em « Ek, supported by continious annihilation of 
the magnetic field, for the smaller radiuses of the flow where the main enegry 
production happens. This hypothesis is usually accepted in the modern picture 
of accretion (see e.g Narayan and Yu, 1994). Schwartzman (1971) considered 
an averaged quasistationary picture with local equipartition, another variant 
of magnetic accretion, where equipartition and magnetic field annihilation is 
accompanied by formation of shock waves was considered by Bisnovatyi-Kogan 
and Sunyaev (1972), see also Chang and Ostrikcr (1985). More accurate account 
of the heating of matter by magnetic field annihilation was done by Bisnovatyi- 
Kogan and Ruzmaikin (1974), where exact nonstationary solution for a field 
apmlification in the radial accretion flow was also obtained. 

Consider the general problem of the field amplification in the accretion fiow 
where gas has an angular momentum and field annihilation is approximated by 
a finite "turbulent" conductivity. This approach should reproduce a smooth 
averaged flow, considered by Schwartsman. In the known stationary velocity 
field v(r, 0) the magnetic field B is described by Maxwell equations in the MHD 
approximation (without the displacement current), with a finite conductivity 

47r 1 dB 
VxB = — j, VxE = — , V-B = 0, 
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j=o- E+-vxB 



(246) 



where a is the plasma conductivity, j is the electrical current density, and E 
is the electrical field strength. Consider an axisymmetrical pictire in a spher- 
ical coordinate system r, 0, where all derivative are zero, but all vectors 
V, B, E, j may have nonzero all t hree components. The equation containing 
only magnetic field, following from ( |246| ) 



f = Vx|vxB| + Av'B. 



(247) 



These equations are completed by the condition of zero divergence of the mag- 
netic field, which for the choosen case has a view 



1 d 



1 



d 



• sin d do 



{sine Bg) ^ 0. 



(248) 



4.1.1 Time-dependent amplification of a Magnetic field in Spherical 
Accretion 



For a spherical accretion with v = {vr, 0, 0) the equation (247) reduces to 



dt 
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— r-^T^(sin6'urBe) + 
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J.2 Qj, 



,dB^ 
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sin^ I 



For initial nonstationary stage of the field amplification, when field dissipation is 
negligible, an alyti cal solution may b e obt ained (Bisno vaty i- Kog an & Ruzmaikin, 
1974). From ( ^49|) with account of (HI), and from (|25C|) , (p5ll) we get 



d{r^Br) 



0, 



d{rVrBe) 



0, feM^o, 



dt dt dt 

where the full Lagrangian derivative is determined as 



(252) 
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A solution of ( ^52| )- ( |253| ) is determined by 4 first integrals of the characteristic 
equations 

r dj. 

Ci=t- 0-2 = r^Br, C3 = rVrBe, C4 = rVrB^. (254) 

J Vr 



For a free fall case with Vr = —\ 2^iM -^ye get 



2 r3/2 

Ci=t + -^^. (255) 

3 V2GM 

The initial condition problem at t is solved separately for poloidal and 
toroidal fields. For initially uniform field i?rO = BqCOsO, Bgo — — i3osin6' we 
obtain the integrals 

2 



Ci = - C2 = r^Bo cos 61, C3 = \/2GMrBo sin 61. (256) 

3 v2GM 



The relation between the integrals found at initial moment after excluding r{C 



is valid at any time, so finally we get from (254) the solution (Bisnovatyi-Kogan, 
Ruzmaikin, 1974) 

Br = ^^{ry^ + ltV^y\ (257) 

Sosin6' / 3 , 

Bg = f r'^/^ + -tV2GM\ . 

The radial component of the field is groving most rapidly. It is ^ for 
sufficiently large times, '-^ t*/^ at given sufficiently small radius, and is groving 
with time everywhere. 

For initially dipole magnetic field Bro = ^^^^1^, Bgo = - %lf^ we obtain 
the following solution 



Br = — ( r^/' + -W2GM ] , (258) 

-5/3 



Bg = -f^^iry' + ltV2GM] 



Here the magnetic field is decreasing everywhere with time, tending to zero. 
That describes a pressing of a dipole magnetic field to a stellar surface. The 
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azymuthal stellar magnetic field if confined inside the star. When outer layers of 
the star are compressing with a free-fall speed, then for initial field distribution 
B^o = BQr"f{6) the change of magnetic field with time is described by a relation 

B^ = - (^r3/2 + ^-tV2GM^ . (259) 

4.1.2 Amplification of magnetic Field in Rotating Accretion Flows 

In a quasi-spherical accretion flow of a perfectly conducting plasma, rf(B • 
dS)/dt = 0, where dS is a surface area element and d/dt is the derivative 
following the flow (v = VrV + v^cf) , in spherical coordinates). With dS = r^dflr, 
Vr ~ ?;r(|r|), and dO, a flxed solid angle increment, one finds that Br oc 
independent of the toroidal velocity v^. Writing Bj. = Bo (''out /^)^, the radius 
at which = B^/{STr) is 



' eqm 



\VgmmJ 



2/3 



(260) 



where Vk = {GM/ry^^ and the accretion speed is assumed Vr = —gvK, with 
g = const < 1. For example, for a black hole mass M = 10^ Mq, M = O.lM0/yr, 
Bq = 10~'^G, To ~ Ipc, and g ~ 0.1, we find requi ~ 6 x IQ-"^^ cm, which is much 
larger than the Schwarzschild radius, rs « 3 x lO^'^ cm. Thus the magnetic field 
resulting from fiux-freezing is dynamically important in quasi-spherical accretion 
fiows (Shwartsman 1971), and this is independent of dynamo processes or MHD 
instabilities. Further accretion for r < r^qui is possible only if magnetic flux 
is destroyed by reconnection and the magnetic energy is dissipated. Figure 7 
shows a sketch of the instantaneous poloidal magnetic field configuration. 

For example, for a flow with time-averaged velocity v = Vri + v^(f) , an ap- 
proximate solution for the time-averaged magnetic field for infinite conductivity 
can readily be obtained in the vicinity of the equatorial plane. In this region it 
is reasonable to assume B = B^r + B^<p ; that is, Bg is neglected. Thus 

V-B = 0= L^(r'Br) + ^^{B,) . 
Near the equatorial plane, sin^ w 1. This equation is then satisfied by taking 

Br = exp (ij dr kr + im0 J , (261) 



= f{r) exp ^i J dr kr + imcf) ) , (262) 
with the actual fields given by the real parts and with 
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k • B = krBr + k4,B4, = , 



(263) 



where k = rkr+<p k^,k^ = m/r is the azimuthal wavenumber, and m = ±1, ±2, 
etc. The stationary magnetic field satisfies Vx(v x B) = 0, and this requires 

V X B = 6{v^Br - VrB^) = , or B^ = {v^/Vr)Br . 
That is, V cx B and k • v = 0. Owing to equation (3), 

K = -{v0/vr)k^ . (264) 

In ADAF accretion models, w^/wr = const, so that kr oc 1/r. Thus, \Br\ = 
\vr/v^\\B^\ CX l/r^, and /(r) cx l/r^. 

The nature of the field is most easily obtained from the flux function rAe 
(with A the vector potential), where 



Br 



1 dirAg) 



1 djrAe) 
r dr 



for sin^ « 1. Note that (B • W){rAe) = 0. Thus the frozen-in field is described 
by 

rAg = const exp \ im — In r + ) , (265) 

where it is assumed that > 0, t^r < 0, and v^/vr =const. The fact that 
IrAffl =const corresponds to the conservation of the ingoing or outgoing flux in 
a given tube. A given fleld line satisfles 



r = const exp 



(266) 



Figure 2 shows the nature of the equatorial field for m = 2. 

Consider now the problem of magnetic field amplification in an accretion 
fiow which is not infinitely conducting. The flow is described by the magneto- 
hydrodynamics (MHD) equations. 



^-7- = -Vp + pg + -J X B 
at c 



An 

V X B = — J 

c 



V X E = - 

J = cr(E + V X B/c) 



IdB 

"c'dt ' 



(267) 



where v is the flow velocity, B the magnetic field, p the plasma pressure, a 
the plasma conductivity, r] the dynamic viscosity (with y = r]/ p the kinematic 
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viscosity), J the current density, and E the electric field. These equations can 
be combined to give the induction equation, 

OB 

— =Vx (vxB)-Vx(77™VxB) , 

V X (v X B) + r],n\7^B , (268) 

where rjm = /{Aira) is the magnetic diffusivity, and the approximation involves 
neglecting V?7„ 



In equations (267) and ( |268| ) both the viscosity v and the magnetic diffusivity 
rjm have the same units and both are assumed to be due to turbulence in the 
accretion flow. Thus, it is reasonable to express both transport coefficients using 
the "alpha" prescription of Shakura (1972) and Shakura and Sunyaev (1973), 



(269) 



where it is the outer scale of the turbulence, and Cs = \/p/p is the isothermal 
sound speed. Bisnovatyi-Kogan and Ruzmaikin (1976) introduced am and pro- 
posed that a„i ^ a. Note that a and am characterize a turbulent MHD flow in 
which there is a Kolmogorov cascade of energy from large scales {£t) to much 
smaller scales where the actual (microscopic) dissipation of energy occurs. 

We obtain stationary B field solutions for accretion flows with time-averaged 
flow V = rvr + 4> and with finite diffusivity r]m in the form of equations ( ^61[ ) 
and (^62[ ) , but with complex, = + ikj.. Equation (262) still applies, 
but (264) does not. In its place we have from equation (^68[), 



« ik X (v X B) - ?7„ik^B 



(270) 



where the approximation involves neglect of terms of order l/(|k|r) ^ 1 or 



smaller. Solution of (270) gives 



k^Vr + k^v^ , and k^. 



(271) 



where in this equation k-^ = (^i?)^ + ^0- Consistent with the first part of (12), 
we take k^ = rir/r, where rir — const. For the AD AF fl ows, we write the second 
part of (p69h as rjm = amCsT. The second part of (271) then gives 



(272) 



where the quantity amCs/\vr\ — const. 

The imaginary part of kr gives an additional radial dependence of Br and 
(and rAg), multiplying the earlier expressions (261),( p62|) , and ( ^65|) by a factor 
r^, where S = (amCs/|wr|)(?ir + ™^)- 



For radial distances r < r^qui, we have 
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equipartition with magnetic energy density Em oc r~^/^ so that |Bp| oc r~^/^. 
This requires 5 = 3/4 so that {rAel oc r^/^, which corresponds to the flux within 
a given tube decresing as r decreases. The accretion speed is 

l^rl = ^{nl + m'^)amCs ■ (273) 

Note that for the present solution the accretion depends on the magnetic dif- 
fusivity r]m but not the viscosity v; more general solutions have Vr dependent 
on both viscosity and diffusivity. Validity of this solution requires |A;^|/|k| = 
(3/4)/-\/n^ + <C 1. The field line pattern is similar to that in Figure 2, but 
the number of field lines in a tube a r^/^. 

The magnitude of the magnetic field follows from the requirement that 
the inflow of angular momentum carried by the matter, L^ech 
sin^ > pVrV^ « AiTr^pVrV^ < 0, equal the outflow carried by the fleld, Lmag = 
—r^BrB^ > 0, which gives pw^ w B'^/Att. More generally, turbulent viscosity 
carries away part of the angular momentum in which case the magnetic field 
strength is reduced. For r < Vequi, we have Bj. = BQ{requi/f)^^^-, and zero 
angular momentum flux gives 

1 /2 

This gives So w 55 G for r^qui = 10^^ cm, M = IQ^Mq, M = O.lMe/yr, 
am = 0.1, Cg/v^ = 0.1, and + = 10^. 

4.2 Stationary picture of the magnetic field distribution 
in quasi-spherical accretion: turbulence, equipartition 
and field annihilation; two-temperature advective 
discs 

In the optically thin accretion discs at low mass fluxes the density of the matter 
is low and energy exchange between electrons and ions due to binary collisions 
is slow. In this situation, due to different mechanisms of heating and cooling for 
electrons and ions, they may have different temperatures. First it was realized 
by Shapiro, Lightman, Eardley (1976) where advection was not included. It was 
noticed by Narayan and Yu (1995) (see also Ichimaru, 1977), that advection in 
this case is becoming extremely important. It may carry the main energy flux 
into a black hole, leaving rather low efficiency of the accretion up to 10^^ — 10~'* 
(advective dominated accretion fiows - ADAF). This conclusion is valid only 
when the effects, connected with magnetic field annihilation and heating of 
matter due to it are neglected. 

In the ADAF solution the ion temperature is about a virial one fcTj ~ 
GMnii/r, what means that even at high initial angular momentum the disk 
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becomes very thick, forming practically a quasi-spherical accretion flow. It is 
connected also with an "alpha" prescription of viscosity. At high ion temper- 
atures, connected with a strong viscous heating, the ionic pressure becomes 
high, making the viscosity very effective. So, due to suggestion of "alpha" vis- 
cosity in the situation, when energy losses by ions are very low, a kind of a 
"thermo-viscous" instability is developed, because heating increases a viscos- 
ity, and viscosity increases a heating. Developement of this instability leads to 
formation of ADAF. 

A full account of the processes, connected with a presence of magnetic field 
in the flow, is changing considerably the picture of ADAF. It was shown by 
Schwarzman (1971), that radial component of the magnetic field increses so 
rapidly in the spherical flow, that equipartition between magnetic and kinetic 
energy is reached in the flow far from the black hole horizon. In the region 
where the main enegry prodiction takes place, the condition of equipartition 
takes place. In presence of a high magnetic field the cfficiecy of a radiation 
during accretion of an interstellar matter into a black hole increase enormously 
from ~ 10~^ up to ^ 0.1 (Schwarzman, 1971), due to efficiency of a magneto- 
bremstrahlung radiation. So possibility of ADAF regime for a spherical accre- 
tion was noticed long time ago. To support the condition of equipartition a 
continuous magnetic field reconnection is necessary, leading to annihilation of 
the magnetic flux and heating of matter due to Ohmic heating. It was obtained 
by Bisnovatyi-Kogan and Ruzmaikin (1974), that due to Ohmic heating the 
efficiency of a radial accretion into a black hole may become as high as ~ 30%. 
The rate of the Ohmic heating in the condition of equipartition was obtained 
in the form 

rr^dS 3 

^^=-7T7^ • (275) 

dr 2 87rpr ^ ^ 

In the supersonic flow of the radial acrretion equipartition between magnetic 
and kinetic energy was suggested by Schwarzman (1971): 

|! « ^ = (276) 

For the disk accretion, where there is more time for a field dissipation, almost 
equipartition was suggested (Shakura, 1972) between magnetic and turbulent 
energy, what reduces with account of "alpha" prescription of viscosity to a 
relation 

(277) 

where am characterises a magnetic viscosity in a way similar to the turbulent a 
viscosity. It was suggested by Bisnovatyi-Kogan and Ruzmaikin (1976) the sim- 
ilarity between viscous and magnetic Reynolds numbers, or between turbulent 
and magnetic viscosity coefficients 
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i?e„ = ^, (278) 

V Vm 

where the turbulent magnetic viscosity r/m is connected with a turbulent con- 
ductivity a 

a = (279) 
Taking T]m = ^'^i S^t a turbulent conductivity 

a = t'? = ^ (280) 

in the optically thin discs. For the radial accretion the turbulent conductivity 
may contain mean free path of a turbulent element It , and turbulent viscocity Vt 
in ( |28C| ) instead of h and Vg- In ADAF solutions, where the ion temperature is 
of the order of the virial one the two suggestions ( |276| ) and (277) almost coinside 
at am ~ 1- 

The heating of the matter due to an Ohmic dissipation may be obtained 
from the Ohm's law for a radial accretion in the form 

T— = = — (281) 

what coinsides with ( ^75[ ) when am = or It — ^^^J^ ■ Here a local 

electrical field strength in a highly conducting plasma is of the order of£ ^ , 
Ve ~ Vt ~ aVs for the radial accretion. 

Equations for a radial temperature dependence in the accretion disk, sepa- 
rate for the ions and electrons are written as 



dEi Pi dp 
~dr ~ 'f^'dt 



— Ti-jji + Ti-Bi — Qie , (282) 



,, 9 7, — '^Tje ^ "^Be ^ Qie ^brem ^cyc i (283) 

dt p'^ dt 

Here ^ = ^ + ''^r ^ • A rate of a viscous heating of ions Ti^i is obtained from 



as 

M 2 r V 

The equality in the last equation is related to binary collisions. Combining ([70|), 
(|2l),(||),(^, we get I 

vkJ vk AnaV^ r^vi 
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where vk = rflK, J = 1 — The rate of the energy exchange between 
ions and electrons due to the binary colhsions was obtained by Landau (1937), 
Spitzer (1940) as 

milTLe \me nii J 

with £nA = 0(20) the Coulomb logarithm. Thermodynamic functions must 
take into account a possible relativistic effects for electrons, which energy may 
exceed rrieC^. Neglecting pair formation for an low density accretion disk, we 
may write an exact expression for a pressure 

Pg=P,+P,= n,kT, + n.kTp = n,k{Te + Ti), (287) 

and an approximate expression for an energy, containing a smooth interpolation 
between nonrelativistic and relativistic electrons 

E.,E.^E..t^4±^!^^l'Jl. (288) 

The electron bremstrahlung Cbrem and magneto-bremstrahlung Ccyc cooling 
of maxwellian semi-relativistic electrons, with account of free-bound radiation 
in nonrelativistic limit, may be written as by inrerpolation of limiting cases 
(Bisnovatyi-Kogan and Ruzmaikin, 1976) 

777-gC ~r f^J- e 



„ 1.7xlO»„..c^ -H.7 xlO'^(tr.)^ ^^^ 

TOgC + kle 



Effects of the cyclotron self-absorption may be important for nonrelativistic 
electrons, where radiation itself is low. With increasing of the electron temper- 
ature self-absorption effects decrease, because the black-body emission increase 
with a temperature (~ T^) more rapidly then magneto-bremstrahlung volume 
losses. Account of self-absorption was investigated by Trubnikov (1973). 

In the case of a disk accretion there are several characteristic velocities, vk, 
Vs, and Vt = avs, all of which may be used for determining " equipartition" 
magnetic energy, and one characteristic length h. Consider three possible choices 
with = vj(, v^, and for scaling — Anpv^, note that choice of cross 
products of different velocities and of Vs , is also possible. We get the following 
expressions for magnetic field strengths in different cases, with account of ( p85[ ) 

a. vb^vk, t;— = ^7^, ' 



8tt 2 ' \aV2 
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pvl „ / Ma Vs\^''^ 



c. - wt, — = — B^a„A — j=— — 



1/2 



The expression for an Ohmic heating in the turbulent accretion disk also may 
be written in different ways, using different velocities ve in the expression for 
an effective electrical field 

£ = (292) 
c 

A self-consistency of the model requires, that expressions for a magnetic heating 



of the matter TIb, obtain ed fr om the condition of stationarity of the flow (275), 



and from the Ohm's law (281), should be iden tical . That gives some restric tions 



for t he c hoice of a charac teris ti c vel ocity in (292). Comparison between ( 275) 



and (281) with account of ( |28C| ),( |285| ) shows the identity of these two expressions 
at 



a 



JOLr, 



(293) 



That implies w r ~ ~ °"l-r — < Vt So, the model is becoming self- 

vkvJ vkVJ 

consistent at the reasonable choice of the parameters. Note, that in the advective 
models J7 is substituted by another function which is not zero at the inner edge 
of the disk. The heating due to magnetic field reconnection Hb in the equations 



(284) as 



(282), (|283|), obtained from (|275|) and (291), may be written with account of 



Hb - ——^r ^ —^v^ ( 7^ ) • (294) 



So, at vb = vk the expressions for viscous and magnetic heating are almost 
identical. The distribution of the magnetic heating between electrons and ions 
has a critical influence on the model, if we neglect the influence of a plasma 
turbulence on the energy relaxation, and take into account only the energy 
exchange by binary collisions from ( ^86|) . Observations of the magnetic field 
reconnection in the solar flares show (Tsuneta, 1996), that electronic heating 
prevails. 

It follows from the physical picture of the field reconnection, that transfor- 
mation of the magnetic energy into a heat is connected with the change of the 
magnetic flux, generation of the vortex electrical held, accelerating the parti- 
cles. This vortex field has a scale of the turbulent element and suffers rapid 
and chaotic changes. The accellerating forces on electrons and protons in this 
fields are identical, but accelerations themselfs differ ~ 2000 times, so during 
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a sufRciently short time of the tubulent pulsation the electron may gain much 
larger energy, then the protons. Additional particle acceleration and heating 
happens on the shock fronts, appearing around turbulent cells, where reconnec- 
tion happens. In this process acceleration of the electrons is also more effective 
than of th e pro to ns. I n the paper of Bisnovatyi-Kogan and Lovelace (1997) the 
equations (282), (283) have been solved in the approximation of nonrelativistic 
electrons, vb=vk, what permitted to unite a viscous and magnetic heating into 
a unique formila. The combined heating of the electrons and ions were taken as 



Tie = (2- g)n,r,, He - gH^j^. (295) 

In the expression for a cyclotron emission self-absorption was taken into account 
according to Trubnikov (1973). The results of calculations for g — 0.5—1 show 
that almost all energy of the electrons is radiated, so the relative efficiency of 
the two-temperature, optically thin disk accretion cannot become lower then 
0.25. Note again that accurate account of a plasma turbulence for a thermal 
relaxation and corresponding increase of the term Qie may restore the relative 
efficiency to its unity value, corresponding to the optically thick discs. 

Studies of the solar corona have led to the general conclusion that the energy 
build up in the chaotic coronal magnetic field by slow photospheric driving is 
released by magnetic reconnection events or flares on a wide range of scales and 
energies (Benz 1997). The observations support the idea of "fast reconnection," 
not limited to a rate proportional to an inverse power of the magnetic Reynolds 
number (Parker 1979, 1990). The data clearly show the rapid acceleration of 
electrons (e.g., hard X-ray flares, Tsuneta 1996) and ions (e.g., gamma ray line 
events and energetic ions, Reames et al. 1997), but the detailed mechanisms of 
particle acceleration are not established. 

An essential aspect of the coronal heating is the continual input of energy 
to the chaotic coronal magnetic field due to motion of the photospheric plasma. 
As emphasized here, there is also a continual input of energy to the chaotic 
magnetic in a quasi-spherical accretion flow due to compression of the flow. 
The fact that the ratio of the plasma to magnetic pressures is small in the 
corona but of order unity in the accretion flow may of course affect the details 
of the plasma processes. 

In magnetic field reconnection plasma flows into the neutral layer from above 
and below with a speed of the order of the turbulent velocity vt ■ Magnetic flux 
is destroyed in this layer as required for the reasons discussed earlier. Con- 
sequently, there is an electric field in the direction = —{l/c)dAz/dt = 
0{vtBo/c), where Bq is the field well outside of the neutral layer. This is 
the same as the estimate of the electric field given by Bisnovatyi-Kogan and 
Lovelace (1997). In the vicinity of the neutral layer {\y\ ^ £t) this electric 
field is typically much larger than the Dreicer electric field for electron runaway 
Ed — 47re^ (rie/fcTg) In A (Parail & Pogutse 1965), and this leads to streaming 
motion of electrons in the z— direction. Note that it is much more difficult to 
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have ion runaway because of the large gyro radii of ions (Bisnovatyi-Kogan and 
Lovelace, 1997; Lesch, 1991). Of course, in the almost uniform magnetic field 
assumed by Quatacrt (1998) there is no nmaway of particles. 

The streaming of the electrons can give rise to a number of different plasma 
instabilities the growth of which gives an anomalous resistivity. Galeev and 
Sagdccv (1983) (see also Li et al. 1999) discuss the quasi-linear theory for such 
conditions and derive expressions for the rate of heating of electrons, 

dT, 1 f d^k k-Ue 

7k Wu , (296) 



dt nj (27r)3 ' Wk 
and ions, 

dTi 1 f d^k 

-dF = -nJi^^^'^^' ^'''^ 

where is the wavenumber energy spectrum of the turbulence, 7k is the linear 
growth rate of the mode with real frequency lJ]^, n is the number density, and 
Ue is the electron drift velocity. If the quantity k • Ue/wk is assumed constant 
and taken out of the integral sign in (19), then the ratio of heating of electrons 
to that of ions is 

m = Rfc) ■ 

This result is independent of the instability type. Galeev and Sagdeev (1983) 
point out that for most instabilities, relation (21) predicts faster heating of 
electrons than ions. Earlier, Lesch (1991) and Di Matteo (1998) emphasized the 
role of reconnection in accelerating electrons to relativisitic energies in accretion 
flows of AGNs. 



5 Accretion on to Magnetized Stars 

5.1 Disk Accretion on to a Rotating Star with an Aligned 
Dipole B-Field 

5.1.1 Introduction 

The problem of matter accretion on to magnetized stars has been of continued 
interest since the discovery of X-ray pulsars in binary systems (Giacconi et al. 
1971; Schrcicr et al. 1972; Tananbaum et al. 1972). The X-ray pulsars were 
interpreted as magnetized neutron stars, accreting matter from a companion 
star (Pringlc and Rees 1972; Davidson and Ostriker 1973; Lamb, Pethick, and 
Pines 1973; Rappaport and Joss 1977). Of the more than 40 of known X-ray 
pulsars, most are in high-mass {Mtot > 15) binary systems, while only few are in 
low-mass {Mtot < 3) binary systems (see, for example, review by Nagase 1989). 

From the early observations, many X-ray pulsars were found to show secular 
decreases of pulse period (spin-up of the neutron star rotation) with relatively 
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large rates of change of —PjP = 10~^ — 10~^ yr"^ (Gursky and Schreier 1975; 
Schreier and Fabbiano 1976; Rappaport and Joss 1983). The long-term moni- 
toring of X-ray pulsars during the last two decades have shown a wide variety 
of pulse-period evolution on different time-scales. Some X-ray pulsars show a 
steady spin-up evolution during many years which rapidly changes to a steady 
spin-down evolution (e.g., GX H-4, SMC X-1, 4U 1626-67). Others show the 
wavy variations of P on time-scales of a few years on a background of systematic 
spin-up or spin-down (e.g., Cen X-3, Her X-1, Vela X-1) (Nagase 1989; Sheffer 
et al. 1992; Bildsten et al, 1997). 

Short-term fluctuations of pulse-period on time-scales of days to months, 
including clear evidence of spin-down episodes, were found for a number of 
sources, for example. Her X-1 (Giacconi 1974), Cen X-3 (Fabbiano and Schreier 
1977), Vela X-1 (Nagase et al. 1984). Detailed observations of Vela X-1 with 
Hakucho and Tenma satellites, and with HEAO-1 have revealed fluctuations of 
large amplitude with time-scales about 3-10 days, which were estimated to be 
random in both amplitude and sign (Nagase 1981). BATSE observations have 
shown that in some cases (Cen X-3, Her X-1) the timescale of fluctuations is 
less than one day with the change of sign from spin-up to spin-down and vice- 
versa (Bildsten et al. 1997). Note that some X-ray pulsars (usually the long- 
period ones) show no evidence of regular spin-up or spin-down (for example, 4U 
0115-^63, GX 301-2). 

In some cases, like Her X-1, the existence of an accretion disk is indicated 
from various observations in the X-ray and optical bands (for example, Middled- 
itch and Nelson 1976; Bisnovatyi-Kogan et al. 1977; Middleditch, Puetter, and 
Pennypacker 1985). In several other cases, such as Cen X-3, SMC X-1, GX 1-1-4 
(Chakrabarty et al. 1997; Cui 1997), and 4U 1626-67, there arc also different 
kinds of evidence, that matter from companion star forms a disk (for example, 
Tjemkes et al. 1986; Nagase 1989). Most of these are short-period pulsars with 
period about several seconds or less (excluding GX 1+4 with period 122 s.). 
However, most of long-period pulsars, are supposed to be powered by the cap- 
ture of stellar wind from companion. However if the companion is the star of 
Be type with slow rate of outflow, then it is probable that once again the disk 
first forms, opposite to the case when the companion star is of 0-B kind with 
high speed of matter outflow (sec review of Nagase 1989). 

Some aspects of matter accretion by magnetized neutron stars were con- 
sidered for the first time by Bisnovatyi-Kogan and Friedman (1969), and by 
Shvartsman (1971). Ideas on the spinning-up of pulsars by accretion were first 
proposed in the case of disk-fed pulsars by Pringlc and Rces (1972), Lamb, 
Pethick and Pines (1973), and Lynden-Bell and Pringle (1974), and for wind- 
fed pulsars by Davidson and Ostriker (1973), and lUarionov and Sunyaev (1975) 
(see also Lipunov 1993). The spin-down of a rapidly rotating magnetized star 
due to "propellar" action of the magnetosphere was suggested by Shvartsman 
(1971) and Illarinov and Sunyaev (1975), but a detailed model has not been 
worked out. Pringle and Rees (1972) considered matter accretion from a vis- 
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cous disk to a magnetized rotating neutron star. Their theory considered the 
ideaUzed case of an aligned rotator where the magnetic moment of pulsar is 
aligned parallel or anti-parallel with the pulsar's spin-axis and the normal to 
the plane of the accretion disk. They supposed that accreting matter stops at 
the point where the magnetic pressure of the magnetosphere becomes equal to 
pressure of matter. At this distance, matter transfers its angular momentum to 
the star through a thin boundary layer. Another idea regarding the width of 
the transition zone between the unperturbed accretion disk and magnetosphere, 
was considered in series of papers by Ghosh and Lamb (1978, 1979a, b). They 
assumed that magnetic field of the neutron star threads the accretion disk at 
different radii in a broad "transition" zone, in spite of predominance of the disk 
stress compared with magnetic stress in this region. Their work is based on the 
idea of anomalous resistivity connected with the supposed fast reconnection of 
toroidal component of the magnetic field lines across the disk. 

Lovelace, Romanova, and Bisnovatyi-Kogan (1995, hereafter LRBK95; Li 
and Wickramasinghe 1997) proposed an essentially different picture for disk 
accretion onto a star with an aligned dipole magnetic field. The idea is based on 
the fact that when there is a large difference between the angular velocity of the 
star and that of the disk, the magnetic field lines threading the star and the disk 
undergo a rapid inflation so that the field becomes open with separate regions of 
field lines extending outward from both the star and the disk. As a result, the 
magnetosphere consists of an open field line region far from the star and a closed 
region approximately corotating close to the star. This is shown schematically 
in Figure 9. The LRBK95 model does not assume an anomalous resistivity, 
but rather a turbulent magnetic diffusivity proposed by Bisnovatyi-Kogan and 
Ruzmaikin (1976), analogous to the turbulent a viscosity of Shakura (1972) 
and Shakura and Sunyaev (1973). This corresponds to a resistivity about 10^ 
times smaller than that assumed by Ghosh and Lamb (1978, 1979a, b). Campbell 

(1992) has earlier discussed closed field line models of disk accretion onto an 
aligned dipole star assuming a magnetic diffusivity comparable to the turbulent 
a viscosity. 

The existence of an open magnetic field region of the disk leads to the possi- 
bility of magnetically driven outflows. LRBK95 developed a model for magne- 
tohydrodynamic (MHD) outflows and their back influence on the disk using the 
work on MHD outflows and magnetized disks by Blandford and Payne (1982), 
Lovelace. Berk, and Contopoulos (1991), Lovelace, Romanova, and Contopoulos 

(1993) , and Lovelace, Romanova, and Newman (1994). In §5.1.1 -5.1.8 we review 
the arguments of LRBK95. In §5.2 we discuss the magnetic propeller regime of 
accretion first considered by lUarionov and Sunyaev (1975) and more recently 
by Lovelace, Romanova and Bisnovatyi-Kogan (1999; hereafter LRBK99). 
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5.1.2 Theory 

The basic equations for an assumed stationary plasma configuration are 

47r 

V • (pv) = , V X B = —J , V-B = 0, VxE = 0, 

c 

J = cre(E + V X B/c) , p(v V)v = -Vp + pg+ ij X B + F™ . (299) 

c 

Here, v is the flow velocity, p is the density, Ce is the effective electrical conduc- 
tivity, F^"* is the viscous force density, p = pkBT/m is the gas pressure (with ks 
the Boltzmann constant, m the mean particle mass, and with the radiation pres- 
sure assumed negligible) , and g is the gravitational acceleration. Outside of the 
disk, dissipative effects are considered to be negligible (ere oo, F^"' = 0, etc.). 
We neglect the self-gravity of the disk and relativistic effects so that g = — V$g 
with $g = — G'M/(r^ -I- z^)^/^, where M is the mass of the central star. Equa- 
tions ( |299| ) are supplemented later by an equation for the conservation of energy 
in the disk. 

A general axisymmetric B-field can be written as B = Bp + where 
Bp = V X ((/()\['/r) is the poloidal field, is the toroidal field, and \1/ is the 
flux function (see, for example, Mestel 1968, or Lovelace et al. 1986). We use 
a non-rotating cylindrical coordinate system so that Bp = {Br,0, Bz). Notice 
that 4'(r, z) = const labels a poloidal field line, (Bp • V)\E' = 0, or a flux-surface 
if the poloidal field line is rotated about the z-axis. For the present problem, the 
Bp field can be represented as ^' = + 5*', with the star's field-assumed 
to be a dipole = pr^ir'^ + z'^)~i and with due to the non-stellar toroidal 
currents. 

5.1.3 Inflation and Opening of Coronal B-Field 

Figure 10 shows the poloidal projections of two nearby field lines connecting 
the star and the disk. Because d/dt = 0, the E field is electrostatic and the 
poloidal plane line integral of E around the loop, 1^2^3^4^ 1, is 
zero. Because of axisymmetry and the fact that E-|-vxB/c = outside of 
the disk, the line integrals of E along the two curved segments in Figure 10 
vanish separately. That is, the electrostatic potential is a constant on any given 
flux surface (see, for example, Lovelace et al. 1986). The potential difference 
between the points (1,2) on the star's surface is — (Sri2 • E* — w^^^'/c, with 
S'^ = 5*1 — 5*2, where we have assumed that the star is perfectly conducting 
[E* = — (v X B)*/c], and where is the angular rotation rate of the star. Thus, 
we have 5r^4^{Er)d = —ijJ*5'^/c — —Lu^,rBz{r,0)Sr34/c, where the path element 
6r34 is in the mid plane of the disk, and the "d" subscript indicates that the 
quantity is evaluated in the disk. In turn, we have Er\^ = {Jr/<^e — v<j>Bz/c)d- 
From Ampere's law, Jr = —{c/4:TT){dB^/dz), the fact that B^ is necessarily 
an odd function of z, and the approximation dB^/dz — (B^)h/h, we have 
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Jra « ~{c/4:TT)(Bcj,)h/h, where h is the half-thickness of the disk, and the 
notation (• • ■)h indicates that the quantity is evaluated dX z = h. Combining 
these results gives for stationary conditions 

hv 

{B^)h{r) = [uJd{r) - io,]B,{r) , (300) 

m 

where rjt = c^/(47rtTe) is the magnetic diffusivity of the disk, and oJd is the 
angular rotation rate of the disk. The fractional variation in Bz between the 
midplane and surface of the disk is negligible [0{h/r)] if the disk is thin in the 
sense that h/r <^ 1. For the conditions of interest here, it follows from ( [70| ) 
that h/r < Cs/vk^ where Cg = {ksT /m)^ the Newtonian sound speed based 
on the mid-plane temperature of the disk, and vk = {GM/r)^ is the Keplerian 
velocity at r. Campbell (1992) has independently derived equation ( [300D but 
does not place any limit on the ratio \B^/Bz\. 

In contrast with the work of Ghosh and Lamb (1978, 1979a, b), we assume 
that the magnetic diffusivity of the disk rjt is of the order of magnitude of the 
disk's effective viscosity (see Bisnovatyi-Kogan and Ruzmaikin 1976; Parker 
1979; Campbell 1992). Additionally, we assume that the disk viscosity is the 
turbulent viscosity as formulated by Shakura (1972) and Shakura and Sunyaev 
(1973), vt ~ {'^/'i)oLCsh^ where a is a dimensionless quantity less than unity. 
In this work, we assume that a is in the range 10~^ to 10"^. The possible 
contribution to the turbulent momentum flux due to small-scale magnetic field 
fluctuations is assumed included in a (Eardley and Lightman 1975; Coroniti 
1981; Balbus and Hawley 1992; Kaisig, Tajima, and Lovelace 1992; Rastatter, 
L. & Schindler, K. 1999). We write rjt = Dvt where D is dimensionless and of 
order unity. 



With rjt the turbulent diffusivity, equation (|300|) gives \{B,f,)h/ Bz\ 
= {i/2)r\u)d{r) — |/(a£'cs). This ratio is a measure of the twist of the B field 
between the star and the disk. For the outer part of the disk u)d{r) is expected 
to be close to the Keplerian value lok — {GM/r^)^ . Thus, it would appear that 
the twist \{B^)h/Bz\ can be very much larger than unity, say, > 10. We argue 
here that such large values of the twist do not occur. 

For a discussion of the limitation on the twist \{B,j,)h/ Bz\ we consider that 
the plasma outside of the disk is force-free, J x B = 0, which is the coronal 
plasma limit of Gold and Hoyle (1960). This limit is valid under conditions 
where the kinetic energy density of the plasma is much less than B^/87r. The 
general response of force-free coronal magnetic field loops to stress (differential 
twisting) applied to the loop foot points (in the solar photosphere) has been 
studied intensely by Aly (1984, 1991), Sturrock (1991), and Porter, Khmchuk, 
and Sturrock (1992). The conclusion of these studies is that a closed field loop 
with a small twist evolves into an open field-line configuration as the twist is 
increased. The related problem of the twisting of a force free magnetic field 
configuration with foot points at different radii in a Keplerian disk has been 
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studied by Newman, Newman and Lovelace (1992), and Lynden-Bell and Boily 
(1994). The twisting due to the differential rotation of the disk acts to increase 
the total magnetic energy of the coronal field and this in turn acts to "inflate" 
the field. 

An alternative measure of the field twist between the star and the disk is 
simply the difference in the azimuthal location, A0, of the stellar and disk foot 
points of a given flux tube. We have rd(/)/dlp — Bif,{r, z)/|Bp(r, z)|, where dip is 
the length element along a poloidal field line ^'(r, z) — const. Because rB^{r, z) 
depends only on ^' in a force-free plasma, we have 

A0-[r(i3^),],,,, r'-^ ■ (301) 

J star ^1 V I 

For small A(f> (say, < 1), the poloidal field is close to that of a dipole, and the 
integral gives Ai^ — Ct[{B^)ii/ Bz]d with C* = 8/15. For increasing A(f), the 
value of Ct increases because of the longer path length and the smaller value of 
jV^*!. The present situation is analogous to a current-carrying plasma column 
with a longitudinal B field. For the plasma column there is the well-known 
Kruskal-Shafranov stability condition (see, for example, Bateman 1980) on the 
twist of the field (at the column's surface) over the length of the column. 
This condition, A^ < 27r, is required for stability against symmetry changing 
kink perturbations. A related, non-axisymmetric instability may accompany 
the above mentioned inflation and opening of the star/disk field in response to 
increasing A0. 

LRBK95 argued that there is a definite upper limit on the twist Atp of any 
field lines connecting the star and the disk. This limit implies a corresponding 



limit on \{B^)h/ Bz\d in equation (30C). If A(f> is larger than this limit, then the 



field is assumed to be open as indicated in Figure 9. 
5.1.4 Equations for Disk 

The flow velocity in the disk is v(r) = —u{r)f + v^{r)(j), where u is the accretion 
speed. The surface density of the disk is a. The main equations are: 

1. Mass conservation: 

M = 2nr<Tu, = const. (302) 

We assume that the change of AI with r due to possible MHD outflows is small. 

2. Radial force balance: 



(Tvl _ GMa {BrB,)h ^ d 
r 27r dr 



h 



— dzp . (303) 



3. Angular momentum conservation: 



— {MF) = ~r^{B^B,)h , F = r^u:+ , (304) 

dr u dr 



74 



where lo = v,p/r. The term —r^{B^Bz)h represents the outflow of angular 
momentum from the ±z surfaces of the disk (that is, a torque on the disk). 
We assume that the angular momentum carried by the matter of the outflow 
is small compared with that carried by the field. If the B field at r is open 
as discussed in § 5.1.3, then the angular momentum outflow from the disk is 
carried to infinity. On the other hand, if the field is closed, then the angular 
momentum outflow from (or inflow to) the disk is carried by the coronal B field 
to (or from) the star. 
4. Energy conservation: 



h 

n4 



-h 

The first term is the viscous dissipation, the second is the Ohmic dissipation, 
while the third term is the power per unit area carried off by radiation (in the 
±z directions) from the disk which is assumed optically thick. Here, k is the 
opacity assumed due to electron scattering, a and ob = acjA, are the usual 
radiation constants, T is the midplane temperature of the disk, and Tg// is 
disk's effective surface temperature. 

5. Conservation of magnetic flux for a general time-dependent disk: 



5 . „ . d 



— (rB,) = — rB,u - ^{Br)h 
at or v h 

where the generally small radial diffusion of the B^ field has been neglected (see 
LRN). The first term inside the square brackets represents the advection of the 
field while the second term represents the diffusive drift. In a stationary state, 

Mr) ^i^^^. (306) 
Bz Vt 

6. Vertical hydrostatic equilibrium: The ^-component of the Navier-Stokes 



equation (299) gives the condition for vertical hydrostatic balance which can be 
written as 

where Cs is the Newtonian sound speed based on the midplane temperature of 
the disk, and b = r (-Br)^ + (-B^)^ / (Airavj^) (Wang, Lovelace, and Sulkanen 
1990). Radiation pressure is negligible for the conditions of interest. For b ^ 
'2cs/vk, this equation gives the well-known relation h/r = Cs/vk (Shakura 
and Sunyaev 1973), while for b » 2cs/vk it gives h/r — b^^{cs/vK)^ which 
is smaller than Cs/vk owing to the compressive effect of the magnetic field 
external to the disk (Wang et al. 1990) It is useful to write b — e /3'^, where e = 
rB2(r,0)/(27ro-w|f) and = [{Br)l + {B^)l]/{2Bl). In § 5.1.6, we discuss that 
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(5 = 0(1) in order to have outflows from the disk, and /3 < 1 with no outflows. 
Thus, the Alfven speed in the midplane of the disk is va — vk^^ (h/r)^ . We 
may term the magnetic field as weak if e < Cs/vk ~ h/r. In this hmit, va/cs ~ 
(er/h)^ < 1, and the magnetic compression of the disk is always small. In the 
opposite, strong field limit, e > Cs/vk- If at the same time /3 = 0(1), then the 
disk is magnetically compressed, and for e ^ Cs/vk, va/cs ~ 1//3. 

Note that with the magnetic field terms neglected these equations give the 
Shakura-Sunyaev (1973) solution for region "6". 



5.1.5 Magnetically Driven Outflows 

Consider the outer region of the disk where the value of the field twist Acf) given 



by equation (|30l| ) [or | {B(j,)h/ \ of equation (|300| )] is larger than a critical value. 
In this region the B field threading the disk is open as discussed in § 5.1.3. The 
angular momentum flux carried by the disk is MF. If the angular rotation of 
the disk is approximately Keplerian, cj(r) « a;;f(r), then the viscous transport 
contribution to F is {v / u){dijj / dr) « ~{2>/2)r'^u)[h/ {Prf)]- As discussed in 
detail below a necessary condition for MHD outflows is th at Pr — 0{1). Con- 
sequently, the viscous contribution to F in equation ( |304|) is smaller than the 
bulk transport term by the small factor h/r « 1, and F w bjr^ . 

In equation ( |304| ) we therefore have d{MF)/dr « Muixf/'^.- In equation 
( ^04| ) we let {B^)h = -fi^B^{r,Q). Studies of MHD outflows (Blandford and 
Payne 1982, LBC, LRC) indicate that (3^ ~ const — 0(1). As a result, equation 



( P04D implies that 

[B.(r,0)]„ = fc/r5/4, with k = [M(GAf) ^(Z/?^)]* , (308) 

where the w subscript indicates the wind region of the disk. Of course at a 
sufficiently small radius, denoted r^i , the outfiow ceases, and the dependence of 
Bz{r, 0) reverts to approximately the stellar dipole field, /i/r^. As a consistency 
condition we must have 



£ = {k|r''^t)/i^^|rL) 



< 1 . (309) 



The transition between the dipole and outflow field dependences is handled by 
letting Bz{r,0) = {^J'/r'^)g{r) + (/c/r^/^)[l — .g(r)], where, for example, gir) = 
{1 + exp[(r - r „,i)/ Ar] |~^ and /^rjr^i < 1. 

Equations ( |302| ) - (|307|) can be used to obtain the disk parameters in the 
region of outflow. As shown below the viscous dissipation is much larger than 



the Ohmic, and then Eq. ( |30^ ) gives T = \%\M'^{GMfl'^Kl{\2%'K'^aDI3lac)f^^ 
r^^l^. For illustrative values, M = \Mq, M = lO^'^g/s « 1.6 x IQ-^Mo/yr, a = 
0.1, Pr — 0.58, and D ~ l,we find at the representative distance r — lO^cm that 
T « 0.44 X lO^K, VK = 1.15 x lO^cm/s, = 7.6 x lO^cm/s, u = 0.88 x lO'^cm/s, 
a = 180g/cm^, and h/r ^ 6.6 x 10^'^. The ratio of the ohmic dissipation to 
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that due to viscosity is: (4/27)aD/3j? <C 1. The ratio of the viscous dissipation 
to the power output of the outflows (the ± z Poynting fluxes) (per unit area 
of the disk) is: {Si /2){h/ r){f3rD)~^ <C 1. Thus, most of the accretion power 
for r > ryji goes into the outflows. However, the fraction of the total accretion 
power in outflows is small if r^i ^ r*, where is the star's radius 10^ cm for 
a neutron star). Note that M — lO^^g/s and M — IMq correspond to a total 
accretion power or luminosity Lg = GMM/r^ « 1.33 x 10^^(erg/s)(10^cm/r«). 

Conservation of the poloidal flux threading the disk implies that there is 
an outer radius, denoted r^o, of the outflow from the disk. That is, we have 
/J"™" rdr{k/r^/'^) — JJ"™" rdr{fi/r^). So, the outer radius is r^o ~ ''toi(/f )~*^^, 
where / «lforl — £^1 and / = 3/4for£^l. For r > r^o, Bz{r, 0) rapidly 
approaches ze ro. 

Equation (307) gives the disk thickness. From equation (308), we have e = 
rBl/{2TTav\) — [u / Cs){cs/vk)/ {2(3(1,) ■ The magnetic compression of the disk is 
small if e < Cs/vk or equivalently if u/c^ < 1. Our numerical solutions have 
u/cs < 1 in the region of outflow. 



5.1.6 Necessary Condition for MHD Outflows 

A necessary condition on f3r for MHD outflows from the disk can be obtained by 
considering the net force J^p on a fluid particle in the direction of its poloidal mo- 
tion above the disk. Notice that for stationary flows, the poloidal flow velocity 
Vp is parallel to Bp owing to the assumed axisymmetry and perfect conductiv- 
ity. However, in general v x B 7^ so that the fluid particles do not move like 
"beads on a wire". Above, but close to the disk (/i < z <C r), we assume that 
the poloidal fleld lines are approximately straight. Thus, the poloidal position 
of a fluid particle above the disk is r = {vo + S sm{9))r + S'cos(6')z, where S 
is the distance along the path from the starting value So — h/cos{9) at a ra- 
dius To, tan(0) = {Br)h/Bz = (3r, and z = Scos{9). The effective potential is 
U{S) = -{l/2)uj^{ro + S sm{9))^ -GM/\r\ (LRN), and the force in the direction 
of the particle's poloidal motion is !Fp = —dU/dS. For distances z much less 
than the distance to the Alfvcn point, uJq ~ wk(''o)- In this way we find 

•^P = -K-^'>(^) +^1^(3/3.' -1)(^) , (310) 

for h < z r, where we assume B^ > 0. 

The slow magnetosonic point of the outflow occurs where J-p — (LRC) at 
the distance 

The factor 1 — {uj/luk)'^ can be obtained from the radial force balance equation 
(^031). The radial pressure force is small compared with the magnetic force for 
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the conditions we consider, a(3^/{3l3^) > h/r ^ 1. Thus, 



h ~\ 2,MujK / 3/32 - 1 



(312) 



For r ^i < r < Two, the quantity in brackets is simply a/ (S/?^) owing to equation 
( ^08|) . Note that the minimum of /3^/(3/32 - 1) is 1/2 at 13^ = 1. If the gas near 
the surface of the disk {h < z < z^) is assumed isothermal with temperature 
Ti <C T and sound speed Csi = (fcsTi/m) 2 , then the density at the slow magne- 
tosonic point ps can be obtained using the MHD form of Bernoulli's equation. At 
the slow magnetosonic point, the poloidal flow speed is Vp = Csi|Bp|/|B| = Vsm, 
the slow magnetosonic speed, and ps — phsxp { — 1/2 — [U{zs) — U{h)]/cli} , 
where U is now regarded as a function of z. In turn, the mass flux density from 
the +z surface of the disk is PsVsm cos{9) — psCsiBz/\^\- 

The consistency of the magnetically driven outflow solutions requires Zg > 
h. If this were not the case, the outflows would be matter rather than field 
dominated near the disk. There are two possible regimes having Zg > h: one 
is with aD/{3P^) < 1 [or aD/{3P^) <^ 1] and (3r > l/\/3; another is with 
aD/{Qj3^) > 1 and Pr > 1- In the present work, we consider the first regime 
which is consistent in the following respect: The magnetic pinching force on the 
inner part of the outflow (r < Vwo) increases as the mass flux in the outflow 
increases (LBC). In turn, an increase in the pinching force acts to decrease 
Pr = {Br)h/Bz thereby increasing z^/h and decreasing the mass flux in the 
outflow. 

In addition to the condition Zs > h, the slow magnetosonic point must be at 
a distance Zs not much larger than h, say 2h, in order for the outflow to have a 
non-zero mass flux. At the inner radius of the outflow, r^,,; , the Bz{r,0) held 
in the disk is /i/(2r^j) as discussed in § 5.1.5. We then have 



> 



24M(GAf)- 



(313) 



0.75 X iC^cm 



Oi 



1030Gcm'^ 



M 



\ M ) 



from equation ( |312| ). Note that from equation ( 309 ), we get E > [aZ)/(48/30)] 2 
01 £ > 0.046 and r^o < QOr^i for a = 0.1, D = 1 and /3(A = 1. 



5.1.7 MHD Outflows and Spin-Up {rto < Vcr) 



We have numerically integrated equations ( |302| ) - (307) assuming magnetically 
driven outflows from the disk. We integrate the equations inward starting from 
a large radial distance r < r^o ■ The inner radius of the region of outflow is 
determined as discussed above. For r < r^i , dF/dr = 0, and the solutions in 
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this range of r all exhibit a "turn-over radius", rto , where du/dr = 0. This 
"turn over" results from the radially outward magnetic force in equation ( |303[ ) 
which becomes stronger as r decreases. In the region close to the turnover, the 
magnetic compression of the disk becomes strong in that e = 0{1). From a least 
squares fitting of many integrations, we find 

no « 0.91 X 10^cm(^)° V ^ 

VO.l/ Vi03OGcm'V V Af y V M / ' 

The turn-over radius is less than but in all cases close to Vwi- Thus, equation 
( ^14| ) is compatible with ( [313| ). Our radius rto has a role similar to that of the 
Alfven radius rA of Ghosh and Lamb (1978, 1979a, b). The dependences we 
find of rto on /i, M, and M are close to those of r^. However, our analysis 
shows an important dependence on aD which is proportional to the magnetic 
diffusivity of the disk. Furthermore, our rto is smaller than by a factor of 
about {aD/l2)°-^, which is w 0.24 for aD = 0.1. 

The turn-over radius is important in the respect that the inward angular mo- 
mentum flux carried by the disk (MF) is Mcotorf^ because dcu/dr = at r = rto. 
[Note that we do not consider here the possibility of significant poloidal current 
flow along the open field lines extending from the polar caps of the star {B^ > 
for z > in figure 9) which would act to remove angular momentum from the 
star to infinity.] From numerical integrations, we find to a good approxima- 
tion that ojto — (GM/rfg)^ J so that the influx of angular momentum to the 
star is MFto — M{GMrto)2 , Thus, the rate of increase of the star's angular 
momentum J is: 

^=M{GMrto)i. (315) 

With / the moment of inertia of the star, J = loj^, and dJ/dt = {dl /dM)Mu!^ + 
I{dut,/dt). For the situation of interest, the term proportional to dl/dM is 
negligible because rto is much larger than the star's radius. Thus, we have 
"spin-up" of the star, 

duj^ _ M{GMrto)i 

~dr ~ I ' 

or 

P dt yrVlJVlOiVs^V / )\McJ VlO^cmJ ' 

. (^^^^ 

where P — I-kIlo.!, is the pulsar period. Because the total accretion luminosity 
L = GMM/r.^ with the star's radius, P cx _p2iO-85^o.285 f^j. constant 
aZ?, r*, and M . 

Our outflow solutions are consistent only in the case where ujto = {GM/rf^)^ 
> uj^,. In this case, with r decreasing from rto^ the rotation rate of the disk a;(r) 
decreases and approaches the rotation rate of the magnetosphere (w lj^) from 
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above where it matches onto the magnetospheric rotation through a radiaUy 
thin turbulent boundary layer. We find no consistent stationary solutions when 
uj(r) decreases below tj*. The condition ujto > is the same as 

/GMY ^ - f M \ ^ f P 



..<....^— J .1.5xl0«cm^— J ^-J , (317) 

where Vcr is commonly referred to as the "co-rotation radius" . Note that Vcr 
decreases during spin-up. Thus, r^^ approaches rto if 7'to does not decrease 
rapidly. 

5.1.8 Magnetic Braking of Star by Disk; Spin-Up or Spin-Down 

{rto > rcr) 

For Tto > Tcr, the B field in the outer part of the disk (r ^ Vcr) is open as 
discussed above, but there are no consistent stationary MHD outflows, {B^)fi = 
0, f3r = hu/r]t <^ 1, and the Ohmic dissipation is small compared with that due 
to viscosity. Thus, the outer part of the disk obeys essentially the equations 
of Shakura and Sunyaev (1973) for a Keplerian disk with F ~ const. = Fqo in 



equation (304) undetermined. We show below that Foo can be determined self- 
consistently by considering the region near Vcr where the Keplerian accretion 
flow is brought into co-rotation with the star. Because there are no outflows, 
the angular momentum influx to the star is MF^o which is equal to dJ^/dt. 
In the region of the disk where w(r) is close to lo^, that is, where r ~ Tc,-, 



\{B(j,)h/ Bz\ given by equation (300) is not much larger than miity. In this region, 
closed but twisted held lines link the star and the disk. The twist of the field 
acts to remove (or add) angular momentum from (to) the disk if B^B^ > (or 
< 0). The angular momentum removed from the disk is deposited on the star 



via the B field. In this region of the disk the key equations are (300) and ( pO 
which we rewrite as 

— = — (F-wr2), = {uj-uj,)[—rf\H[ . (318) 

Here, H{x) is a Heaviside function such that H{x) = 1 for \x\ < 1 and H{x) — 
for |a;| > 1; 

/ N _ {B4>)h{r) hr 

is a measure of the field twist; T„iax — const is the maximum value of the twist; 
and 

X 9 r^Vt diu 

F{r) = r^uj + -— . 

u dr 

We consider r^ax to be a universal constant. For |t| > Tmax the field configu- 



ration becomes open. In equations (316), B^ is assumed to be the star's dipole 



field. The other equations of §4.1.4 are still needed. Note that both the viscous 
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and Ohmic heating must be retained in the energy equation (305). The viscous 



dissipation is dominant for r > Vcr, while the Ohmic dissipation dominates for 



r < Tcr- Equations ( p03D and ( p06D can be combined to give u as a function 
of r, T, and w. In turn, this expression for u can be combined with the en- 
ergy conservation equation ( |305| ) to give both u and T as functions of r, and 
duj/dr if viscous heating dominates, or as functions of r and lo if Ohmic heating 
dominates. Furthermore, Vt = {2/3)aCsh in equation ( |31q) ca n be derived from 



T and h/r from equation (307) using /3j. from equation ( |30q ), /J^ = — r, and 
from Bz and u. 



We have solved the equations (|o|, [305|^ , & ITs]) by numerical integration 



oo : 



starting from the outside, at a distance where r > Tmax and F =const= 
and integrating inward through Vcr- The fact that the value of i^oo is not known 
a priori points to the use of a "shooting method" for its determination. Using 
this approach, we find that there is in general a unique value of Foo, denoted 
F^, such that the solution for a;(r) smoothly approaches for r decreasing 
below Tcr- If Foo is smaller than F^, then u;(r) follows the Keplerian law as r 
decreases below Vcr- On the other hand, if Foo is larger than F^, then u!{r) goes 
through a maximum and decreases rapidly to values much less than . The 
only physical solution is that with _Foo = F^. Hereafter, the zero superscript 
on Foo is implicit. 



The behavior of equations (318) can be understood qualitatively by noting 
that if i^{r) were to decrease linearly through Vcr, then the radial width of the 
region of braking, dF/dr < where u; < lu^,, would equal that the width of the 
region of acceleration, dF/dr > where w > w*. Consequently, F would change 
by only a small fractional amount. However, if within the braking region F{r) 
attains a value equal to ojr^, then duj/dr — at this point and dtu/dr remains 
small and www* for smaller r. A rough estimation gives the radial width of 
the braking region as Ar TmaxDcg/uJK, where Cs is based on the temperature 
just outside of this region, and the jump in F as AF ~ Tmax^f'{r'^ Bl / M) . 

Figure 1 1 shows the dependence of Foo on /i and on M for sequences of cases 
where rcr is assumed comparable to rto- This is clearly not a necessary choice; 
that is, Tto could differ significantly from r^r- However, with r^-r = Tto, Figure 11 
allows a direct comparison of the angular momentum influx to the star (MFk) 
for the case of outflows {rto < '"cr) with the influx (A/Foo) in the case of no 
outflows and magnetic braking {rto > Tcr)- Figure 11 suggests that a bimodal 
behavior can occur where the star switches between spin-up and spin-down, 
owing to some extraneous perturbation. 

5.2 Magnetic Propeller Regime 
5.2.1 Introduction 

Observations of some X-ray pulsars show remarkable 'jumps' between states 
where the pulsar is spin-ning-down to one where it is spinning-up. Examples 
include the objects Cen X-3 (Bildsten et al. 1997) and GX 1-1-4 (Chakrabarty 
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et al. 1997; Cui 1997). The theoretical problem of disk accretion on to a ro- 
tating magnetized star has been discussed in many works over a long period 
(Pringle & Rccs 1972; Lyndcn-BcU & Pringle 1974; Ghosh & Lamb 1979a,b; 
Wang 1979; Lipunov 1993; Lovelace, Romanova, & Bisnovatyi-Kogan 1995 (re- 
ferred to as LRBK95); Li & Wickramasinghe 1997). However, except for the 
work by Li & Wickramasinghe (1997), the studies do not specifically address 
the "propeller" regime (lUarionov & Sunyaev 1975) where the rapid rotation of 
the star's magnetosphere acts to expell most of the accreting matter and where 
the star spins-down. Recent computer simulation studies of disk accretion on 
to a rotating star with an aligned dipole magnetic field (Hayashi, Shibata, & 
Matsumoto 1996; Goodson, Winglee, & Bohm 1997; Miller & Stone 1997) pro- 
vide evidence of time-dependent outflows but do not give definite evidence for 
a "propeller" regime with spin-down of the star. Recent computer simulation 
studies of disk accretion on to a rotating star with an aligned dipole magnetic 
field (Hayashi, Shibata, & Matsumoto 1996; Goodson, Winglee, & Bohm 1997; 
Miller & Stone 1997) provide evidence of time-dependent outflows but do not 
give definite evidence for a "propeller" regime with spin-down of the star. 

Lovelace, Romanova, and Bisnovatyi-Kogan (1999, hereafter LRBK99) de- 
veloped a model for magnetic "propeller" -driven outflows which were first con- 
sidered by lUarionov and Sunyaev (1975). Such outflows can cause a rapidly 
rotating magnetized star accreting from a disk to spin-down. Energy and an- 
gular momentum lost by the star goes into expelling most of the accreting disk 
matter. The envisioned geometry is shown in Figure 9. The theory gives an 
expression for the effective Alfvcn radius Ra (where the inflowing matter is ef- 
fectively stopped) which depends on the mass accretion rate, the star's mass and 
magnetic moment, and the star's rotation rate. The model points to a mecha- 
nism for 'jumps' between spin-down and spin-up evolution and for the reverse 
transition, which are changes between two possible equilibrium configurations 
of the system. In for example the transistion from spin-down to spin-up states 
the Alfvcn radius Ra decreases from a value larger than the corotation radius 
to one which is smaller. In this transistion the "propeller" goes from being "on" 
to being "off." The ratio of the spin-down to spin-up torque (or the ratio for 
the reverse change) in a jump is shown to be of order unity. 

5.2.2 Theory 

Consider disk accretion on to a rapidly rotating star so that the disk-star con- 
figuration is as sketched in Figure 12. Consider the flow of mass, angular mo- 
mentum, and energy into and out of the annular region indicated by the box 
A'ACC in this figure, where A' A is at radius ri and CC is at r2- Notice that 
at this point the values of ri and r2 are unknown. They are determined by the 
physical considerations discussed here. 

Consider first the outer surface CC through the disk. The influx of mass 
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into the considered region is 

Maccr = 27rr2 / dz [p {-Vr)]2 > , (319) 

J-h 

where h is the half-thickness of the disk, and the 2 subscript indicates evaluation 
at r = r2. We assume that mass accretion rate for r > r2 is approximately 
constant equal to Maccr- That is, wc consider that outflow from the disk is 
negligible for r > r2- The influx of angular momentum into the considered 
region is 

Laccr = M,,,r{rv^)2 - T2 , T2 = 2^^ / dz {T^l')2 , (320) 

J — h 

where T"^* is the viscous contribution to the stress tensor which includes both 
the turbulent hydrodynamic and turbulent magnetic stresses. The influx of 
energy into the considered region is 

Eaccr = Maccr + " ^S^s , (321) 

where uj = v^/r vk/t with vk = \/GM/r the Kcplcrian speed and M the 
mass of the star, and where W2 is the enthalpy. For conditions of interest here 
the disk at r2 is geometrically thin so that w <C where Cg is the sound 
speed. 

Consider next the fluxes of mass, angular momentum, and energy across the 
surface AA' in Figure 12. For the physical regime considered, where (ri)/ri 
<C {— the angular rotation rate of the star), the mass accretion across the 
A A' surface is assumed to be small compared with Maccr- The reason for 
this is that any plasma which crosses the AA' surface will be 'spun-up' to an 
angular velocity (by the magnetic force) which is substantially larger than 
the Kcplcrian value, and thus it will be thrown outwards. Thus the efflux of 
angular momentum across this surface from the considered region is ii = — Ti . 
The efflux of energy across this surface is Ei = —lu^Ti, where uj^ is the angular 
rotation rate of the star and the inner magnetosphere as shown in Figure 12. 
For the conditions considered, the star slows down and loses rotational energy 
so that Ti > 0. 

We have Ei/Li ^ uj^,. Because the interaction of the star with the accretion 
flow is by assumption entirely across the surface AA' , this is consistent with the 
spin-down of a star with constant moment of inertia 7; that is, AE^/AL^ = 
/»a;*Au;*/(/,Au;*) = w*. 

Next we consider the mass, angular momentum, and energy fluxes across the 
surfaces A'C and AC in Figure 12. As mentioned, accretion on to the star is 
small for rcr ri where rcr is the corotation radius as indicated in Figure 12. 
Thus, the mass accretion goes mainly into outflows, Mo«t ~ Maccr, where "ouf 
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stands for outflows. The angular momentum outflow across the surfaces A'C 
and AC, Lout, must be the difference between the angular momentum lost by 
the star and the incoming angular momentum of the accretion flow. The angular 
momentum carried by radiation from the disk is negligible because (vk /c)^ <C 1, 
where c is the speed of light. That is, Lout = Laccr — Li. The energy outflow 
across the A'C and AC surfaces is Eout + E^ad — Eaccr — Ei, where Erad is 
the radiation energy loss rate from the disk surfaces between r = ri and r2 , and 
Eout is the rate of energy loss carried by the outflows. 
Angular momentum conservation gives 



Lout = Maccr{rV4,)2 - T2 + Ti . (322) 



GMMaccr 

Eout + Erad = 7: i^2T2 + UJ*Ti . (323) 

2^2 



We have 



We can solve equation (p23) for T2 and thereby eliminate this quantity from the 



energy equation (323). This gives 



Eout + Erad — ^ + ^2Lout + {^* — ^2)Tl . (324) 

2r2 

The preceeding equations are independent of the nature of the outflows from 
the disk. At this point we consider the case of magnetically driven outflows 
as treated by Lovelace, Berk, and Contopoulos (1991, hereafter LBC). In the 
LBC model the outflows come predominantly from an annular inner region 
of the disk of radius ~ Ra where the disk rotation rate is ujq — \/ CM / R\ . 
Thus we assume that the outflows come from a region of the disk which is 
approximately in Keplerian rotation. For the present situation, shown in Figure 
1, it is clear that we must have ri < Ra < f2 ■ Further, we will assume 
Ra is close in value to ri with {Ra — ri)/i?^^ 1. For conditions where the 
outflow from the disk is relatively low temperature (sound speed much less than 
Keplerian speed), equations (16) and (18) of LBC imply the general relation 
Eout = Wo Lout — iCM Mout I {"^ Ra ) ■ This equation can be used to eliminate 



Lout in favor of Eout in equation ( 324 ). Recalling that Mout — Maccr we have 



1 - 8-^)Eout + Erad - C^* - W2 Ti ^^££L ^ _ 335 

2r2 

where 6 = i?yi/?'2 < 1- 

The energy dissipation in the region of the disk r = ri to r2 heats the disk 
and this heat energy is transformed into outgoing radiation Erad- Thus we have 

Erad = dr , (326) 
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(Shakura 1973; Shakura & Sunyaev 1973), where Mdisk{r2) = Maccr and Mdisk 
(ri) « 0. The essential change in Mdisk{r) occurs in the vicinity of Ra so that 
Erad ~ {3GMMaccr/2){R^^ - T^^). Thus equation ( p5| ) becomes 



iGMMg 



-(l-<5^) 



(327) 



Thus the power from the spin-down of the star o^^Ti must be larger than a 
certain value in order to drive the outflow. 

For magnetically driven outflows, the value of Eout can be written as 



3 2 1 ^ 1 • i 

Eout = 2 -^0 ^A ^0 Maccr , 



(328) 



(equation 34 of LBC), where J^o is a dimensionless numerical constant ^ 0.234, 
and Bq is the poloidal magnetic field at the base of the outflow at r = Ra- We 
take the simple estimate Bq = pi/ R\ , which omits corrections for example for 
compression of the star's field by the inflowing plasma. 

Next we consider the torque on the star Ti. Because most of the matter 
inflowing in the accretion disk at r — r2 is driven off in outflows at distances r > 
ri , the stress is necessarily due to the magnetic field. The magnetic field in the 
vicinity of ri has an essential time-dependence owing to the continual processes 
of stellar flux leaking outward into the disk, the resulting field loops being 
inflated by the differential rotation (LRBK95), and the reconnection between 
the open disk field and the closed stellar field loops. The time scale of these 
processes is ti ~ 27rri /w^ (j'l ) . We make the estimate of the torque Ti = 
—2TTrf{2Az){BrBij,)i/{4:TT), where Az is the vertical half-thickness of the region 
where the magnetic stress is significant, and where the angular brackets denote 
a time average of the field quantites at r '--^ ri. The magnetic field components, 



Br 



Bif,, with B^ cx —Br necessarily, must be of magnitude less than or of the 



order of the dipole field Bi = p/r\ at r = ri. The fact that B^ has the opposite 
sign to that of Br is due to the fact that the B^ field arises differential rotation 
between the region r < ri, which rotates at rate w*, and the region r > ri, 
which rotates at rate a;x(^i) < Also, it is reasonable to assume Az ~ ri. 
Therefore, 

,2 



^A 



(329) 



where a < 1 (the time average of a{t)) is a dimensionless constant analogous to 
the a— parameter of Shakura (1973) and Shakura & Sunyaev (1973). (Note that 

because ri ^ Ra, a and a\ — (r\j RA)'^a are of the same order of magnitude.) 
Substituting equations (|32§|) and (p29^ into equation (327) gives 



Ra 



--{1-S^). (330) 
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Here, we have introduced two characteristic radii - the corotation radius, 

1/3 



( GM\ g ,,,l/3n2/3 

= 1 — y\ w 1.5 X 10 cm M/ P^' , 



with Pi = (27r/a;*)/ls the pulsar period and M\ = M/MQ. [For a young stellar 
object, rcr ~ 1-36 x lO^^M^'^F^p^^, where Piod is the period in units of 10 days.] 
The second is the nominal Alfven radius 



Ra 



rVGM 



2/7 4/7 

w 3.6 X lO^cm . ^ ,^ , (330a) 



where Mn = Maccr/(10^^g/s) is the accretion rate with lO^^g/s w 1.6 x 
10~^ Mq/ut and fi^o = /i/(10"^°Gcni^) with the magnetic field at the star's 
equatorial surface ^/r^ = 10^^G'(10^cm/r)'^. [For a young stellar object ta ~ 
1.81 X lO^^cm/igg^g /(M^/''My^), where the normalization corresponds to a 
stellar radius of 10^^ cm, a surface magnetic field of 3 x 10"^ G, and an accre- 
tion rate of 1.6 x lO~^M0/yr.] The corotation radius is the distance from the 
star where the centrifugual force on a particle corotating with the star (w^r) 
balances the gravitational attraction (GAZ/r^). The Alfven radius is the dis- 
tance from a non-rotating star where the free-fall of a quasi-spherical accretion 
flow is stopped, which occurs (approximately) where the kinetic energy-density 
of the flow equals the energy-density of the star's dipole field. Note that the 
assumptions leading to equation (330a) require Ra > rcr- 

Notice that Ra (or ri ~ Ra) is the effective Alfven radius for a rotating 
star. It depends on both va and Vcr in contrast with the common notion that 



the Alfven radius is given by even for a rotating star. From equation (329), 
the spin-down rate of the star is I{dLLiifldt) = — a fi^ /R\, where / is the moment 
of inertia of the star (assumed constant). Thus the spin-down rate depends on 
both VA and rcr- 

Figure 13 shows the dependence of Ra on va and rcr for a sample case. 
For conditions of a newly formed disk around a young pulsar, the initial system 
point would be on the upper left-hand part of the curve. Due to the pulsar 
slowing down (assuming fi and M constant), the system point would move 
downward and to the right as indicated by the arrow. In this region of the 
diagram, Ra « rH^/rll^ oc ^/uJ (for 6 <C 1), so that the torque on the star 
is Ti « (GMM)3/2/(^y^ u;3/2). Thus the braking index is n = -3/2, where n 
is defined by the relation u!■^, — —const tj". Numerically, 

1 dP _ 1.55 X 10-^ Pf/'(Mi Mi7)3/2 ^gg^^ 



P dt yr ^ L 



45 



where I45 = //(lO'^^g cm^). For Ra ^ fcr, the mass accretion rate to the star 
is small compared with Maccr, but some accretion may occur due to 'leakage' 
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of relatively low angular momentum plasma across field lines near ri (Arons & 
Lea 1976). 

As Ra decreases, the spin-down torque on the star increases. Over a long 
interval, Ra will decrease to a value larger than rcr but not much larger. In this 
limit, mass accretion on to the star Mi may become significant. Our treatment 
can be extended to this limit by noting that Mi = M2—Mout, ii = Miujirf—Ti, 



and El = Mi{-GM/2ri) - u^Ti, where Ti is given by equation ( |329D . In this 
limit the accretion luminosity is GMMi/r^, where is the star's radius. Figure 
13 is not changed appreciably for Mi < Maccr- 

Further spin-down of the star will cause the system point in Figure 13 to 
approach the right-most part of the curve. Further spin-down of the star is 
impossible. At this point of the evolution, the only possibility is a transition to 
the spin-up regime. In this regime the effective Alfven radius is the 'turnover 
radius' rto of the disk rotation curve calculated by LRBK99, the star spins-up 
at the rate Iduii,/dt = MaccriGMrtoY^"^ , and most of the disk accretion Maccr 
falls onto the star. The dashed horizontal line in Figure 13 indicates rto which 
is necessarily less than Tcr- 

The location of the turnover line rto in Figure 13 suggests the possible evo- 
lution shown by the sequence of points a— >6— s-c^d^a. The system can 
jump down from point a where the star spins-down to point h where it spins-up. 
The spin-down torque at a is —Q.fi^/R\^ (where RAa is the effective Alfven 
radius at point a), whereas the spin-up torque at b is Maccr{GMrtoY^'^ ■ The 
magnitude of the ratio of these torques is 

spin-down -fZA\"" (^\'\ (332) 



spin -up \RAaJ \ rto 

With the system on the rto hne, it evolves to the left. Because rto < ?'cr, there 
must be an upward jump from point c to point d. For this case the torque ratio 



is given by equation (332) with RAa — RacI- From point d the system evolves 
to the right. For the example shown in Figure 2, the torque ratio is ~ 3.45 for 
a ^ b whereas it is « 0.87 for c ^ d. The vertical line c ^ d is at the left-most 
position allowed for the considered conditions, but this transistion could also 
occur if the line is shift to the right. The line a b can be displaced slightly 
to the right or it can be displaced to the left to be coincident with the c ^ d 
line. In the latter case the torque ratio for the spin-down to spin-up jump is 
approximately equal to the torque ratio for the spin-up to spin-down jump and 
is « 0.87. The smaller the horizontal separation of the a — > 6 and the c d 
lines, the shorter is the time interval between jumps. 

Summarizing, we can say that the horizontal locations of the transistions, 
a ^ b, and c ^ d, are indeterminate within a definite range. The locations of 
the jumps in the {Ra-, rcr) plane may in fact be a stochastic or chaotic in nature 
and give rise to chaotic hysteresis in the of the spin-down/spin- up behavior of 
the pulsar. The jumps could be triggered by small variations in the accretion 
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flow {M for example) and magnetic field configuration (the time-dependence 
of a in the torque Ti). Analysis of the accreting neutron star system Her X- 
1 (Voges, Atmanspacher, & Scheingraber 1987; MorfiU et al. 1989) suggests 
that the intensity variations are described by a low dimensional deterministic 
chaotic model. The transistions between spin-down and spin-up and the reverse 
transistions may be described by an analogous model. 

The allowed values in Figure 13 have Tct/ta < const = fc, where k « cP'/^. 
This corresponds to pulsar periods 

/ - \ / 2 \ 3/7 
M30 



'^-'-''W)\^^> ■ ''''' 

For some long period pulsars such as GX 1-1-4 this inequality points to mag- 
netic moment va lues /H30 appreciably larger than unity. Periods much longer 
than allowed by (|333|) can result for pulsars which accrete from a stellar wind 



(Bisnovatyi-Kogan 1991). [For a young stellar object, equation ( |333[ ) gives 
P<8d(a/Mf/^)(^i6.5/Mi8)3/7.] 

5.3 Cyclotron emission from magnetic poles in X-ray pul- 
sars 

Accretion flow into a strongly magnetized neutron star is canalyzed by the 
magnetic fleld and is falling into magnetic poles. Observational spectra of X- 
ray pulsars are formed by combination of black body and comptonization in 
the hot coronal layer. In presence of a strong magnetic fields ^ 10^° — 10^^ 
Gs the cyclotron or synchrotron emission of electrons also falls into an X-ray 
band. There are a lot of observations of the spectral feature, interpreted as a 
cyclotron line in the X-ray spectrum of the pulsar Her X-1 (Triimper et al., 
1978; Tuellcr ct al., 1984; Ubcrtini ct al., 1980; Voges et al., 1982; Gruber et 
al., 1980; McCray at al., 1982; Mihara et al., 1990; Sheepmaker et al. 1981) at 
39-58 KeV (see Table Magnetic field intensity was usually calculated from 
the non-relativistic formula 

H= , (334) 

e 

where u> is the cyclic frequency of the photons, m is the mass of the electron. 
From observations of the line frequency in the Table the magnetic field strengths 
are of the order of (3 — 5) • 10^^ Gs. But as large as this value comes into con- 
flict with some theoretical reasonings like interpretation of the observations of 
pulsar spin acceleration (Bisnovatyi-Kogan and Komberg, 1973), condition for 
the transparency for the outgoing of the directed radiation (Bisnovatyi-Kogan, 
1973, 1974), consideration of the interrelation between radio and X-ray pul- 
sars (Bisnovatyi-Kogan and Komberg, 1974), interpretation of beam variability 
during the 35-day cycle (Sheffer et.al., 1992; Deeter et al., 1998; Scott et al.. 
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2000), see also Lipunov (1987). It seems likely that the reason of this conflict 
is an unsuitability of the non-relativistic formula in this case. According to 
Bisnovatyi-Kogan and Fridman (1969), the temperature of the electrons emit- 
ting a cyclotron line could be ~ lO^^K, and therefore they are ultrarelativistic. 
By this means the mean energy of the cyclotron line is broadened and shifted 
relativistically by a factor of 7 ~ . In the paper of Baushev and Bisnovatyi- 
Kogan (1999) the spectral profile of the cyclotron line had been calculated for 
various electron distributions. Furthermore, the model of the hot spot on the 
pulsar magnetic pole was considered and it was shown that the overall observed 
X-ray spectrum (from 0.2 to 120 KeV) could arise in presence of the surface 
magn etic field (~ 5 • 10^" Gs) which is well below then those, obtaind from 

(il- 

5.3.1 The spectrum of the cyclotron radiation of the anisotropic 
relativistic electrons 

According to Bisnovatyi-Kogan (1973), Gnedin and Sunyaev (1973), in the mag- 
netic field near the pulsar the component of a momentum perpendicular to the 
magnetic field emits rapidly, while the parallel component of a velocity remains 
constant. Hence the momentum distribution of the electrons is anisotropic 

(335) 

where p± <^ mc, 3> mc. Assume for simplicity that the transverse electron 
distribution is two-dimensional Maxwellian 

/ mu'^ \ mv? 
dn=-exp(^-— jd— , (336) 

where Ti <C mc? . Let us calculate the cyclotron emission of N such particles 
that move at a rate V along the magnetic field. 

For a single particle, having the transverse velocity m, we find (Ginzburg, 
1975): 



e^H^u^l - ^)^[(l + cos0)(l + ^) - cos ( 
87rc^m2(l — cos( 

where 9 is an observational angle in a laboratory frame of reference. Integrating 



m - ^ ""l^ . Z V Zn.t (337) 



over the distribution (336), we obtain for A^ particles: 



^e^H^T,{l ~ y^f[{l + cos0)(l + K)-A 



V 



cos f 



J(6») = / j{e)dn = N ^ ^ £ -. (338) 

For the spectrum we find: 
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1 _ 

u;{e)=ujH- y 

1 — cos t 



LOH 



mc 



(339) 



When y ~ c the cyclotron radiation is highly directed and diagram has a pencil 
beam along V (0 ~ 0). Under these conditions {6 = 0, V ~ c) we obtain from 
(p8|),(|339|): 



what gives 



J(0) 



W(0) = LOH] 



1 



'1+1 
1-^ 



F 2lj 



2u!h- 



Eu 



TUeC 



H 



C OJ^ 

Let us consider the parallel momentum distribution of the electrons as: 

dn = f{p\\)dp\\. 
Substituting of dn for N and using 



mc uj 
2 UH 

we obtain for the spectral density: 



1 



V 2lu 



H 



Let us consider two important cases. When / is a relativistic Maxwell: 
noc 



f 



■ exp 



T2>mc^>Ti, 



P||C 

where no is a number of emitting electrons. Then the spectrum is: 



Lo exp 



2ttcloh T2 

This spectrum has a single maximum at 



2ujhT2 



dLO. 



4T2 

r.2 • 



LOH TIIC^ 

In the second case consider the function / as 



(340) 
(341) 

(342) 
(343) 
(344) 

(345) 

(346) 

(347) 

(348) 
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The spectrum of radiation is 

2 / / mc til n'lZ \ 

= J '^"^ 

When (7 ^ a this spectrum has a single maximum at 

2a 

CO ~ — UJH- (351) 
mc 

Notice that in all cases the maximum is shifted to 

^ ^ A. (352) 

It is a common property of relativistic cyclotron line, that is independent of 
the particular form of /. We had approximated experimental spectrum taken 



from Mihara et al. (1990) (solid line) by (347), (350), representing the second 
variant with better coinsidence by the dot line in Figure 14. The spectrum from 
McCray et al. (1982) was approximated by ( |35C| ) only (see Figure 15). 

Setting in accordance with Bisnovatyi-Kogan and Fridman (1969) the lon- 
gitude electron temperature as ~ 2 • 10^^ K, that is T2 = 2 • 10^^ K and 
a = 7 • 10^'* 2X1^ we obtain for the magnetic field strength 5 = 4- 10^° Gs, 

cm ' G G : 

8 • 10^" Gs, and 4-10^° Gs respectively. Here we estimate the spectral form of the 
cyclotron line averaged over the pulsar period, supposing a uniform distribution 
of f{p II) over the polar cap. In this model the beam of the cyclotron feature 
is determined by the number distribution of the emitting relativistic electrons, 
moving predominantly along the magnetic field, over the polar cap. 

5.3.2 Model of the X-ray spectrum of Her X-1. 

In order to obtain the whole experimental spectrum of the Her X-1 the following 
model of the hot spot is considered (see Figure 16). A coUisionless shock wave is 
generated in the accretion flow near the surface on the magnetic pole of a neutron 
star. In it's front the ultrarelativistic electrons are generated having small pitch- 



angle values, so the condition (335) is fuUfiUed automatically. Behind the shock 
there is a hot turbulent zone with a temperature Tg, and optical depth r^, and 
under this zone a heated spot on the surface of the neutron star with a smaller 
temperature is situated. 

The whole X-ray spectrum of pulsar Her X-1 taken from McCray et al. 
(1982) is represented in Figure 15 by the solid line, and in Figure 16 together 
with the schematic picture of the emition region of Her X-1. 
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Table 2: Observations of the cyclotron line in Her X-1 



Date 


References 


Wmax (KeV) 


Line width (KeV) 


May 1978 
September 1977 
February 1978 
April 1980 
May 1980 
September 1980 


Triimper et al. (1978) 
Voges et al. (1982) 
Gruber et al. (1980) 
Gruber et al. (1980) 

Ubertini et al. (1980) 
TuUer et al. (1984) 


58 
51 
48 
54 
49.5 
39 


2ll? 
28 ± 7 

-,1+14 

-•^-"^-11 
181^ 

97+21 
^ ' -20 



There are three main regions in the spectrum: a quasi-Plankian spectrum 
between 0,3 and 0,6 KeV, that is generated near the magnetosphere of the X-ray 
pulsar; power-law spectrum (0.6 20) KeV with a rapid decrease at 20 KeV, 
and the cyclotron feature. 

The power-law spectrum area appears as follows. A surface emits the black- 
body spectrum with a temperature T,- Travelling through the turbulent zone 
this radiation is comptonized. This comptonized spectrum has been calculated 
according to Sunyaev and Titarchuk (1980). Setting the neutron star radius 
equal to 10 km, distance from the X-ray pulsar 6 Kps, hot spot area S = 
2 • 10^^ cm^, we have found the best approximation conditions at = 1 KeV, 
Te = 8 KeV, Te = 14. The best approximation of the X-ray spectrum of the 
pulsar Her X-1 is represented in Figure 15 by the dashed line. It agrees nicely 
with the experimental curve. 

Observation of the variability of the cyclotron line is reported in Mihara 
et al. (1997). Ginga detected the changes of the cyclotron energies from 4 
pulsars. The change is as much as 40 % in the case of 4U 0115-1-63. Larger 
luminosity of the source corresponds to smaller average energy of the cyclotron 
feature. These changes might be easily explained in our model. The velocity of 
the accretion flow decreases with increasing of the pulsar's luminosity because 
localy the luminosity is close to the Eddington limit. As a result the shock 
wave intensity drops as well as the energy of the ultrarelativ istic electrons in 
it's front. The cyclotron energy decreases in accordance with ( |352|) . 

Observations on Beppo-Sax and RXTE satellites (Santangelo et al., 1999; 
Heindl et al., 1999) revealed 4 cyclotron-like spectral features (harmonics) in 
the spectum of the X-ray pulsar 4U 0115+63. Interpretation of these features 
is still not completed. 

6 Magnetohydrodynamic Origin of Jets from 
Accretion Disks 
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6.1 Introduction 



Powerful, highly-coUimated, oppositely directed jets are observed in active ga- 
laxies and quasars (see for example Bridle & Eilek 1984), and in the "micro- 
quasars" - old compact stars in binaries (Mirabel & Rodriguez 1994; Eikcnbcrry 
et al. 1998). Further, highly coUimated emission line jets are seen in young 
stellar objects (Mundt 1985; Biihrke, Mundt, & Ray 1988). Different ideas 
and models have been put forward to explain astrophysical jets (sec reviews 
by Begelman, Blandford, & Rees 1984; Bisnovatyi-Kogan 1993b; Lovelace et al. 
1999b). Recent observational and theoretical work favors models where twisting 
of an ordered magnetic field threading an accretion disk acts to magnetically 
accelerate the jets as first proposed by Lovelace (1976) and Blandford (1976). 
The nature of the ordered magnetic field threading an accretion disk envisioned 
by Bisnovatyi-Kogan and Ruzmaikin (1976) is shown in Figure 17. Figure 18 
shows the outflows from a disk with such an ordered field from Lovelace (1976). 
Two main regimes have been considered in theoretical models, the hydromag- 
netic regime where the energy and angular momentum is carried by both the 
electromagnetic field and the kinetic flux of matter, and the Poynting flux regime 
where the energy and angular momentum outflow from the disk is carried pre- 
dominantly by the electromagnetic field. 

The theory of the origin of hydromagnetic outflows has been developed by 
many authors (Blandford & Payne 1982; Pudritz & Norman 1986; Sakurai 1987; 
Koupelis & Van Horn 1989; Lovelace. Berk & Contopoulos 1991; Pelletier & 
Pudritz 1992; Konigl & Ruden 1993; Lovelace, Romanova, & Contopoulos 1993; 
Cao & Spruit 1994; Contopoulos & Lovelace 1994; Contopoulos 1995; Ostriker 
1997). Physical imderstanding of the hydromagnetic outflows from disks has 
developed from MHD simulations (Uchida & Shibta 1985; Shibata & Uchida 
1986; Bell 1994; Ustyugova et al. 1995, 1999; Koldoba et al. 1995; Romanova 
et al. 1997, 1998; Meier et al. 1997; Ouyed & Pudritz 1997; Krasnopolsky, Li, 
k Blandford 1999; Koide et al. 2000). 

Recent work on hydromagnetic outflows assumes the geometry sketched in 
Figure 19. A strong case for hydromagnetic jets as an explanation of jets in 
protostellar systems emerges because the temperature of the inner regions of 
these systems is insufflcient to permit driving by thermal or radiation pres- 
sure (Konigl & Ruden 1993). Part of the investigations have been analytical 
or semi-analytical and outgrowths of the self-similar solution of Blandford & 
Payne (1982) (Pudritz & Norman 1986; Konigl 1989; Pelletier & Pudritz 1992; 
Contopoulos & Lovelace 1994). The outflows in this model are often referred 
to as "ccntrifugally driven" owing to the driving force close to the disk: Close 
to the disk, MHD fluid particles move along the magnetic flcld lines (sometimes 
described as "beads on a wire"). If the poloidal magnetic field lines threading 
the disk diverge away from the z— axis as shown in Figure 17 (making an an- 
gle with the 2;-axis of more than 30°), then the siun of the gravitational and 
centrifugal forces is outward, away from the disk. The self-similar models are 
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unsatisfactory in the respect that they must be cutoff at smaU cyhndrical radii, 
r < Vmin- This is the most important region of the jet flow. Observations of 
optical stellar jets (Mimdt 1985) reveal jet velocities ^ 200 — 400 km/s, which 
are comparable to the Keplerian disk velocity close to the star's surface. This 
suggests that the jets originate from the inner region of the disk close to the 
star (Shu ct al. 1988: Pringlc 1989). 

The limitations of analytical models has motivated efforts to study jet forma- 
tion using MHD simulations. Simulation studies of hydromagnetic jet formation 
have addressed two main regions: The jet formation region where the matter en- 
ters with sub slow-magnetosonic speed and exits with super fast-magnetosonic 
speed. The second region includes the disk and the problem of the Velikhov 
(1959) - Chandrasekhar (1981) instability ( Balbus-Hawley 1998) and the re- 
sulting 3D MHD turbulence. A number of studies have addressed the coupled 
problem of the disk and near jet regions (Uchida & Shibata 1985; Shibata & 
Uchida 1986; Stone & Norman 1994; Bell and Lucek 1995;). MHD simulations 
of the near jet region have been carried out by several groups (Bell 1994; Ustyu- 
gova et al. 1995, 1999; Koldoba et al. 1995; Romanova et al. 1997, 1998; Meier 
et al. 1997; Ouycd & Pudritz 1997). 

The powerful jets observed from active galaxies and quasars are probably 
not hydromagnetic outflows but rather Poynting flux dominated jets. The mo- 
tion of these jets measured by very long baseline interferometry correspond to 
bulk Lorentz factors of F = 0(10) which is much larger than the Keplerian disk 
velocity predicted for hydromagnetic outflows. Furthermore, the low Faraday 
rotation measures observed for these jets at distances <kpc from the central 
object implies a very low plasma densities. Similar arguments indicate that 
the jets of microquasars are not hydromagnetic outflows but rather Poynting 
jets. Poynting jets have been proposed to be the driving mechanism for gamma 
ray burst sources (Katz 1997) Theoretical studies have developed models for 
Poynting jets from accretion disks (Lovelace, Wang, & Sulkanen 1987; Colgate 
& Li 1998; Lyndcn-Bell 1996; Romanova & Lovelace 1997; Lcvinson 1998; Ro- 
manova 1999; Lovelace et al. 2001; and Li et al. 2001). Stationary Poynting 
flux dominated outflows were found by Romanova et al. (1998) and Ustyugova 
et al. (2000) in axisymmctric MHD simulations of the opening of magnetic loops 
threading a Keplerian disk. 

Here, we first review recent results on MHD simulations of hydromagnetic 
jet formation and later discuss the Poynting flux regime. Section 4.3 discusses 
general considerations of hydromagnetic outflows. Sections 4.4 discusses MHD 
simulations which give non-stationary and stationary hydromagnetic outflows. 
Sections 4.5 describes the Poynting flux regime which remains to be fully ex- 
plored by simulations. Section 4.6 gives the conclusions. 
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6.2 Basics of MHD Outflows from Disks 



The main forces which drive a hydromagnetic outflow from a disk threaded by 
a magnetic field are the centrifugal force and the magnetic pressure gradient 
force. If disk has a hot corona, the pressure gradient may also be important. 
We neglect the radiative force but in this regard see Phinney (1987). Accreting 
matter of the disk carries magnetic field inward thus generating a Br component 
of the magnetic field as sketched in Figure 19. On the other hand, rotation of 
the disk acts to generate a toroidal component of the field (< if i?^ > 0). 

For a sufficiently inclined magnetic field {6 in Figure 19 sufficiently large), 
outflows can result from the centrifugal force (Blandford & Payne 1982) and/or 
the magnetic pressure gradient force (— Vz B"^/ {Sir)) (Lovelace et al. 1989, 1991; 
Koupelis & Van Horn 1989). This depends on the ratio of energy densities at 
the base of the outflow at the inner radius of the disk denoted r^. Thus, the 
main parameters are eth = {cs/vk)^ and Sb = {va/vk)1, where vk is Keplerian 
velocity, Cg the sound speed, va the Alfven speed, and the i subscript indicates 
evaluation on the surface of the disk at its inner radius, r = ri. For Sb ^ ^ 
the outflow is magneto-centrifugally driven, whereas for Sfh ~ 1, the flow is 
thermally driven. 

Processes in the disk arc of course coupled to the outflows (Lovelace et 
al. 1994, 1997; Falcke, Malkan, & Biermann 1995). However, it is difficult to 
simultaneously simulate the disk and the outflow regions because the time scales 
of the accretion and outflow are in general very different. The radial accretion 
speed of the disk is typically \vr\ <C Cs <C vk, where Cs is the sound speed in 
the disk, and vk is the Keplerian velocity. On the other hand the velocity of 
magnetic outflows from disks is of the order of vk- Clearly, it is "easier" (less 
computation time is required) to simulate the outflows. 

Furthermore, processes in the disk may involve the small scale MHD insta- 
bility of Chandrasekhar-Vclikhov instablity emphasized by Balbus, and Hawley, 
and therefore require high spatial resolution. Stone & Norman (1994) attempted 
to simultaneously simulate the internal MHD dynamics of a disk and MHD dy- 
namics of outflows. This proved impractical because essentially all of the spatial 
resolution was needed for treating the unstable dynamics of the disk. Also, there 
was the problem that the initial configuration was far from equilibium. Simula- 
tion of the internal MHD disk dynamics has led several groups to the problem 
of simulating 3D MHD turbulence in a sheared flow of a patch of a disk (for 
example, Hawley et al. 1995; Brandenburg ct al. 1995). This is a much larger 
project than that of understanding MHD outflows. At the same time it is widely 
thought, and observations of cataclysmic variables support the view, that the 
disk turbulence - including MHD turbulcnc;c - can be modeled approximately 
using the Shakura (1973), Shakura & Sunyacv (1973) "alpha" viscosity model 
(Eardley & Lightman 1975; Coroniti 1981; see also section 1.5). In contrast 
with the internal disk dynamics, there is theoretical and simulation evidence 
that the outflows can be treated using axisymmetric (2D) MHD (Blandford & 
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Payne 1982; Lovelace et al. 1991; Ustyugova et al. 1995). Here, we consider 
outflows from a disk represented as a boundary condition. This approach has 
subsequently been adopted by other groups (Meier et al. 1997; Ouyed & Pudritz 
1997). This treatment of the disk is justified for outflows from a disk where the 
accretion speed is small compared with the Keplerian speed (Ustyugova et al. 
1995). 



6.3 Theory of Stationary MHD Flows 

The theory of stationary, axisymmetric, ideal MHD flows was developed by 
Chandrasekhar (1956), Woltjer (1959), Mestcl (1968), KuUkovskyi & Lyubimov 
(1962), and others. Under these conditions the MHD equations can be reduced 
to a single equation for the "flux function" \l/(r, z) in cylindrical (r, (p, z) coor- 
dinates (Heinemann & Olbert 1978; Lovelace et al. 1986). The flux fimction ^' 
labels flux surfaces so that ^'(r, 2;) =const represents the poloidal projection of 
a field line. The equation for ^ is commonly referred to as the Grad-Shafranov 
equation (Lovelace et al. 1986). The discussion of this section follows that of 
Ustyugova et al. (1999). 



6.3.1 Integrals of Motion 

For axisymmetric conditions the flow field can be written as v = Vp-|-t;^e,^ where 
Vp is the poloidal (r, z) component, = Lur > is the toroidal component, and 
e<^ is the unit toroidal vector. Similarly, the magnetic field can be written as 
B = Bp + B^e^. The ideal MHD equations then imply that certain quantities 
are constants on any given flux surface 5'(r, z) =const or equivalently they are 
constants along any given stream line or a given magnetic field line. These 
integrals are functions of ^ (see for example Lovelace et al. 1986), 

= ^^^^^ 

_ ^ (354) 
K 

uj - = n(^) (355) 
Anpr 

S = S'(*) (356) 
0^7.2 1 

«; + $-^ + ^ + -(c^-0)V =£;(*) (357) 

Here, S is the entropy, w is the enthalpy, and $ is the gravitational potential. 
The quantity K corresponds to the conservation of mass along a streamline, A 
to the conservation of angular momentum, Q to the conservation of helicity, S 
to the conservation of entropy, and E (Bernoulli's constant) to the conservation 
of energy. 
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The remaining MHD equation (which cannot be written in the integral form) 
is the Euler force equation across the poloidal magnetic field line (Bogovalov 
1997), 



de 



cos 61 2 2 \ 1 ^ 

r ^ ^ p on 



which is equivalent to the Grad-Shafranov equation. Here, d/dn is the derivative 
in the direction perpendicular to magnetic field lines and directed outward from 
the axis, 9 is the angle of inclination of the poloidal magnetic field line away 
from the 2— axis, s is the distance from the disk along a magnetic field line, 
VAp = |Bp|/^/47rp and VA4, = l^^l/V^^rp are the poloidal and azimuthal Alfven 
velocities. The quantity 89/ ds is the curvature of magnetic field line. The 
first two terms in equation (6) are determined by the non-diagonal (tension) 
part of the stress tensor, pviVk + {p + B^/Stt) Sik — BiBk/^Tr. The third term 
is determined by the total (matter plus magnetic) pressure p + B^/Stt and the 
gravity force d^/dn. 

To clarify the physical sense of the integrals of motion, it is useful to derive 
the fiuxes of mass, angular momentum (about the z— axis), and energy. The 
corresponding conservation laws for stationary conditions are 
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Because Vp || Bp, the vector flux densities are directed along the field lines. 

Consider the fluxes through an annular region with surface area element dS. 
The matter flux through the axisymmetric surface S extending out from the 
z— axis is 

Tm = y" dS • pvp = i- JdS- BpKi^) , (362) 



where we took into account the integral (353). dS ■ Bp is the magnetic flux 



through the annular region bounded by flux surfaces ^ and 5* + d^*. Thus 
we can change from space integration to integration over Because Br = 
— {l/r)d'^/dz and B^ — (l/r)9^/9r, we have 



•^m(*) = \j d^'K{^') , 



(363) 
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where 5* = corresponds to the axis. Similarly, 
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(364) 
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Thus, Kd'i>/2 is the matter flux between the flux surfaces separated by d'^, 
AKd'^/2 is the angular momentum flux, and {E + Afl) K d'^/ / 2 is the energy flux. 
Note that A(^') is the specific angular momentum carried along the magnetic 
field line ^' — const, E{'^) + Ail{^^) is the specific energy, and ^{'^) is the 
angular velocity of the disk at the point where the magnetic field line or flux 
surface 5' = const intersects the disk (for |vp| ^ at the disk). 



6.3.2 Alfven Surface 

Conditions at the Alfven surface are known to be important for the global prop- 
erties of MHD flows (Weber & Davis 1967). Equations (2) and (3) constitute a 
linear system of equations for lo and . The determinant of this system is zero 
if = All p. Under this condition a solution exists if A = r^il (Weber & Davis 
1967) which corresponds to Vp = Bp/ ^JAirp = v^p. This is the condition which 
defines the Alfven surface. Figure 20 shows the Alfven surface for a sample case. 
The radius at which this field line intersects the disk is rd{'^)- The radius at 
which it crosses the Alfven surface is ryi(^'). The density at this point on the 
Alfven surface is pa{^)- Thus, 



Pa(^)^ K^{'1!)/Att , r^(l') = A(1')/17(^') . 



Equations (354) and (355) give 



1 - PAr\lpr'^ 
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Taking into account equations (366) - (36S), one can express the fluxes of mass, 
angular momentum and energy, using only the values of physical quantities on 
the Alfven surface: 









(369) 



(370) 



(371) 



6.3.3 Forces 

For understanding the plasma acceleration, we project the different forces onto 
the poloidal magnetic field lines. As mentioned, in a stationary state, matter 
flows along the poloidal magnetic field lines. The acceleration in the poloidal 
(r, z) plane is 

(vp • V)vp + U0(e0 • V)(w0e0) . (372) 

The last term represents the centrifugal acceleration —{v'^/r)er = —nJ^er. To 
get the force per unit mass along a magnetic field line, we multiply the Euler 
equation by a unit vector b parallel to Bp. This gives 

/ = c.Vsin0-i^-^ + -^b.[(VxB)xB] , (373) 
p OS OS 47rp 

where 9 is the inclination angle of the field line to the z— axis. The final term 
of ( |373D is the projection of the magnetic force in the direction of b, which can 
be transformed to 

/M = ^b.[(VxB)xB]^-, ' ^('^^^^^ 



Anp Snpr^ ds ' 

which is useful for understanding our results. 

When magnetic field lines are inclined outward, away from the symme- 
try axis, the gravitational force fc = d^/ds opposes the outfiow of matter 
from the disk. If the matter is relatively cold then the pressure gradient force 
fp = — (1/yo) (dp/ds) is unimportant. Then matter is accelerated outward by the 
centrifugal force fc = uj^rsin0 and/or the magnetic force /m- This determines 
the driving mechanisms of the outflow, centrifugal and/or magnetic. The cen- 
trifugal force always acts to accelerates matter outward if the distance between 
magnetic field line and the axis increases. However, under the same condition. 
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the magnetic force can act to accelerate or decelerate the flow depending on the 
direction of —W{rB^)'^ (see discussion in Ustyugova et al. 1999). 

Thus, magnetic and centrifugal forces may both accelerate matter, but this 
depends on the configuration of magnetic field and current-density lines. 



6.3.4 Collimation 



Consider now the collimation of the flow. From equation (35S), taking into 
account that cos6 = dr/dn^ we have 



At large distances from the Alfven surface r ^ r^, the density p ^ pA, but 
values pr^ and remain finite (Heyvaerts and Norman 1989). Then v\p = 



{p/pAf'v'i <C Vp, so that the second term on the left-hand side of (|374|) is 



negligible. On the right-hand side, only the first term is important for r ^ r^. 



Then, equation (374) simplifies to 
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[rB^Y ■ (375) 



Thus, the curvature of magnetic field lines in the region r ^ r^i is determined 
by the gradient (ri?^)^. This force can act either to "collimate" the fiow or 
"anticoUimate" it depending on the orientation of the rBff,{r,z) =const and 
\E'(r, z) =const surfaces (Ustyugova et al. 1999). The simulation study of Ustyu- 
gova et al. (1999) in fact finds only very small collimation of the hydromagnetic 
outflows over distances of the order of 100 times the inner disk radius. Recent 
work by Ustyugova et al. (2001) extents this conclusion by showing only small 
collimation out to distances of > 1000 times the inner disk radius. The colli- 
mation observed in protostellar jets at small distances, for example, in HH 30 
observed with the space telescope, may result from the pressure of the external 
medium (Lovelace et al. 1991). 

The conclusion by Heyvaerts and Norman (1989) that the magnetic collima- 
tion of MHD jets is a mathematical rather than physical result. The conclusion 
is based on the premise that the entire half-space above the disk is described by 
a stationary solution ^'(r, z) of the Grad-Shafranov equation. In reality, MHD 
jets result from initial conditions at the disk (as in the simulation studies) which 
result in the outward propagation of the jet into a highly conducting ionized 
external (interstellar or intergalactic) medium. At a given time, the region in- 
side the surface r — TZ{z,t) is filled by the jet flow and the region outside it is 
filled by external plasma. The net axial current carried by the jet 2tt rdrJz 
(for say 2; > 0) is necessarily zero for any {z,t). Therefore, the toroidal mag- 
netic field is zero on the surface 7?., and there is no magnetic confinement of the 
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flow. Of course, it is possible that a fraction of the flux \1/ is collimated along 
the 2;— axis while the remaining flux is "anti-collimated" (that is, pushed away) 
from the z— axis. Exactly such a case is discussed in §1.5 below for the case of 
Poynting jets. This case also has the physical condition that the jet carries no 
net magnetic flux, 27r jj^ rdrBz = 0. As in Landau's (1946) treatment of the 
damping of plasma waves, the physically correct treatment of the coUimation of 
MHD jets requires consideration of an initial value problem. 

6.4 Numerical Simulations of MHD Outflows 

In order to test and more fully understand analytical models of stationary out- 
flows, MHD simulations of flows from a disk treated as a boundary condition 
have been carried out by a number of groups. 

6.4.1 Non-Stationary Outflows 

The initial magnetic field configuration was chosen so that the magnetic field 
was significantly inclined to the disk {9 > 30°) over most of the disk surface 
(Ustyugova et al. 1995; Koldoba et al. 1995). The simulations involve solving 
the complete system of ideal non-rclativistic MHD cqiiations using a Godunov- 
type code assuming axisymmetry but taking into account all three components 
of velocity and magnetic field. Matter of the corona was initially in thermal 
equilibrium with the gravitating center. At t = 0, the disk is set into rotation 
with Keplerian velocity and at the same time matter is pushed from the disk 
with a small poloidal velocity equal to a fraction of the slow magnetosonic 
velocity ( Vp = avsm, with a = 0.1 — 0.9). A relatively high temperature and 
small magnetic field was considered. We found that at the maximum of the 
outflow, matter is accelerated to speeds in excess of the escape speed and in 
excess of the fast magnetosonic speed within the simulation region (^ lOOr,;). 
The acceleration is due to both thermal and magnetic pressure gradients. The 
outflow collimates within the simulation region due to strong amplification, 
'wrapping up' of the toroidal magnetic field and the associated pinching force. 

However, the outflows are not stationary. The matter flux grows to a peak 
and then decreases to relatively small values. The strong collimation of the 
outflow reduces the divergence of the field away from the z— axis {6 < 30°) and 
this "turns off" the outflow of matter and leads to flow velocities less than the 
escape speed. Thus, this simulation is an example of a temporary outburst of 
matter to a jet. Unfortunately, this type of flow has a significant dependence 
on the initial conditions. 

6.4.2 Stationary Outflows 

More recently, stationary magnetohydrodynamic outfiows from a rotating accre- 
tion disk have been obtained by time-dependent axisymmetric simulations by 
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Romanova et al. (1997) and systematically analyzed by Ustyugova et al. (1999). 
The initial magnetic field in the latter work was taken to be a split-monopole 
poloidal field configuration (Sakurai 1987) frozen into the disk. The disk was 
treated as a perfectly conducting, time-independent density boundary [p{r)] in 
Keplerian rotation which is different from the earlier specification of a small 
velocity outflow (§4.3.1, Ustyugova et al. 1995). The outflow velocity from the 
disk is determined self-consistently from the MHD equations. The temperature 
of the matter outflowing from the disk is small in the region where the magnetic 
field is inclined away from the symmetry axis (c^ ^ but relatively high 
{cj. ~ vj^) at very small radii in the disk where the magnetic field is not inclined 
away from the axis. We have found a large class of stationary MHD winds. 

Figure 20 shows the nature of the stationary MHD outflows found by Ustyu- 
gova et al. (1999). The outflows are approximately spherical, with only small 
coUimation within the simulation region. The collimation distance over which 
the flow becomes coUimated (with divergence less than, say, 10°) is much larger 
than the size of our simulation region. Because the flow is super fast-magneto- 
sonic (see below), the flow is not expected to magnetically coUimate at larger 
distances. Notice for example that the solar wind is also super fast-magnetosonic 
and essentially uncoUimated. Collimation may however result from the pressure 
of external plasma (Lovelace et al. 1991). 

Figure 21 shows the variation of the flow speed and other characteristic 
speeds along a sample streamline. The outflow accelerates from thermal velocity 
(~ Cs) to a much larger asymptotic poloidal flow velocity of the order of one-half 
\fGMjri where M is the mass of the central object and rj is the inner radius of 
the disk. This asymptotic velocity is much larger than the local escape speed and 
is larger than fast magnetosonic speed by a factor of '-^ 1.75. The acceleration 
distance for the outflow, over which the flow accelerates from ~ to, say, 90% 
of the asymptotic speed, occurs at a flow distance of about SOr^. Close to the 
disk the outflow is driven by the centrifugal force while at all larger distances 
the flow is driven by the magnetic force which is proportional to — V(ri?0)^, 
where B^f, is the toroidal field. 

The stationary MHD flow solutions allow us (353) to compare the results 
with MHD theory of stationary flows, ( [354[ ) to investigate the influence of differ- 
ent outer boundary conditions on the flows, and (355) to investigate the influence 
of the shape of the simulation region on the flows. The ideal MHD integrals 
of motion (constants on flux surfaces discussed by Lovelace et al. 1986) were 
calculated along magnetic fleld lines and were shown to be constants with accu- 
racy 5 — 15%. Other characteristics of the numerical solutions were compared 
with the theory, including conditions at the Alfven surface. 

Different outer boundary conditions on the toroidal component of the mag- 
netic field can significantly influence the calculated flows. The commonly used 
"free" boundary condition on the toroidal field leads to artificial magnetic forces 
on the outer boundaries, which can give spurious collimation of the flows. New 
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outer boundary conditions which do not give artificial forces have been proposed 
and investigated by Ustyugova et al. (1999). 

The simulated flows may also depend on the shape of the simulation region. 
Namely, if the simulation region is elongated in the z— direction, then Mach 
cones on the outer cylindrical boundaries may be partially directed inside the 
simulation region. Because of this, the boundary can have an artificial influ- 
ence on the calculated flow. This effect is reduced if the computational region 
is approximately square or if it is spherical as in Figure 20. Simulations of 
MHD outflows with an elongated computational region can lead to artificial 
coUimation of the flow. 

Recent simulation studies have treated MHD outflows from disks with more 
general initial B field configurations, for example, that where the poloidal field 
has different polarities as a function of radius (Hayashi, Shibata, & Matsumoto 
1996; Goodson, Winglee, & Bohm 1997; Romanovaet al. 1998). The differential 
rotation of the foot-points of B field loops at different radii on the disk surface 
causes twisting of the coronal magnetic field, an increase in the coronal magnetic 
energy, and an opening of the loops in the region where the magnetic pressure 
is larger than the matter pressure (/3~1) (Romanova et al. 1998). In the region 

where /3~1, the loops may be only partially opened. Current layers form in the 
narrow regions which separate oppositely directed magnetic field. Reconnection 
occurs in these layers as a result of the small numerical magnetic diffusivity. 
In contrast with the case of the solar corona, there can be a steady outflow of 
energy and matter from the disk surface. We find that the power output mainly 
in the form of a Poynting flux. Opening of magnetic field loops and subsequent 
closing can give reconnection events which may be responsible for X-ray flares 
in disks around both stellar mass objects and massive black holes (Hayashi et 
al. 1996; Goodson, Winglee, & Bohm 1997; Romanova et al. 1998). 

6.5 Poynting Jets 

In a very different regime from the hydromagnetic flows discussed in §2 - 4, 
a Poynting outflow or jet transports energy and angular momentum mainly 
by the electromagnetic fields with only a small contribution due to the matter 
(Lovelace et al. 1987; Colgate & Li 1998; Romanova et al. 1998; Ustyugova et 
al. 2000; Li et al. 2001; Lovelace et al. 2001). A steady Poynting jet can be 
characterized in the lab frame by its asymptotic (large distance) magnetic field 
= —Bo[ro/rj{z)] at the jet's edge, r = rj{z), where ro is the jet's radius at 
z = and Bq is the poloidal field strength at this location. The electric field in 
the jet is E = — v x B/c and consequently the energy flux or luminosity of the 
jet is Lj = vB^r^/8 ~ 4x 10^^eTg/s{v/c){ro/10^\mf{Bo/10^Gf. Propagating 
disturbances in such field dominated jets provide a simple, but self-consistent 
physical model for the gamma ray flares observed in Blazars (Romanova & 
Lovelace 1997; Levinson 1998; Romanova 1999). Owing to pair production 
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close to the black hole, the main constituent of a Poynting flux jet may be 
electron-positron pairs. 



6.5.1 Theory of Poynting Jets 

Consider the coronal magnetic field of a differentially rotating Keplerian ac- 
cretion disk for a given poloidal field threading the disk. That is, the disk 
is perfectly conducting with a very small accretion speed. Further, consider 
"coronal or force-free" magnetic fields in the non-relativistic limit. We use 
cylindrical (r, 4>^ z) coordinates as in §3 and consider axisymmetric field con- 
figurations. Thus the magnetic field has the form B = Bp + , with 
Bp = BrV + Bzi ■ Because V-B = 0, B=:VxA with A the vector po- 
tential. Consequently, Br — —{l/r){d'^/dz) ,Bz — {\ / r){d'^ / dr) , where 
\E'(r, z) = rA^{r,z). The \l/(r, z) = const lines label the poloidal field lines; 
that is, (Bp-V)* = = B • V^* = 0. Note that 27r*(r, z) is the magnetic flux 
through a horizontal, coaxial circular disk of radius r. 

The magnetic field threading the disk at z = is assumed to evolve slowly so 
that it can be considered approximately time-independent, ^'(r, z = 0) = ^'o(r). 
However, the magnetic field above the disk will in general be time-dependent, 
\E' — \E'(r, z,t), due to the differential rotation of the disk. Figure 22 shows a 
possible vacuum poloidal magnetic field threading the disk. This field may be 
regarded as the "initial" field for the considered case where the space exterior 
to the disk is filled with a low density plasma. 

The non-relativistic equation of plasma motion in the corona of an accretion 
disk is 



pdv/dt -Vp + pg + J X B/c, (376) 

where v is the fiow velocity, p is the pressure, and g is the gravitational ac- 
celeration. The equation for the B field is V x B — ^J, because the dis- 
placement current is negligible in the non-relativistic limit. In the "coronal or 
force-free plasma limit," the magnetic energy density B^/(87r) is much larger 
than the kinetic or thermal energy densities; that is, sub-Alfvenic fiow speeds v. 



<C f ^ = B^/47rp, where va is the Alfven velocity. Equation (376) then sim- 
plifies to « J X B; therefore, J = AB (Gold & Hoyle 1960). Because V • J = 0, 
(B • V)A = and consequently A — A(^'), as well-known. Thus Ampere's law 
becomes V x B = 47rA(*)B/c. 

The r and z components of this equation imply rB^ — H{'^), and dH{'^)/<Ni 
= 47rA(^')/c, where H{'^) is another function of . Thus, H{'^) = const are lines 
of constant poloidal current density; Jp — (c/47r)(di?/d^')Bp so that [3p-V)H = 
0. The toroidal component of Ampere's law gives 

A** = -i/(*)_L^ , (377) 
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which is the Grad-Shafranov equation for \1/ (see for example Batemann 1980; 
Lovelace et al. 1986). Here, A* = d'^/dr'^ - {l/r){d/dr) + d'^ /dz^ is the adjoint 
Laplacian operator. Note that A** = r(V^ - 1/r^)^^ and that HidH/d"^) = 
AwrJ^/c. From Ampere's law, § d\-B = {4tt/c) J dS-3, so that rB^{r, z) = H{^/) 
is (2/c) times the current flowing through a circular area of radius r (with 
normal z) labeled by ^{r,z)= const. Equivalently, — -ff[^'(r,0)] is (2/c)x the 
current flowing into the area of the disk < r. For all cases studied here, — i?(\l/) 
has a maximum so that the total current flowing into the disk for r < is 
Itot = \{-H)max, where r„j is such that -H[*(r„,0)] = {-H)max so that 
is less than the radius of the O— point, tq. The same total current Itot flows out 
of the region of the disk r = to ro- 

The fimction H{'^) must be determined before the GS equation can be 
solved. Lynden-Bell and Boily (1994) determined the form of ff(\E') by requiring 
that the solution ^'(r) be a self-similar function of |r|. Their stationary solution 
has no energy or angular momentum outflow from the disk. However, in general 
the magnetic field threading a disk will have definite length-scales and the field 
solution will not be self-similar. 

For definiteness we consider the initial value problem where the disk at t = 
is threaded by a poloidal field such as shown in Figure 22. Such a field can arise 
from a disk-dynamo driven by star-disk collisions as discussed by Pariev, Finn, 
& Colgate (2001). We assume for simplicity that the field threading the disk 
has the symmetry of a dipole although it is predicted in the reference above 
that the quadrupole field grows somewhat faster. 

The form of H{'^) is then determined by the differential rotation of the disk: 
The azimuthal twist of a given field line going from an inner footpoint at ri to 
an outer footpoint at r2 is fixed by the differential rotation of the disk. The field 
line slippage speed through the disk due to the disk's finite magnetic diffusivity 
is estimated to be negligible compared with the Keplerian velocity. For a given 
field line we have rd(})/B^ = dSp/Bp, where dSp = V dr'^ + dz^ is the poloidal 
arc length along the field line, and Bp = ^/B^ + B^. The total twist of a field 
line loop is 

A,m = j\s,:^ = -Hml'^^. (378) 

with the sign included to give A(j) > 0. For a Keplerian disk around an object 
of mass M the angular rotation is lok = ^jGMjr'^ so that the field line twist 
after a time t is 

A<^(*)=a;ot f-1 = (wo i) ^(*/*o) (379) 



where ro is the radius of the O— point, wq = \/ GM/r^, and is a dimensionless 
function (the quantity in the square brackets). At sufficiently small r\ one 
reaches the inner radius of the disk r, (<C ro) where we assume wk saturates 
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at the value coxi — ^/GM/rf. For the dipole-like field of Figure 22, w 
39/8(vE,g/4,)3/4 foj. ^/^^ ^ ^Yiile T w 3.64(1 - */*o)^/^ for 1 - ^/fo < 1- 

The quantity T = u>ot (in radians) we refer to as the "twist" of the disk. 
For the dashed field lines of Figure 22, T is the twist A(/)(^'i) of the field line 
with ^i/^Q « 0.932 which has footpoints at ri w 0.73ro and r2 ~ 1.4ro. For 
\l//'i'o = 0.1, where the footpoints are at ri « 0.14ro and r2 ~ T.Srp, the field 
line twist is Acj) « 18.3 rad. for T 1. 

For example, for an accretion disk around a massive M — M^IO^Mq black 
hole in the nucleus of galaxy, the twist is T = {t/3.17d){ro/l0^^cmf/^/y/M^. 
The inner radius of the disk is « Mg9 x lO'^'^ cm for a Schwarzschild black 
hole. 

The Grad-Shafranov equation ( j377| ) is an elliptic equation and can be readily 
solved by the method of Successive Over-Relaxation (SOR) (see for example. 
Potter 1973). Figure 23 shows a sample high twist solution exhibiting a Poynting 
jet (Lovelace et al. 2001). These high-twist solutions consist of a region near 
the axis which is magnetically collimated by the toroidal field and a region 
far from the axis, on the outer radial boundary, which is anti- collimated in the 
sense that it is pushed against the outer boundary. The field lines returning to 
the disk at r > tq are anti-collimated by the pressure of the toroidal magnetic 
field. 

Figure 24 shows a three dimensional view of two sample field lines for the 
same case as Figure 23. 



6.5.2 Analytic Model of Poynting Jets 

Figure 24 shows that most of the twist A(j> of a field line occurs along the jet 
from z = to Z„iax- Because —r^d(t>/H{'^) = dz/Bz, we have 

A0(^) _ MTj-^/-^^) ^ Zrna. 

where r^Bz{r, z) is evaluated on the straight part of the jet at r = r(^). In the 
core of the jet \l/ <C ^"0, we have « 3^/^(5'o/5')'^^^, and in this region we can 
take 

^.C*o(^)\ and H^-K[^){£j , (381) 

where C, (Z,/C, and s are dimensionless constants. The GS equation for the 
straight part of the jet, '^rr-'^r/r = -HH' implies q = l/(l-s) and C^(i-'*) = 
s(l - s^K? /{I - 2s). In turn, equation ( |38l| ) implies s = 1/4 so that q = 4/3, 
C = [9/32]2/3/C4/3, and /C = 'i^/H{rou:ot/ Z^^,). 

In order to have a Poynting jet, we must have K, larger than « 0.5. For the 
case of uniform expansion of the top boundary, Z^^x = ^ztj this condition is 
the same as Vz < 9.2(roWo)- For the case of Figures 23 and 24, JC ~ 0.844. 
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The field components in tiie straigiit part of tlie jet are 

= -^B. = -^(A) (I) (I2L)^'\ (382) 

Tliese field components of course satisfy the radial force balance equation 
d{Bl)/dr+{\/r'^)d{r'^B'^)/dr = which is equivalent to the GS equation. These 
dependences agree approximately with those found independently in the nu- 
merical simulations of Poynting jets by Ustyugova et al. (2000). On the disk, 
^ Rj 3^/^*0 (?'/?'o)^ for r < Using this and equation (16) for ^'(r) gives 

the relation between the radius of a field line in the disk, denoted r^, and its 
radius in the jet, r/ro = 2^/^3^/*(rd/ro)^/^/C~^. Thus the power laws apply for 
n < r < r2, where n = 'i^/^2^''^rQ{ri/rQf/'^/K. and ra = (25/2/3)ro//C, with r, 
the inner radius of the disk. The outer edge of the Poynting jet has a transition 
layer where the axial field changes from Bz{r2) to zero while (minus) the toroidal 
field increases from —B^{r2) to {—H)max/r2- Using equations (19), which are 
only approximate at r2, gives {—H)max ~ (3/2)^/^/C'S'o/ro. This expression 
agrees approximately with our numerical results. 



6.5.3 MHD Simulations of Poynting Jets 

Ustyugova et al. (2000) did full, axisymmetric MHD simulations of a Keplerian 

disk initially threaded by a dipole-like magnetic field as shown in Figure 22. For 
these simulations the outer boundaries at r = Rmax and z = Z^ax were treated 
as free boundaries following the methods of Ustyugova et al. (1999). These 
simulations established that a quasi-stationary coUimated Poynting jet arises 
from the inner part of the disk while a steady uncoUimated hydromagnetic 
outflow arises in the outer part of the disk. Figure 25 shows the evolution 
towards a stationary Poynting jet. 

Quasi-stationary Poynting jets from the two sides of the disk within tq give 
an energy outflow per unit radius of the disk dEs/dr = rvK{—B^Bz)h, where 
the h subscript indicates evaluation at the top surface of the disk. This out- 
flow is ~ rodEe/dr ~ VK{ro){^-maxlrof ~ lQ''''{eTg/s)rlil^s[B,{{))/(SkGf 
where rois is in units of lO^^cm, and Mg in units of lO^M©. This formula 
agrees approximately with the values derived from the simulations. The jets also 
give an outflow of angular momentum from the disk which causes disk accre- 
tion (without viscosity) at the rate Msir) = —27r'Evr = —2{r'^ /vK){B^Bz)h ^ 
2^mai/[^o''^<(''o)], where S is the disk's surface mass density. (Lovelace et al. 
1997). The Poynting jet has a net axial momentum flux = (1/4) / rdr{B^ — 
B1) ~ Q.^i^max/roY , which acts to drive the jet outward through an external 
medium. Further, the Poynting jet generates toroidal magnetic flux at the rate 

max- 

For long time-scales, the Poynting jet is of course time-dependent due to 
the angular momentum it extracts from the inner disk (r < Tq). This loss of 
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angular momentum leads to a "global magnetic instability" and collapse of the 
inner disk (Lovelace et al. 1997). An approximate model of this collapse can 
be made if the inner disk mass Md is concentrated near the O— point radius 
ro{t), if the field line slippage through the disk is negligible (Lovelace et al. 
1997), ^max =const, and if {-rB^)max ~ ^max/ro{t) (as found here). Then, 
Mddro/dt = -2*2j„^(G'Mro)"^/^. If U denotes the time at which ro(ti) = ri 
(the inner radius of the disk), then ro{t) = ri[l — {t — ti)/tcoiiY^^ , ior t < ti, 
where tcoU = VGM Mdrf^^/(3^^„2.) is the time-scale for the collapse of the 
inner disk. (Note that the timc-scalc for tq to decrease by a factor of 2 is 
~ Uifo/riY^"^ » ti for ro » r^.) The power output to the Poynting jets is 
E{t) = (2/3)(A£tot/tco«)[l - (i - ^i)/^coH]-^/^ where AE^t = GMMij^n 
is the total energy of the outburst. Roughly, tcoU ^ 2dayM|(Md/Af0)(6 x 
10^^Gcm^/\E'maa;)^ for a Schwarzschild black hole, where validity of the analysis 
requires tcoU > ti- 

6.6 Summciry Regarding Jets from Accretion Disks 

MHD simulations carried out by a number of groups over the last several years 
support the idea that an ordered magnetic field of an accretion disk can give 
powerful outflows or jets of matter, energy, and angular momentum. Most of the 
studies have been in the hydromagnetic regime and find asymptotic flow speeds 
of the order of the maximum Keplerian velocity of the disk, vk%- These flows 
are clearly relevant to the jets from protostellar systems which flow speeds of 
the order of vxi- In contrast, observed VLBI jets in quasars and active galaxies 
point to bulk Lorentz factors F ~ 10 - much larger than the disk Lorentz factor. 
The large Lorentz factors as well as the small Faraday rotation measures point 
to the fact that these jets are in the Poynting flux regime. These jets may involve 
energy extraction from a rotating black hole (Blandford and Znajek 1977; Livio 
et al. 1999). 

7 Discussion 

Observational evidence for massive black holes at the center of our Galaxy and 
in the nuclei of other galaxies (Cherepashchuk, 1996; Ho, 1998; Haswell, 1999; 
Richstone et al. 1998; Gebhardt et al. 2000; Merritt & Ferrarese 2001) make 
it necessary to revise theoretical models of the disk accretion. A large part of 
the high energy radiation originates close to the black hole, where the standard 
accretion disk model is not appropriate. Improvements of the accretion models 
have been made by taking into account the advective terms and including the 
important consequences of ordered and disordered magnetic fields. The conclu- 
sions obtained from ADAF solutions - for optically thin accretion flows at low 
mass accretion rates - are open to question because of omission of all magnetic 
fleld effects. Heating due to the annihilation of the magnetic fleld in the accre- 
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tion flow is predicted to be important. Further, the energy exchange between 
electrons and ions may be increased substantially by the turbulent electric fields 
in the plasma. Account of the magnetic field leads to the conclusion that the 
radiative efiiciency of the accretion is not smaller than about 0.25 the efficiency 
for standard disk accretion (Bisnovatyi-Kogan & Lovelace, 1997, 2000) It is 
expected that more accurate treatment of the relaxation connected with the 
plasma turbulence will increase the efficiency further, making it close to unity 
(see also Fabian and Rees, 1995). 

Some observational data which were interpreted as an evidence for the ex- 
istence of the ADAF regime have disappered after additional accumulation of 
data. The most interesting example of this sort is connected with the claim 
of the proof of the existence of event horizon of the black holes due to mani- 
festation of the ADAF regime of accretion (Narayan et al., 1997). Analysis of 
a more complete set of observational data (Chen et al, 1997) had shown the 
disappearance of the statistical effect claimed as an evidence for ADAF. This 
example shows how dangerous is to base a proof of the theoretical model on the 
preliminary observational data. It is even more dangerous, when the model is 
physically not fully consistent. Then even a reliable set of the observational data 
cannot serve as a proof of the model. The classical example from astrophysics of 
this kind gives the theory of the origin of the elements presented in the famous 
book of Gamov (1952), where the model of the hot imiverse was developed. In 
addition to rich advantages of this model, the author also wanted to explain the 
origin of heavy elements in the primordial explosion, neglecting the problems 
connected with an absence of the stable elements with the number of barions 
equal to 5 and 8. Gamov considered the good coincidence of his calculations, 
where the mentioned problem was neglected, and the observational curve, as a 
proof of his theory of the origin of the elements. Further developments have 
shown that his outstanding theory explains many things, but not the origin of 
the heavy elements which are produced due to stellar evolution. 

It appears difficult explain under-luminous AGNs with ADAF solutions 
where the radiative efficiency is much smaller than for standard accretion disk 
models. Possible explanations include the fact that the mass accretion rates to 
the black holes may have been overestimated. More importantly, other mecha- 
nisms of energy loss may be involved in the accretion flows, for example, in the 
form of accelerated particles, as in the radio-pulsars, where the energy carried 
by the relativistic wind far exceeds the electromagnetic radiation (e.g., the Crab 
pulsar). Winds and/or jets are likely to happen in the presence of a large scale 
magnetic field which can have a key role in their the formation. To extend this, 
we may suggest that under-luminous AGN loose significant part of their energy 
to the formation of winds and/or jets. A search for a correlation between the 
occurrence of outflows or jets and low luminosity could be very informative. 
[Note that the ADAF model of Blandford and Begelman (1999) takes into ac- 
count outflows, but it requires the same assumption as earlier ADAF models 
that there is no magnetic field and no heating of electrons by field anihilation. 
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As we have argued, these assumptions are implausible.] 

The behavior of accretion disks in the presence of ordered and disordered 
magnetic fields was discussed. The interaction of an accretion disk with a ro- 
tating magnetized star was reviewed including the formation of magnetohy- 
drodynamic (MHD) outflows, magnetic braking of the star's rotation, and the 
propeller effect. The discussed problems of the MHD origin of jets in both in 
the hydromagnetic and Poynting flux regimes are currently under active study. 
Both regimes are being studied by different groups using axisymmetric and more 
recently 3D MHD simulations. Many simulation papers have been published on 
MHD outflows and jets, but the results often represent the time evolution for 
only a short interval (typically the period of rotation at the inner edge of the 
disk). These results are therefore strongly dependent on the initial conditions 
which are unknown. Furthermore, the boundary conditions used in the simula- 
tion codes can give false results, for example, coUimation of an MHD flow when 
there actually is no coUimation (see e.g. Ustyugova et ai, 1999). The difficulty 
of obtaining generally valid results from MHD simulations is even more chal- 
lenging in the case of relativistic MHD (Komissarov, 1999) including the Kerr 
metric near a rapidly rotating black hole (Koide et al. 2000). Energy losses 
connected with magnetohydrodynamic (MHD) outflows and jets from the ac- 
cretion disk may be very important and represent an alternative explanation of 
under-luminous AGN black holes. 
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Figure captions 



Fig.l 

Geometry of a thin accretion disk from Shapiro and Tcukolsky (1983). Here, 
M is the accretion rate, j is the total angular momentum flux, and = 
{l/r'^)d{r'^tr^) is the viscous force with — rjdrt/dr the main component 
of the viscous stress tensor and rj the dynamic viscosity coefficient. 

Fig.2 

The dependences of the optical depth tq on radius, r« = r/rg, for the case 
Mbh = 10^ Mq, a = 1.0 and different values of m. The thin solid, dot-triple 
dash, long dashed, heavy soHd, short dashed, dotted and dot-dashed curves 
correspond to m = 1.0,3.0,8.0,9.35,10.0,11.0,15.0, respectively. The upper 
curves correspond to the optically thick family, lower curves correspond to the 
optically thin family. 

Fig.3 

The specific angular momentum jin = iin as a function of the mass accretion 
rate M for different viscosity parameters a = 0.01 (squares), 0.1 (circles) and 
0.5 (triangles), corresponding to viscosity prescription (fz^). The solid dots 
represent models with the saddle-type inner singular points, whereas the empty 
dots correspond to the nodal-type ones, from Artcmova ct al. (2001). 

Fig.4 

The specific angular momentum jin = iin as a function of the mass accretion 
rate M for different viscosity parameters a — 0.01 (squares), 0.1 (circles) and 
0.5 (triangles), corresponding to viscosity prescription ( |11S| ). The solid dots 
represent models with the saddle-type inner singular points, whereas the empty 
dots correspond to the nodal-type ones, from Artemova et al. (2001). 

Fig.5 

The position of the inner singular points as a function of the mass accretion 
rate M for different viscosity parameters a — 0.01 (squares), 0.1 (circles) and 
0.5 (triangles), corresponding to viscosity prescription (fz^). The solid dots 
represent models with the saddle-type inner singular points, whereas the empty 
dots correspond to the nodal-type ones, from Artemova et al. (2001). 

Fig.6 

The position of the inner singular points as a function of the mass accretion 
rate M for different viscosity parameters a = 0.01 (squares), 0.1 (circles) and 
0.5 (triangles), corresponding to viscosity prescription ( |lig| ). The solid dots 
represent models with the saddle-type inner singular points, whereas the empty 
dots correspond to the nodal-type ones, from Artemova et al. (2001). 
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Fig.7 

Poloidal projection of the instantaneous magnetic field in an accretion flow. 



Fig.8 

Equatorial projection of the magnetic field given by equation (265) for an in- 
finitely conducting accretion flow. 

Fig.9 

Schematic drawing of the magnetic field configuration considered by LRBK95. 
The magnetosphere consists of an inner part, where the magnetic field lines are 
closed and an outer part where the field lines are open. The outer region where 
open field lines thread the disk gives rise to a MHD wind from the disk. 

Fig.lO 

This shows two poloidal field lines used in the derivation of equation (pOG). 



Fig. 11 

The figure shows the dependence of F^o and Fk = {GMrcr)^^^ on the pulsar 
period P assuming = for M = lO-'^^g/s and 4 x 10^^ g/s. The other 
parameters have been taken to be a = 0.1, D = 1, T,nax — 5, and M — Mq. 
Thus, each M = const, curve corresponds to different magnetic moments as 
jjL oc P^-^^ . Fao and Fk are measured in units of {GMtoY^"^ with Tq — lO^cm. 
This plot allows a comparison of the angular momentum influx to the star 
(MFk) in the case of outflows and spin- up {rto < Tcr) with the case of magnetic 
braking where the influx is MFoo{rto > Trc)- 

Fig.l2 

Geometry of disk accretion on to a rapidly rotating star with an aligned dipole 
magnetic field considered by LRBK99. Here, ri and r2 are the boundaries 
of the region; Ra is the effective radius where the outflow begins; is the 
star's rotation rate; ojk is the Keplerian rotation rate of the accretion disk; and 
Tcr = {GM /ujIY^^ is the corotation radius. The magnetic field in the vicinity of 
Ti has an essential time-dependence owing to the continual processes of stellar 
flux leaking outward into the disk, the resulting field loops being inflated by the 
differential rotation (LRBK95), and the reconnection between the open disk 
field and the closed stellar field loops. The dashed lines marked by the letters 
'n' indicate neutral surfaces along which reconnection occurs. 

Fig.l3 

Effective Alfven radius Ra (normalized by the nominal Alfven radius ta) as a 
function of the normalized corotation radius Tct/ta for a — 0.2 obtained from 



equation (33C). We have neglected the 5'^/'^ — (-Ra/''2)'^^^ terms compared with 
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unity and taken the value Tq — 0.234 from LBC. Only the part of the solid 
curve above the dashed line Ra — rcr is consistent with our assumptions. The 
dashed horizontal line rto indicates the "turnover radius" of the disk rotation 
curve or effective Alfven radius in the regime discussed by LRBK95 where Maccr 
falls onto the star and the star spins-up. For this line the turbulent magnetic 
diffusivity of the disk is taken to be aD = 0.1. The points abed and the 
associated vertical lines represent possible transitions between spin-down and 
spin- up of the pulsar as discussed in the text. 

Fig. 14 

Comparision of the observational and computational spectra of the cyclotron 
line. The solid curve is the observational results taken from Mihara et al. (1990), 
the dot curve is the approximation by the comptonized spectrum, and feature 
(17) with a = 7 • 10^"* sY:^, ct = 2 • 10"'' from Baushev and Bisnovatyi- 

V / cm ' cm ' ^ 

Kogan (1999). 
Fig.15 

Comparison of the observational and computational X-ray spectra of Her X-1. 
The solid curve is the observational results taken from McCray et al. (1982), the 
dot curve is the approximation by the comptonized spectrum, and the cyclotron 
feature (|5^) with = 0.9 KeV, Te = 8 KeV, = 14, a = 7 • lO"-* 
^ = 10"" S = 4 X 10^° Gs; from Baushev and Bisnovatyi-Kogan (1999). 

Fig.16 

Schematic structure of the accretion column near the magnetic pole of the neu- 
tron star (top), and its radiation spectrum (bottom); from Bisnovatyi-Kogan 
(1999a). 

Fig.l7 

Sketch of the magnetic field threading an accretion disk shown increase of the 
field owing to flux freezing in the accreting disk matter from Bisnovatyi-Kogan 
and Ruzmaikin (1976). 

Fig. 18 

Sketch of the electromagnetic outflows from the two sides of the disk owing to 
the Faraday unipolar dynamo action of a rotating magnetized disk (Lovelace 
1976). 

Fig.l9 

Sketch of an accretion disk threaded by a magnetic field for conditions which 
may lead to hydromagnetic jet formation. 

Fig.20 
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Simulations results for stationary MHD outflow obtained in a spherical coor- 
dinate system (Ustyugova et al. 1999). The solid lines represent the poloidal 
magnetic field, and the arrows the velocity vectors. The dashed lines in the 
left-hand plot represent the slow magnetosonic surface (the lowest dashed line), 
the Alfven surface (the middle line), and the fast magnetosonic surface (the top 
line) . In the right-hand panel, the dashed lines are surfaces of constant toroidal 
current density, while the lines on the outer boundary are the projections of the 
fast magnetosonic Mach cones. 

Fig.21 

Dependences of different velocities on distance s measured in units of rj from 

the disk along the second magnetic field line away from the axis in Figure 20 
(Ustyugova et al. 1999). This field line "starts" from the disk at r 6ri where 
it has an angle w 28° relative to the 2;— axis. The velocities are measured in 
units of GM/vi . Here, Vp is the poloidal velocity along the field line and 
is the poloidal velocity perpendicular to the field line. Also, VAp is the poloidal 
Alfven velocity, c/m is the fast magnetosonic velocity, and Vesc is the local escape 
velocity. 

Fig.22 

"Initial" dipole-likc vacuum magnetic field. In this and subsequent plots, r 
and z are measured in units of the radius ro of the O— point in the disk plane. 
The solid lines are the magnetic field linos for the case where the flux function 
on the disk surface is ^'(r, z = 0) = c£'r'^K/{a? + r^)''/^ with a = ro/\/2. In 
this and subsequent plots \E' is measured in units of \l/o = rgif/S'^/-^. Note that 
B^(0) = 6\/3*o/r-o « 10.4*0/^0- The dashed lines are the field lines for the case 
where the outer boundaries {Rmax = 10, Zmax = 10) are perfectly conducting; 
for this case *(r, z = 0) ^ *(r, 0)[1- (r/ro- 1)^/81] so that the O-point is stiU 
at ro and ^'o is unchanged. Because of axisymmetry and reflection symmetry 
about the z = Q plane, the field need be shown in only one quadrant. 

Fig.23 

Poloidal field lines for Poynting jet case for twist T = 1.84 rad. and {~H)max = 
l.lS^'o/ro with \1/ =const contours measured in units of \E'o which is the maxi- 
mum value of The outer boundaries at r = Rmax and z = Z^ax are perfectly 
conducting and correspond to an external plasma. This external plasma will ex- 
pand in response to the magnetic pressure of the jet so that Rmax and Zmax will 
increase with time. The initial poloidal magnetic field is shown by the dashed 
lines in Figure 22. The dashed contour is the scparatrix with the X-point in- 
dicated. Note that the radial width of the upgoing field lines along the axis is 
about one-half the width of downgoing field lines at the outer wall as required 
for equilibrium. 

Fig.24 
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Three dimensional view of two field lines originating from the disk at a; = 
±0.32 ro = 0.4*0) for the Poynting jet of Figure 23. Each field hne has a 
twist of « 8.22 rad. or about 1.31 rotations about the z— axis from its beginning 
at ri and end at r2- The z— axis is tilted towards the viewer by 30°. 

Fig.25 

Time evolution of dipole-like field threading the disk from the initial configura- 
tion t = Q (bottom panels) to the final quasi-stationary state t = 1.2tout, where 
tout is the rotation period of the disk at the outer radius Rmax of the simulation 
region. The left-hand panels show the poloidal field lines which are the same as 
^{r, z) =const lines; ^ is normalized by '^max and the spacing between lines is 
0.1. The middle panels show the poloidal velocity vectors Vp. The right-hand 
panels show the constant lines of — ri?^(r, z) > in units of "^max/fo and the 
spacing between lines is 0.1. 
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